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ABSTRACT 


The  theory  and  technique  of  estimating  power  spectral 
densities  by  Fraunhofer  diffraction  are  investigated  analyti¬ 
cally  and  experimentally.  A  technique  for  estimating  the  spec¬ 
tral  density  is  developed  by  which  the  estimate  is  independent 
of  changing  light  levels  in  the  diffraction  system.  It  is 
shown  that  the  estimated  spectrum  is  biased  and  in  the  case  of 
a  rectangular  aperture,  the  error  can  be  written  in  terms  of 
the  Hilbert  transforms  of  the  spectrum.  The  variability  of  the 
estimate  is  computed  and  the  ratio  of  the  variance  to  square 
mean  is  shown  to  approach  unity  for  large  apertures,  indicating 
the  need  for  smoothing  by  additional  filtering.  Reasonable 
bounds  on  the  bandwidth  of  the  scanning  filter  are  established 
theoretically.  The  experimental  results  are  found  consistent 
with  the  theoretical  evaluations . 

Techniques  for  reducing  the  zero-order  light  level  are  in¬ 
vestigated.  It  is  found  that  the  use  of  large,  apodised  aper¬ 
tures  are  useful  in  estimating  the  spectrum  at  low  frequencies. 
The  application  of  a  half-wave  phase  grating  for  measuring  the 
spectrum  of  grain-noise  at  DC  is  investigated  theoretically  and 
experimentally.  A  serious  limitation  to  measuring  the  grain- 
noise  at  DC  is  the  presence  of  residual  phase  noise. 

The  spectral  densities  of  several  common  photographic 
films  are  experimentally  determined  and  compared  with  models 
of  grain-noise  and  total  film-noise.  It  is  found  that  a  modi¬ 
fied  form  of  the  overlapping  circular  grain  model  of  film-grain 
noise  furnishes  an  excellent  fit  to  the  measured  spectra.  The 
range  of  Silberstein ' s  quantum  theory  of  photographic  exposure 
is  experimentally  explored.  An  analysis  of  grain-noise  in  the 
presence  of  signal  indicates  that  film-noise  cannot  be  regarded 
as  additive  or  multiplicative. 
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1. 


INTRODUCTION 


1.1  GENERAL 

The  research  described  in  this  thesis  has  been  concerned 
with  the  theoretical  and  experimental  aspects  of  power  density 
measurements  by  diffraction.  The  theoretical  design  considera¬ 
tions  and  the  experimental  techniques  that  evolved  from  the 
study  furnish  a  procedure  by  which  the  power  spectrum  of  trans¬ 
lucent  objects  can  be  obtained.  Because  of  the  outstanding  role 
played  in  modern  optical  systems  by  photographic  film,  it  was 
chosen  as  the  subject  of  the  experimental  investigation.  The 
motivation  for  measuring  film  noise  spectra  is  a  consequence 
of  the  following  recent  developments,  described  below. 

The  discovery  of  tne  laser  as  a  powerful  source  of  co¬ 
herent  light  has  enabled  the  construction  of  optical  signal 
processing  and  filtering  systems  which  offer  significant  ad¬ 
vantages  over  their  electronic  counterparts  in  certain  applica¬ 
tions.  These  advantages  stem  from  the  following  properties  of 
coherent  optical  systems: 

(1)  The  two  dimensions  available  in  optical  systems, 
offer  an  additional  degree  of  freedom  over  elec¬ 
tronic  systems. 

(2)  Coherent  optical  systems  can  be  made  to  generate 
successive  Fourier  transform  pairs  by  simply 
adding  lenses . 

(3)  Independent  control  over  the  amplitude  and  phase 
of  spatial  filters  is  easily  effected  and  im¬ 
plies  that  a  wide  variety  of  filter  functions 
can  be  synthesized  (Cutrona,  et  al.,  i960). 

In  these  systems,  the  filters  and  the  storage  media  are  almost 
always  photographic  transparencies. 


Another  activity  which  was  given  fresh  impetus  because 
of  the  development  of  the  laser  is  waveform  reconstruction  im¬ 
aging  or  holography  (Gabor  1948,  Gabor  1951 >  Stroke  1966).  In 
holography,  the  waveform  recording  media  are  generally  photo¬ 
graphic  transparencies . 

A  fundamental  limitation  in  the  performance  of  both  op¬ 
tical  signal  processing  and  wavefront  reconstruction  systems 
is  the  noise  introduced  by  the  photographic  film  (Goodman  1967 ) . 

The  noise  produced  by  photographic  film  results  from  two 
causes:  (1)  the  turbidity  of  the  emulsion  and  (2)  the  varia¬ 

tions  in  the  optical  path  length  in  the  cross  section  of  the 
transparency.  The  turbidity  is  a  consequence  of  the  suspension 
of  microscopic  silver  grains  in  the  emulsion.  These  grains  pro¬ 
duce  random  amplitude  fluctuations  in  light  and  this  effect  is 
known  as  grain  noise.  The  randomness  of  the  optical  path  length, 
which  results  from  variations  in  the  thickness  and  refractive  in¬ 
dex  of  the  emulsion,  produces  phase  errors  in  coherent  light. 

This  effect  is  similar  to  that  produced  by  random  imperfection 
in  lenses  (Rowe  I966). 

Both  of  these  effects  can  be  described  conveniently  by 
defining  a  complex  transmittance  of  a  transparency  as  the  ratio 
of  the  complex  amplitudes  of  the  transmitted  to  incident  light. 
The  complex  transmittance,  R(x,y),  can  be  written  as  the  pro¬ 
duct  of  a  real  amplitude  transmittance  t(x,y)  representing  the 
grain  noise  and  a  complex  exponential  phase  term  representing 
the  phase  noise. 

The  power  spectrum  of  a  photographic  transparency  fur¬ 
nishes  an  excellent  way  of  describing  the  noise  characteristics 
of  the  transparency.  In  optical  correlators,  where  the  signal- 
to-noise  ratio  is  approximately  inversely  proportional  to  the 
power  spectrum  of  the  transparency  noise,  knowledge  of  the  power 
spectrum  is  useful  for  obvious  reasons. 


The  possibility  of  spectral  density  measurements  by  dif¬ 
fraction  has  been  recognized  earlier  (Shack,  et  al.,  1961). 
The  advantages  offered  by  this  method  are  a  consequence  of  the 


following: 

(1) 

The  entire  three-dimensional  sample  is  involved 
simultaneously.  A  large  volume  is  thus  simul¬ 
taneously  analyzed  instead  of  a  single  trace. 

(2) 

The  spectrum  can  be  obtained  directly  without 
any  intermediate  operations  and  with  high  res¬ 
olution  . 

(3) 

The  spectrum  can  be  obtained  for  positive  and 
negative  frequencies.  This  is  useful  in  de¬ 
tecting  non-isotropic  noise. 

(4) 

The  diffraction  method  offers  the  possibility 
of  measuring  the  spectra  of  colloids,  turbid 
fluids,  etc.  for  which  the  application  of 
standard  techniques  is  difficult. 

1.2 


INVESTIGATIONS  AND  CONTRIBUTIONS 


The  conclusions  enumerated  in  this  section  relate  to  a 
coherent  optical  spectral  analyzer  which  in  simplest  form  con¬ 
sists  of  a  coherent  source,  two  lenses,  an  aperture,  a  slit  capable 
of  scanning  in  a  plane  and  the  object  to  be  analyzed  which  must 
be  able  to  transmit  light.  The  coherent  source  which  in  practice 
is  a  laser,  generates  a  spherical  beam  which  is  collimated 
(made  parallel)  by  a  lens  and  incident  thr  ugh  an  aperture 
on  the  object  located  in  the  front  focal  plane  of  the  second 
lens.  The  light  transmitted  by  the  object  is  passed  through  a 
second  lens  which,  according  to  Fraunhofer  diffraction  theory, 
produces  a  back  focal  plane  field  proportional  to  the  two- 
dimensional  Fourier  transform  of  the  front  focal  plane  field. 

The  scanning  slit,  which  is  a  spatial- frequency  bandpass  filter, 
enables  a  small  portion  of  the  focal  plane  field  to  impinge  on 
a  photomultiplier  tube  which  generates  a  current  proportional 
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to  the  power  of  the  beam  transmitted  by  the  slit.  The  total 

back  focal  plane  field  is  the  sum  of  the  field  produced  by  , 

the  diffraction  by  the  aperture  and  the  field  produced  by  the 

random  transmittance  of  the  object.  The  latter  generates  phase-  , 

front  perturbations  through  the  injection  of  random  phases  and 

amplitudes  resulting  from  inhomogeneities  in  the  body  of  the 

object . 

With  the  above  as  a  background,  the  theoretical  investi¬ 
gation  led  to  the  following  results: 

(l)  The  expected  value  of  the  field  generated  by  the 
random  transmittance  of  the  object  in  the  back  focal  plane  of 
the  second  lens  will  be  a  reasonable,  albeit,  biased,  estimate 
of  the  spectral  density  of  the  transmittance,  provided  that  the 
aperture  dimensions  are  large  compared  with  the  maximum  correla¬ 
tion  intervals  of  interest. 

(?)  There  are  two  contributions  to  the  bias,  both  re¬ 
sulting  from  the  use  of  a  finite  aperture.  Measurements  made 
with  an  aperture  of  transverse  dimension  D  cannot  furnish  infor¬ 
mation  about  correlation  intervals  greater  than  D.  This  pro- 

« 

duces  what  might  be  called  a  truncation  error  and  represents  a 
loss  of  information.  Even  when  D  is  large  enough  to  include 
all  significant  correlation  intervals,  there  is  still  a  contribu¬ 
tion  to  the  bias  because  the  optical  system  cannot  resolve  spec¬ 
tral  components  in  the  po//er  spectrum  to  better  than  D  1 .  This 
represents  a  resolution  limitation  and  contributes  to  the  bias. 

In  the  absence  of  a  truncation  error,  the  effects  of  finite  res¬ 
olution  can  be  negated  by  further  filtering  and  thus  do  not  con¬ 
stitute  a  loss  of  information. 

When  a  rectangular  aperture  is  used,  the  bias  error  can 
be  shown  to  be  directly  proportional  to  the  partial  derivatives 
of  the  Hilbert  transforms  of  the  true  power  spectrum  and  in-  * 

versely  proportional  to  the  dimensions  of  the  aperture  (Eq. 

(3.4-11)),  provided  that  the  truncation  error  is  zero.  If  the 
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truncation  error  is  not  zero,  the  above  is  still  true  provided 
that  the  true  spectrum  is  replaced  by  the  spectrum  of  the  trun¬ 
cated  process . 

(3)  Under  wide-sense  stationary  Gaussian  circumstances, 
the  variance  of  the  spectral  density  estimator  as  obtained  by 
diffraction  was  derived  for  a  rectangular  aperture.  It  was 
found  that  the  variance  depends  on  the  average  transmittance 
of  the  sample,  the  true  power  spectrum  of  the  process,  the 
transverse  fpatial  frequencies  in  the  back  focal  plane  and  the 
aperture  dimensions.  For  a  large  aperture,  the  ratio  of  the 
standard  deviation  to  the  mean  becomes  unity  everywhere  except 
at  the  origin.  This  result  shows  that  the  estimate  of  the  spec¬ 
tral  density  has  much  variability  which  must  be  reduced  by  means 
other  than  enlarging  the  aperture. 

Control  of  the  variability  of  the  power  spectrum  esti¬ 
mator  is  best  effected  by  the  scanning  filter  which  has  the 
effect  of  averaging  the  back  focal  plane  field  over  frequency. 
The  best  choice  of  slit  requires  a  compromise  between  stability 
and  resolution.  For  an  aperture  of  length  L  and  maximum  sig¬ 
nificant  correlation  interval  a  (it  is  assumed  that 

max 

L  o'  „)•  reasonable  bounds  on  the  width,  d,  of  the  slit 
max 

are  furnished  by 


where 

■\  =  wavelength  of  the  radiation 

f  =  focal  length  of  lens 

(4)  With  the  random  checkerboard  configuration  as  a 
basis  for  representing  the  granularity  of  photographic  film, 
it  was  shown  that  the  grain  noise  cannot  be  treated  as  uncor¬ 
related  additive  or  multiplicative  noise. 
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(5)  The  principal  interference  with  the  observation  of 
-  the  spectrum  comes  from  the  zero-order  light  which  represents 
the  diffraction  by  the  aperture.  Apodisation  by  aperture  shap¬ 
ing  may  be  used  to  reduce  the  sidelobe  level  of  the  zero-order 
for  selected  directions  in  the  back  focal  plane,  without  signi¬ 
ficantly  contributing  to  the  bias  of  the  measured  spectrum. 

The  peak  of  the  zero-order  occurs  at  the  origin  and  can  be 
many  orders  of  magnitude  above  the  spectral  density,  hence,  not 
allowing  its  observation  there.  A  null  in  the  zero-order  at  the 
origin  can  be  generated  by  using  odd-symmetric  aperture  functions, 
the  simplest  being  the  one-cycle  phase  grating.  The  effect  of 
using  this  grating  on  the  measurement  of  the  spectrum  is  to  de¬ 
crease  the  frequency  resolution  by  approximately  a  factor  of  two. 

The  experimental  phase  of  the  research  was  restricted 
to  measurements  on  photographic  film,  but  the  techniques 
which  were  used  could  be  applied  equally  well  to  other  translu¬ 
cent  objects. 

The  major  aims  of  this  phase  of  the  research  were  (1)  the 
development  of  techniques  by  which  reliable  and  accurate  spectral 
density  measurements  could  be  made,  (2)  the  verification  of  sev¬ 
eral  theoretical  results  and  (3)  the  measurement  of  the  proper¬ 
ties  of  film  noise  and  a  comparison  of  these  with  theoreti¬ 
cal  predictions  where  possible.  Specifically,  the  following 
experimental  results  were  obtained: 

(l)  A  technique  was  developed  by  which  the  estimate 
of  the  power  spectrum  of  film  grain  noise  could  be  obtained  di¬ 
rectly  and  continuously  with  high  reliability  down  to  frequen¬ 
cies  very  near  the  origin.  In  the  case  of  total  film  noise,  i.e., 
the  noise  resulting  when  random  amplitude  and  phase  effects  are  both 
present,  the  technique  enabled  the  observation  of  the  noise  to 
frequencies  as  low  as  5  cycles/mm.  The  technique  offered  sev¬ 
eral  advantages,  including  ease  of  calibration  directly  in  units 
of  mm"  and  furnishing  results  which  were  independent  both  of 
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system-scattered  light  and  changes  in  the  light  level  of  the 
system. 

(2)  The  dependence  of  the  spectral  estimate  on  system 
»  parameters  such  as  aperture  size,  shape  and  bandwidth  of  the 

scanning  slit  was  experimentally  obtained  and  the  results  were 
found  to  be  consistent  with  theoretical  predictions.  It  was 
found  that  the  general  shape  of  the  average  spectral  estimate 
was  relatively  independent  of  aperture  shape,  provided  that  the 
aperture  dimensions  were  large  compared  with  the  largest  signi¬ 
ficant  correlation  intervals  of  the  film  noise  process.  The 
size  and  shape  of  the  aperture  did  influence  the  variability  of 
the  estimate  in  a  way  which  was  theoretically  predicted,  namely 
that  the  variance  was  reduced  for  high  spatial  frequencies  and 
large  apertures. 

Although  increasing  the  aperture  size  reduced  the  vari¬ 
ability,  it  was  found  that  best  control  over  the  variability 
could  be  effected  by  increasing  the  bandwidth  of  the  scanning 
slit.  The  bandwidth  ultimately  used  corresponded  to  a  slit  size 
well  within  the  theoretical  bounds  given  earlier. 

,  (3)  Apodisation  by  aperture  shaping  was  found  to  be  a 

useful  technique  for  reducing  the  zero-order  light.  The  de¬ 
gree  of  sidelobe  suppression  was  not  as  great  as  that  theo¬ 
retically  predicted,  primarily  because  the  aperture  shape 
did  not  reduce  the  component  of  light  scattered  from  lenses, 
reflections,  etc. 

A  half-cycle  phase  grating  was  built  with  the  object  of 
reducing  the  zero-order  light  at  the  origin  in  order  to  measure 
the  spectrum  there.  The  performance  of  the  grating  corresponded 
closely  to  theory  and  the  feasibility  of  its  use  in  power  spec¬ 
tral  measurements  at  the  origin  was  firmly  established. 

1  (4)  Good  agreement  was  found  between  the  measured  spec¬ 

trum  and  the  theoretical  spectrum  of  the  3-parameter  overlapping 
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circular  grain  (OCG)  model  which  includes  the  effect  of  the 
film  gelatin.  The  2-parameter  OCG  model  and  the  random  checker¬ 
board  (RC)  models  furnished  power  spectra  which  were  somewhat 
higher  than  the  measured  spectrum,  hence  enabling  S/N  estima¬ 
tions  based  on  these  models  to  be  on  the  "safe"  side.  A  serious 
drawback  of  the  RC  model  is  that  it  does  not  have  isotropic 
statistics  and  hence,  the  goodness  of  fit  depends  on  the  orien¬ 
tation  of  the  model  "grains." 

The  measurements  enabled  the  determination  of  the  param¬ 
eters  of  these  models.  The  knowledge  of  these  parameters  proved 
useful  in  determining  additional  information  about  the  relation 
of  film  granularity  structure  to  average  transmittance  and  the 
effect  of  developer  on  grain  size. 

The  variation  of  the  grain  noise  spectra  with  average 
amplitude  transmittance  was  experimentally  determined  for  two 
film  types.  In  each  case,  the  maximum  was  observed  in  the  spec¬ 
trum  when  the  average  transmittance  was  near  0.6.  The  existence 
and  location  of  such  a  maximum  are  predicted  by  both  tho  RC  and 
OCG  models . 

Examination  of  the  spectrum  of  non-isotropic  film  grain 
noise  was  found  to  be  a  useful  technique  for  detecting  direc¬ 
tional  anomalities  such  as  granularity  streaking. 

(4)  Measurements  indicated  that  the  total  noise  of  film 
was  a  low-frequency  phenomena.  Typically,  the  total  noise  de¬ 
creased  to  one-half  its  maximum  value  in  a  frequency  span  one- 
fifth  that  of  grain  noise.  At  the  origin,  the  total  noise  was 
much  larger  than  the  grain  noise  -  an  order  of  magnitude  not 
being  unusual. 

The  magnitude  of  the  complex  transmittance  was  found  to 
vary  linearly  with  the  average  amplitude  transmittance.  This 
result  furnished  the  variance  of  the  random  phase  within  the 
framework  of  a  model  which  assumes  independence  between  ampli¬ 
tude  and  phase  and  Gaussian  statistics  for  the  latter. 


-8- 


(5)  Spectral  density  measurements  can  be  used  to  de- 

« 

termine  the  frequency  response  of  optical  systems  because  the 
transition  from  object  to  image  in  optics  is  equivalent  to  the 
action  of  a  linear  filter.  In  addition,  imaging  and  magnifica¬ 
tion  with  high-resolution  optics  can  be  used  to  obtain  the  spec¬ 
trum  over  extended  frequencies  without  generating  unsuitably 
large  diffraction  angles.  Experiments  designed  to  demonstrate 
the  feasibility  of  these  techniques  were  performed  with  a  meas¬ 
ure  of  success . 

1.3  REVIEW  OF  FILM-NOISE  MEASUREMENTS 

Photographic  granularity  of  grain  noise  has  been  a  sub¬ 
ject  of  intense  investigation  for  at  least  fifty  years.  Con¬ 
tributions  have  come  from  scientific  workers  in  many  fields, 
including  chemistry,  physics,  mathematics  and  engineering.  A 
great  and  diverse  fund  of  knowledge  is  scattered  in  the  liter¬ 
ature  in  several  languages  and  many  journals  .  The  investigator 
can  at  best  only  hope  to  refer  to  those  works  which  occupy,  in 

* 

his  opinion,  salient  positions  in  the  continuing  investigation 
of  this  subject. 

e 

One  of  the  earliest  attempts  to  view  photographic  granu¬ 
larity  as  a  random  process  was  the  work  by  Silberstein  in  1922. 
Silberstein  viewed  the  grain-development  process  from  a  quantum- 
mechanical  viewpoint.  According  to  his  theory,  the  probability 
that  a  silver-halide  grain  becomes  exposed  (i.e.,  is  ultimately 
reduced  to  metallic  silver)  is  a  function  of  two  numbers:  (l)  the 
probability  e  that  a  photon  incident  upon  a  grain  is  absorbed 
by  the  grain,  and  (2)  the  quantum  threshold  m  which  is  the 
number  of  photons  required  to  make  a  grain  developable  (Silber¬ 
stein  1922).  Silberstein  later  amended  his  theory  to  include, 
in  any  emulsion,  grains  of  varying  quantum  sensitivities  (Silber- 

>  stein  1941).  The  Silberstein  theory  was  experimentally  justi¬ 

fied  by  Trivelli  and  Righter  (Trivelli  dc  Righter  1922). 
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Attempts  to  measure  photographic  granularity  through  micro¬ 
densitometer  traces  date  back  to  1932  [Mees  and  James (1 )* 1966 1 . 
Attempts  at  evaluating  these  traces  did  not  occur  until  the 
following  year,  when  it  was  observed  that  the  random  optical 
density  of  uniformly  exposed  film  was  in  good  agreement  with  the 
Gaussian  probability  law  [Mees  and  James (2 ) , (3 ) , (4 )  19661.  How¬ 
ever,  it  was  expected  and  ultimately  verified  that  the  variance  of 
the  optical  density,  c^2  ,  varies  with  the  scanning  aperture. 
Bricout  found  a  variation  approximately  as  d0,e,  where  d  is 
the  aperture  diameter.  These  results  limited  the  usefulness 
of  using  crD2  alone  as  a  measure  of  granularity. 

In  1935 s  a  theoretical  study  based  on  the  examination 
of  three  sets  of  assumptions  about  the  nature  of  photographic 
transparencies  led  to  a  new  measure  of  granularity  called  the 
Selwyn  granularity  coefficient  or,  simply,  the  Selwyn  granu¬ 
larity.  The  assumptions  leading  to  this  new  measure  were  the 
following;  (l)  the  number  of  grains  within  a  scanning  aperture 
is  a  random  sample  from  a  binomial  distribution  and  the  optical 
density  is  strictly  proportional  to  the  number  of  grains  in  the 
aperture;  (2)  the  emulsion  can  be  described  as  an  aggregate  of 
elementary  cubes  that  consist  of  opaque  silver  grains  or  trans¬ 
parent  gelatin;  and  (3)  the  density  distribution  involves  only 
one  parameter  in  addition  to  the  area  of  the  scanning  spot, 
densities  of  adjacent  neighboring  areas  are  uncorrelated  and 
fluctuations  of  density  are  small .  Consideration  of  these  as¬ 
sumptions  led  to  the  conclusion  that  the  observed  variance  is 
inversely  proportional  to  the  square  root  of  the  scanning  area 
and  therefore,  that  the  product  of  the  two  must  be  a  constant 
[Mees  ana  James  (5)  1966],  The  Selwyn  granularity  is  defined 


by 


Number  in  parenthesis  refers  to  reference  cited  in  the  work 
by  C.E.K.  Mees  and  T.  H.  James.  See  2  references  at  end. 


where  observed  variance,  and  a  -  area  of  the  scanning 

aperture . 

Experimental  attempts  to  verify  the  above  led  to  mixed 
results.  Departures  from  the  square-root  law  were  observed  by 
several  workers  (Bricout,  van  Kreveld,  Debot )  and  in  one  case 
fMees  and  James  (6)  I966],  it  was  found  that  the  square-root 
law  fitted  experimental  observations  quite  well  if  observations 
were  made  on  the  transmittance  T  instead  of  the  density  D 
(D  1o9  10  1/T)  . 

Careful  studies  on  granularity  revealed  that  the  Selwyn 
measure  for  granularity  was  a  function  of  optical  density  and 
failed  for  both  very  small  and  very  large  apertures  (Jones  <fc 
Higgins  1946).  Furthermore,  it  was  found  that  the  Selwyn  co¬ 
efficient  did  not  always  produce  results  that  corresponded  to 
subjective  visual  evaluation.  To  eliminate  the  latter,  the 
Selwyn  measure  was  modified  to  give  way  to  a  syzygetic  Selwyn 
granularity  coefficient.  This  measure  attempts  to  duplicate 
the  effect  of  granularity  on  retinal  receptors  (Jones  dc  Higgins 
1947).  To  obtain  the  syzygetic  coefficient,  granularity  is 
measured  in  terms  of  the  magnification  for  which  the  average 
density  difference  between  two  microdensitometer  apertures  at¬ 
tain  a  certain  value. 

With  the  advent  of  more  rapid  and  accurate  instruments, 
the  Selwyn  granularity  was  measu  red  fo'-  a  number  of  samples 
as  functions  of  both  the  density  and  the  square  root  of  scan¬ 
ning  area  (Higgins  &  Stultz  1959).  The  results  may  be  summar¬ 
ized  as  follows:  in  general,  the  granularity  increases  with 
density  (or  decreases  with  transmittance).  For  an  aperture  in 
the  range  between  7  and  200  microns,  the  Selwyn  coefficient  is 
fairly  constant.  Above  200  microns,  S  increases  quite  sharply 
with  the  square  root  of  the  aperture.  Results  for  apertures 
smaller  than  7  microns  are  not  available. 
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Granularity  has  been  linked  to  entropy.  If  the  en¬ 
tropy,  H,  of  a  transparency  is  defined  by 

H  =  -Tlog2T  -  (l-T)log  (1-T) 

where  T  is  the  intensity  transmittance,  then  it  has  been 
shown  that  the  function  H  vs  T  closely  follows  the  subjec¬ 
tive  impressions  of  "graininess"  as  a  function  of  transmittance 
(O'Neill  1958). 

During  the  nineteen-fifties,  the  mathematical  tools 
of  communication  theory  were  applied  to  grain-noise  analysis. 

The  usefulness  of  autocorrelation  methods  in  studying  granu¬ 
larity  became  apparent  (Fellgett  1953)*  The  relation  of 
Selwyn' s  granularity  measure  to  the  characteristics  of  the 
autocorrelation  function  was  investigated  by  H .  J.  Zweig 
(Zweig  1956),  who  computed  autocorrelation  functions  from 
microdensitometer  traces  and  found  that  the  size  of  the  scan¬ 
ning  aperture  had  a  profound  influence  on  the  shape  of  the  auto¬ 
correlation  curve.  By  varying  the  diameter  of  the  scanning 
aperture  from  6  microns  to  50  microns,  Zweig  obtained  for  the 
same  sample  very  peaked  autocorrelation  functions  with  signi¬ 
ficant  autocorrelation  intervals  never  larger  than  8  microns 
in  the  case  of  the  6  micron  aperture  to  almost  flat  auto¬ 
correlation  curves  with  significant  intervals  extending  to 
20  microns  in  the  case  of  the  50-micron  aperture. 

Reciprocal  relationships  between  the  autocorrelation 
function  and  the  Selwyn  granularity  were  derived  (Marriage  & 
Pitts  1956).  In  a  widely  quoted  paper,  R.  C.  Jones  argued 
that  the  only  fully  adequate  way  to  describe  the  granularity 
of  photographic  materials  is  by  means  of  the  noise  power  spec¬ 
trum  (R.  C.  Jones  1955)* 

The  power  spectrum  of  the  film  grain  noise  goes  be¬ 
yond  any  previous  measure  in  that  it  leads  to  the  solution  of 


many  of  the  signal- to-noise  problems  that  arise  in  connection 
with  the  detection  of  signals  on  photographic  transparencies. 

The  power  spectrum  can  easily  be  related  to  the  Selwyn  coef¬ 
ficients  . 

A  very  common  technique  for  measuring  the  power  spec¬ 
trum  of  the  grain  noise  consists  of  scanning  a  uniformly  ex¬ 
posed  sample  of  film  with  a  microdensitometer  and  passing  the 
electrical  output  of  the  densitometer  through  a  narrow  band¬ 
pass  filter.  The  spectrum  may  then  be  explored  by  varying 
the  center  frequency  of  the  pass  band.  This  technique  has 
been  applied  to  obtain  the  spectra  of  various  transparencies 
(Jones  1955,  Doerner  1961,  Steel  1966). 

The  method  just  described  has  several  disadvantages. 

Since  this  technique  really  measures  the  power  spectrum  of 
the  electrical  noise  signal  produced  by  the  microdensitometer, 
a  relation  between  the  electrical  noise  spectrum  and  the  grain- 
noise  spectrum  must  be  established.  The  general  relation  as 
computed  by  R .  C.  Jones  (Jones  1955)  depends  upon  the  relative 
velocity  between  the  film  sample  and  the  densitometer,  the  shape 
and  size  of  the  aperture,  the  frequency  response  of  the  densi¬ 
tometer,  and  the  noise  spectrum.  It  is  not  generally  possible 
to  determine  the  grain-noise  spectrum  from  the  electrical  spec¬ 
trum  except  for  certain  densitometer  apertures.  The  spectra 
obtained  by  this  method  tend  to  show  a  characteristic  droop 
which  reflects  the  filtering  action  of  the  aperture  and  not 
any  property  of  the  random  noise  process.  Also,  for  reasons 
that  probably  pertain  to  the  electrical  system,  very  low  spatial 
frequencies  (/l  line/mm)  are  usually  not  investigated  with  the 
above  method . 

With  the  advent  of  lasers,  high-intensity  coherent 
optical  systems  became  useful  as  signal  processing  and  filter¬ 
ing  systems.  The  two-dimensional  Fourier  transforming  property 
of  lenses  could  be  exploited.  Successive  linear  operations  of 


an  integral- transform  nature  could  be  implemented  by  simply 
adding  additional  lenses  to  a  system.  In  a  basic  paper, 

Cutrona  showed  that  most  of  the  filtering  and  signal  process¬ 
ing  operations  performed  by  electronic  equipment  could  be  dupli¬ 
cated  by  using  elements  of  coherent  optical  systems  (Cutrona 
i960).  Electro-optical  spectrum  analyzers  using  the  Debye-Sears 
effect  were  constructed  and  proved  useful  in  processing  and  spec¬ 
trum  analyzing  electrical  signals  (Lambert  1965). 

With  respect  to  film-grain  noise,  it  was  discovered 
that  when  a  photographic  transparency  was  placed  in  the  col¬ 
limated  region  of  a  coherent  optical  system,  scattered  light 
appeared  in  the  back  focal  plane  of  the  Fourier  transforming 
lens.  Explained  in  terms  of  Fraunhofer  diffraction  theory, 
the  scattered  light  was  adjudged  to  be  the  spectral  distribu¬ 
tion  of  granularity  (Shack,  et  al .  1961,  Leith  1962,  Thiry  1963). 

It  thus  appeared  that  the  spectrum  of  grain  noise  could  be  meas¬ 
ured  by  diffraction  techniques.  Attempts  to  measure  the  scat¬ 
tered  light  intensity  proved  to  be  a  difficult  task  chiefly  be¬ 
cause  the  diffracted  source  light  (often  called  the  zero-order 
light),  was  several  orders  of  magnitude  more  intense  than  the 
light  of  interest. 

One  technique,  which  in  principle  could  reduce  the 
sidelobe  level  of  the  zero-order  light,  would  be  to  weight 
the  aperture  with  an  absorbing  pupil  function  (jacquinot  1966). 
However,  the  fabrication  of  accurately  controlled  pupil  func¬ 
tions  is  a  difficult  task, and  such  pupil  functions  would  them¬ 
selves  be  granular  transparencies.  An  indirect  technique  in¬ 
volving  a  combination  of  diffraction,  spatial  filtering  and  im¬ 
aging  can  be  used  to  measure  the  light  scattered  by  the  grains 
in  a  transparency  (Leith  I962).  This  method  involves  a  three- 
step  operation  in  which  the  combined  zero-order  light  and  scattered 
light  is  spatially  filtered  in  the  frequency  plane  and  then  re-imaged 
on  a  fixed  photocell.  According  to  the  author,  if  the  bandpass 
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region  of  the  spatial  filter  is  appropriately  chosen,  the 
zero  order  light  can  be  displaced  from  the  region  where  the 
scattered  light  is  being  measured.  Unfortunately,  this  tech¬ 
nique  involves  such  low  light  levels  that  a  reasonably  wide 
spatial  filter  must  be  used.  A  wide  spatial  filter  implies 
observing  a  smoothed  version  of  the  scattered  light  variation 
and  the  fine  structure  in  the  scattered  light  is  consequently 
lost . 

Attempts  to  measure  grain-scattered  light  by  observ¬ 
ing  the  perturbation  in  the  diffraction  pattern  of  a  very  small 
circular  aperture  were  made  using  incoherent  light  sources  (Thiry 
1963).  A  combination  of  diffraction  and  electronics  techniques 
enabled  Burckhardt  to  measure  the  scattered  light  from  holo¬ 
graphic  plates  down  to  200  lines/millimeter.  Data  were  not 
obtained  for  spatial  frequencies  lower  than  200  lines/mm 
(Burckhardt  19^7  )  • 

In  addition  to  the  noise  resulting  from  the  granularity 
of  a  transparency,  a  coherent  optical  system  will  also  sense 
noise  resulting  from  random  phase  misalignments  in  the  complex 
light  amplitude  emerging  from  such  a  transparency.  The  random 
phases  result  from  variations  in  both  the  thickness  of  the 
transparency  and  the  refractive  index  of  the  medium.  It  has 
been  shown  that  in  order  to  stay  within  the  Rayleigh  limit 
(maximum  optical  path  difference  across  aperture  <(  \/k )  the 
transparency  must  be  inserted  in  a  liquid  gate  whose  refractive 
index  is  equal  to  the  refractive  index  of  the  transparency 
(Ingalls  i960 )  . 

Several  theoretical  models  exist  for  describing  the 
grain  distribution  in  an  emulsion  (O'Neill  1963,  Picinbono 
1955,  Savelli  1958)*  Some  of  these  models  have  been  used  to 
predict  signal- to-noise  ratios  in  wavefront  reconstruction  im¬ 
aging  (Goodman  1967).  Two  commonly  used  theoretical  models  are 
the  overlapping  circular  grain  model  and  the  random-checkerboard 
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model.  The  former  is  the  more  realistic  of  the  two  and  models 
the  grains  as  opaque  disks  which  may  overlap.  If  a  point  is 
covered  by  a  disk,  the  transmittance  of  the  point  is  zero  and 
if  the  point  is  not  covered  by  a  disk,  the  transmittance  is 
unity.  The  random-checkerboard  model  is  the  two-dimensional 
extension  of  the  random  binary  waveform  in  electrical  engin¬ 
eering.  Comparisons  between  the  theoretical  models  and  exper¬ 
imental  data  have  led  to  mixed  results  (Jones  1955)* 
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2 .  BASIC  RELATIONS 


In  this  chapter  the  basic  equations  relating  to  the 
measurement  of  power  spectra  by  the  diffraction  of  coherent 
light  are  derived  for  ordinary  and  astigmatic  optical  systems. 
Ordinary  optical  systems  perform  two-dimensional  operations 
while  the  astigmatic  systems  considered  here  are  one-dimensional. 
Astigmatic  optical  systems  are  simulated  in  practice  by  using 
cylindrical  instead  of  spherical  lenses  and  scanning  along  a 
specific  axis  instead  of  in  a  plane.  Both  ordinary  and  astig¬ 
matic  optical  systems  are  assumed  to  be  free  from  geometrical 
aberrations . 

2.1  COMPLEX  TRANSMITTANCES 

A  convention  that  appears  in  the  literature  on  coherent 
optics  (Born  and  Wolf  1964)  and  is  used  here  is  the  following: 

A  wave  of  the  form 


e  ( x  ,y  ,  z  ,  t )  A  (x  ,y  )  cos  [£  z-cut  f  7(x,y)j  (2.1-1) 


is  denoted  as  the  complex  function  E(x,y)e 


j (pz-^t ) 


where 


E(x,y)  A(x,y)  e-^X’^ 


(2.1-2) 


E(x,y)  is  commonly  called  the  complex  amplitude  of  a  coherent 
light  wave.  The  other  symbols  have  the  following  meanings: 

A(x,y)  real  amplitude,  a  function  of  (x,y) 

cn  temporal  radian  frequency  of  the  wave 
fi  the  propagation  constant  -  2 tt/A 
A  the  wavelength  of  the  wave  in  the  medium 
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-  the  axis  of  propagation,  also  the  position 
of  the  wavefront 
t  -  time  parameter 

7  (x,y)  the  relative  phase,  a  function  of  (x,y) 
j  =  v  -1 

The  principal  motivation  for  adopting  the  representation 
of  Eq.  (2.1-1)  of  a  wave  is  that  (i)  all  the  significant  fea¬ 
tures  of  the  optical  systems  considered  here  are  time  invariant, 
(ii)  the  notation  is  very  useful  when  defining  the  complex 
transmittance  of  an  object,  and  (iii)  all  physical  sensors 
sense  the  power  or  energy  in  a  wave  and  these  quantities  are 
proportional  to  the  time  average  of  e2  or  E (x ,y )E* (x,y ) 
where  *  denotes  conjugate.  In  the  following  work,  the  in¬ 
tensity  associated  with  a  wave  of  complex  amplitude  E(x,y)  is 
defined  as 

I (x,y )  -  I  E(x,y)  i  2  (2.1-3) 

A  transparent  or  semi-transparent  object  illuminated  by 
an  incident  coherent  light  wave  as  shown  in  Fig.  2.1-1  will 
generally  alter  the  amplitude  and  phase  of  the  light  as  the 
latter  passes  through  the  object.  The  effect  of  the  object  on 
the  transmission  of  light  is  conveniently  described  by  defining 
a  complex  transmittance  of  the  object  by  the  ratio  of  the  trans¬ 
mitted  to  incident  light.  The  complex  transmittance,  denoted 
as  R(x,y)  is  then  given  by 

ET(x,y ) 

R(x’y)  ^  etrtt  (2-1'4) 

where  EI(x,y)  and  ET(x,y)  denote  incident  and  transmitted 
light  respectively.  When  arg  R(x,y)  -  0,  the  transmittance 
is  purely  an  amplitude  transmittance  and  is  denoted  by  t(x,y). 
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More  generally,  both  amplitude  and  phase  are  affected  and  it 
is  convenient  to  write 

R(x,y)  =-  t (x,y )  ej^(x’y)  (2.1-5) 

Observe  that  t(x,y)  is  a  real  number  such  that 
0  {  t(x,y)  (  1  and  V'(x,y)  is  a  real  number  resulting  from 
variations  in  the  optical  path  length  across  the  object.  If 
the  thickness  of  the  object  at  (x,y)  is  rr(x,y)  then  the 
phase  retardation  V^(x,y)  at  (x,y)  when  the  object  is  assumed 
immersed  in  air  is  'given  by 

( x , y  )  =  (2ttA)  (n-1)  cr(x,y)  (2.1-6) 

where  n  is  the  index  of  refraction  of  the  object. 

The  values  of  R(x,y)  may  be  represented  by  points  on 
or  within  a  circle  of  unit  radius  in  the  complex  plane,  as  shown 
in  Fig.  2.1-2.  When  arg  R(x,y)  =  0,  implying  a  uniform  phase 
retardation  over  the  entire  object,  the  domain  of  R(x,y)  is 
the  segment  Oa.  When  the  object  has  a  uniform  amplitude  trans¬ 
mittance  with  varying  phase,  the  range  of  R(x,y)  consists  of 
points  lying  on  a  circle  concentric  and  within  the  unit  circle. 

In  many  physical  situations,  the  complex  transmittance 
of  an  object  is  a  randomly  varying  quantity,  the  randomness 
being  inherent  in  t(x,y),  ^(x, y)  or  both.  To  cite  some 

examples,  consider  the  situation  that  arises  when  uniformly 
exposed  photographic  film  is  inserted  into  a  liquid  gate  whose 
refractive  index  is  equal  to  that  of  the  film's;  the  optical 
path  length  at  all  points  is  the  same  and  the  relative  phase 
is  zero.  The  granularity  of  the  film,  however,  produces  ran¬ 
dom  amplitude  fluctuations  in  the  transmitted  light.  In  this 
case  R(x,y)  =  t(x,y),  and  the  latter  can  be  considered  as  a 
real  random  variable  for  any  point  (x,y).  On  the  other  hand. 
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FIG.  2.1-2  THE  RANGE  OF  THE  COMPLEX  TRANSMITTANCE  R(x,y) 
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a  thin  lens  with  random  surface  imperfections  produces  only  a 
random  phase  shift  ip(x,y)  while  maintaining  t(x,y)  at  es¬ 
sentially  unity(Rowe  1966).  Finally,  photographic  film  which 
is  not  in  a  liquid  gate  contains  random  variations  in  both 
t(x,y)  and  ^(x,y).  The  random  variables,  R(x,y)  and  t(x,y) 
do  not  have  zero  mean.  For  subsequent  work,  it  is  convenient 
to  define  the  auxiliary  random  variables  r(x,y)  and  T(x,y) 
as  follows 


r (x ,y )  -  R(x,y)  -  Rq 

/  \  ,  .  (2.1-7) 

t (x,y )  -  t (x,y )  -  tQ 

where  Rq  -  /R(x,y)y  and  t(x,y)  =  7tQ)>  and  7  )>  denotes 

expectation.  Clearly  /r(x,y)N  -  7x(x,y)  ).  =■  0. 

2.2  SPATIAL  RANDOM  PROCESSES 

Any  point  in  the  xy  plane  can  be  considered  as  a  vector 
with  initial  point  at  the  origin  and  terminus  at  the  point. 

Such  a  point  can  be  defined  by  either  the  pair  (x,y)  or 

^  A  A  A  A 

the  vector  s  =  xi  +  yi  where  i  and  i  are  orthogonal 

x  y  x  y 

unit  vectors  parallel  to  the  coordinate  axes.  When  s  is  any 
given  vector,  R(s)  is  the  random  variable  describing  the  com¬ 
plex  transmittance  at  s.  Note  that  r(s)-  is  a  scalar  function 
of  the  vector  s.  When  s  is  used  to  index  R,  a  collection 
of  random  variables  is  formed,  one  for  each  point  in  the  plane. 
This  collection  of  random  variables  constitutes  the  random  pro¬ 
cess  | R 1  s  ) } .  A  similar  argument  holds  for  |t(s)|  and  j '/  ( s  )  | . 

In  all  of  the  present  work,  the  complex  process  |r(s)| 
and  the  real  processes  |t(s)|  and  |^(s)|  are  assumed  to 


be  wide  sense 

stationary. 

This  requires  that 

(1) 

SR(s)  >  = 

Rq  independent  of  s 

(2.2-1) 

(2) 

<r(s  +  7 

)  R*(S)  >  =  *R(i  +  i, s) 

=  «R(i) 

(2.2-2) 
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i.e.,  the  autocorrelation  function  of  |r(s)  J  is  only  a  func¬ 
tion  of  the  difference  of  the  two  vectors,  and 

(  -  )  /  1  R  ( s’ )  I  s>  oo  (2.2-3) 


The  covariance  of  ^  R T s  'l  ^  is  the  autocorrelation  function  of 
jr(s  |  and  is  given  by 


(  ( 


«Ra>  - 


I  2 


(2.2-4) 


From  the  definition  of  0  ft),  it  follows  that 

_s  K 

•  (  •  '  *(-()•  Hence,  in  rectangular  coordinates,  there 

K  R 

results 


0Rfrt,e)  -  0R*  (- a,-F)  (2.2-F) 

from  which  it  follows  that 

R(>  )  =-  4r  '  ( -a  ,  \i  ) 

R  -  '  ^  '■  :R  (  *  ' 

In  certain  physical  situations,  it  is  reasonable  to  as¬ 
sume  that  the  covariance  function  0  (t )  is  a  function  only 
of  the  magnitude  of  t,  i.e.,  the  covariance  function  is  a 
function  of  a  scalar.  In  this  case  the  statistics  are  inde¬ 
pendent  of  the  direction  of  the  coordinate  axes  and  the  covar¬ 
iance  depends  only  on  the  distance  from  the  origin  of  the 
rt.p  coordinate  axes.  Hence,  one  may  write  for  this  special 
process , 

0rf-v.p. )  -  .'r  +  ((  '  where  ,  (c,p  i  h’’  ) 

(2.2-3) 

t  The  symbol  will  be  used  throughout  to  denote  correlation 

or  covariance  functions.  A  subscript  will  be  added  only 
when  there  is  a  need  to  differentiate  among  various  processes. 


■  c  y 


and  since 


0r(a,£)  -  0r(  a, -3)  =  </>/(«, 6)  (2.2-7) 

the  covariance  function  0r(o,p)  is  a  real  function.  When 
the  covariance  may  be  written  as  a  function  of  the  single  argu¬ 
ment  p  as  in  Eq.  (2.2-6),  the  random  process  is  said  to  be 
isotropic .  The  covariance  of  an  isotropic  process  possesses 
circular  symmetry  in  the  sense  that  the  loci  of  constant  <M|  ) 
are  circles  about  the  origin. 

When  the  spatial  random  process  is  the  noise  associated 
with  uniformly  exposed  film,  it  is  reasonable  to  assume  that 
the  process  is  stationary.  The  statistical  properties  of  the 
granularity  and  surface  corrugations  are  the  same  from  one  por¬ 
tion  of  the  film  to  another  (Jones  1955). 

There  is  far  less  justification  for  assuming  on  an  a- 
priori  basis  that  film  noise  is  isotropic.  This  is  equivalent 
to  assuming  that  there  is  no  preferred  direction  on  the  surface 
of  the  film  -  that  film  noise  is  not  streaked  in  any  particular 
direction.  That  the  isotropic  noise  condition  may  indeed  be  un¬ 
justified  can  be  seen  from  the  fact  that  when  emulsion  is  spread 
on  the  film  base,  there  is  a  likelihood  that  the  granularity  ac¬ 
quires  a  directional  character,  related  to  the  direction  in 
which  the  emulsion  is  spread  on  the  film  base. 

The  questions  concerning  the  isotropy  of  film  noise  and 
for  that  matter,  any  other  spatial  noise  process  may  be  resolved 
by  examination  of  the  power  spectrum  of  the  process.  Before 
dealing  with  this  problem,  it  is  noteworthy  that  in  the  case  of 
film  noise,  the  process  |t(x,y)|  is  closely  associated  with 
the  granularity  of  the  film  (the  dichotomous  state  of  grains), 
while  the  noise  process  \ip{x,y)^  is  primarily  a  consequence 
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of  the  surface  irregularities  of  the  emulsion, which  random¬ 
ize  the  optical  path  length  and  produce  random  phase  errors 
in  coherent  light.  More  will  be  said  about  film  noise  in 
Chapter  4 . 

2.3  POWER  SPECTRA 

The  two-dimensional  spectral  density  of  a  stationary 
random  process  is  the  Fourier  transform  of  its  covariance 
V-(o  .r').  Thus  if  W(u,v)  denotes  the  spectral  density. 


W(u,v) 


u.  ,--f  00 

J 

—  oo 


0  (a, a ) 


-  j27i(ua+v3  ) 

e 


dad(3 


(2.3-1) 


where  u  and  v  are  the  transverse  frequencies  of  the  " fre¬ 
quency  domain." 

The  power  spectrum  is  a  real  function  of  the  arguments 
u,v.  This  is  easily  seen  by  considering  the  conjugate  of 
Eq.  (2.3-1)  and  using  the  fact  that  0*(a,3)  =  0(- a,-p).  W(-u,v) 
is  given  by 

,  ,  -  j2)T  (-ua+vp ) 

W ( — u , v )  =  /+0C/+00  0( a,P)  e 


+00 


—  00  —00 


-  j27r(ua+v^ ) 

0(-a,0)  e  dadp 


(2.3-2) 


Since  0(-a  ,0)  is  not  generally  equal  to  0(a,(3),  it  may  be 
concluded  from  the  uniqueness  property  of  the  Fourier  transform 
that  W(-u,v)  is  not  equal  to  W(u,v).  An  identical  argu¬ 
ment  can  be  applied  to  show  that  W(u,v)  is  not  equal  to 
W(u,-v). 
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An  exception  occurs  in  the  isotropic  case, for  which 

0(a,0)  =0  (7a2  +  32  )  r  0(f  ).  Let 


a 

=  p 

cos 

s'* 

6 

=  1 

sin 

0 

u 

=  w 

cos 

c 

V 

=  w 

sin 

c 

Equation  (2.3-1)  may  then  be  written  as 

W(w  cos  C,w  sin  £ )  =  2rr  /°°  0(|  )  J  (2ir  wp  )p  dp  -  W(w) 

0  ° 

(2.3-4) 

Clearly  in  the  isotropic  case,  W(u,v)  is  an  even  function 
with  respect  to  both  arguments . 

Equation  (2.3-4)  is  recognized  as  the  Hankel  transform 
of  0 (p  ) .  Hence 


0(f  )  =  2tt  J°°  W(w)  J  (27rpw)w  dw  (2. 3“ 5) 

0 

2.4  SPECTRAL  MEASUREMENTS  BY  DIFFRACTION 

Consider  the  system  shown  in  Fig.  2.4-1.  The  source 
of  light  is  a  monochromatic*  coherent  point  source,  meaning  that 
the  relative  phases  of  the  light  at  any  two  points  in  the  opti¬ 
cal  system  are  invariant  with  respect  to  time.  If  the  light  in 
plane  Pi  (Pi  being  any  plane  between  Ll  and  L2 )  is  denoted  by 
E(x,y)  and  the  light  in  plane  P2  is  given  by  Ef(x' ,y'  ),  then 
by  application  of  the  Fresnel-Kirchhoff  diffraction  integral, 
it  may  be  shown  that**  (Born  and  Wolf  1964) 

*  This  requires  that  the  source  has  been  on  for  all  time. 

**  With  the  usual  assumptions  about  small  diffraction  angles. 
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where 


Ef(x',r  )  C  //  E(x,y  )  e~i(2n/‘  f )  (xx1  -t  yy'  )  dx(Jy 

A  (2.4-1 ) 

x,y  coordinates  in  plane  Pi 

x' ,y'  coordinates  in  plane  P2 

C  a  complex  number  with  phase  dependence  on 
x'  ,y 1 

wavelength  of  light  in  free  space 
f  focal  length  of  the  lens 

A  surface  of  aperture 

Pi  any  plane  between  Ll  and  L2 

When  Pi  coincides  with  the  front  focal  plane  (P-', )  of  the  lens, 
C  is  independent  of  x'  and  y' .  However,  in  the  work  pre¬ 
sented  here,  the  phase  of  C  is  of  no  interest  and  only  the 
magnitude  of  C,  denoted  by  K,  is  required  for  the  calibra¬ 
tion  of  experimentally  obtained  data. 

It  is  convenient  to  define  the  transverse  spatial 
frequencies  u , v  by 


u  -  x'  /A f 
v  y'/>f 


(2. 4-2 ) 


Spatial  frequencies  have  the  units  of  cycles  per  unit  distance. 

A  convention  that  sometimes  appears  in  the  literature  is  to 
ascribe  units  of  "lines"  per  millimeter  to  u  and  v.  Plane 
P2  is  referred  to  as  the  frequency  domain.  For  the  work  here 
it  is  most  convenient  to  express  the  light  in  P2  in  terms  of  the 
frequencies  (u,v),  giving 

E  (u,v)  -  K  ff  E(x,y)  e”  j2rr  (ux-fvy )  dxdy  (2.4-3) 

A 
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It  will  now  be  demonstrated  that  with  properly  chosen 
parameters ,  the  system  shown  in  Fig.  2.4-2a  can  be  used  to 
measure  the  power  spectral  density  of  the  complex  process 
j  Rfx.y''  |  associated  with  the  transparency  T.  If  the  total 
power  of  the  coherent  beam  is  J  ,  the  complex  amplitude  of 
the  coherent  wave  E^/x.y)  incident  upon  the  object  T  is 
given  by 

EI(x,y)  =  ';i  rect  £  rect  ^  ,  (2.4-4) 


where  L  and  W  are  the  aperture  dimensions,  and 

rect  f  ^  1  |f|  ^  1/2 

-  0  otherwise 


(2.4-5) 


The  light  ET(x,y)  transmitted  by  the  object  is  given  by 


E  (x,y)  =  E_(x,y)  R(x,y), 


(2.4-6) 


and  the  light  Ef(u,v)  in  the  uv  plane  is,  according  to 
Eq .  (2.4-3),  given  by 

Ef(u,v)  -  K  J  J  E  (x,y)  e"^'r^UX+Vy'  dxdy. 


(2.4-7) 


where  the  aperture  surface  comprises  a  rectangle  of  area  LW. 

From  Parseval ' s  theorem, 

/"  /"  |Ef(u,v)|  2  dudv  -  K2  /+“  /+"  iBjtx.y)!  2  dxdy 


—  00  —00 


(2.4-8) 


Substituting  the  expressions  for  u  and  v  in  terms  of  x'  and  y' , 
Eq.  (2.4-2)  gives 


(l/>f)  ;",c'  I  Ef  (x*  /?.f,y'  /Af  )|  "  dx’dy’ 


,f+“  J+00  |  ET(x,y  )|  dxdy 


—  oo  —  oo 


(2.4-9) 


However,  if  the  lens  is  lossless,  the  integrals 


/+°°  /<(  |  Ef (x'  /'Af ,y'  /?, .f  |  dx'dy'  and 


/+“  /+0°  |  ET(x,y)  |  "  dxdy 


—  ro  —00 


must  be  equal  since  they  represent  the  total  power  on  each 
side  of  the  integrating  lens  L2 .  Hence, 


K  =  l/>f 


(2.4-10) 


and  the  light  amplitude  in  the  uv  plane  is  given 


!u,v)  =  lAf  (j-y 7;  R(*,y) 


(2.4-11) 


where  from  Eq.  (2.1-7) 


R(x,y)  -  RQ  +  r (x ,y )  Rq  ~  <  R(x,y)  > 


The  intensity  in  the  uv  plane  is  given 


I' f(u,v)  =  E  (u,v  )  E  *  ( u , v ) 


(2.4-12) 


Observe  that  I'(u,v)  is  a  random  variable?  its  expected  value  is 


TTW  2  2  2 

<f  I'  f  (u,v)>  =  2  —  |R  |  (sine  Lu )  (sine  Wv ) 

>  f 


+  <I'n(u,v)  > 


(2.4-13) 
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where 


sine  x 


sin  /T  x 


The  term  <'I’n(u,v)  y  is  the  contribution  of  the  noise 
and  is  given  by 


I'n(u,v) 


fW/2  l/2j»/2  0  (x-#.y-M) 
A  f  LW  -L/2  -W/2  -L/2  -W/2 


e-j27T  [u(x-P.)  +  v(y-n)]dxdyd(d1i 


(2.4-14 ) 


With  the  transformation 


=  x  -  r 

=  y  -  T) 


(2.4-15) 


the  integral  (I^UjV)  y  can  be  reduced  to  a  double  integral 
The  Jacobian  of  the  transformation  is 


(g.g.’v  _ 

x,y,f  ,t)  j 


(2.4-16) 


and  the  region  of  integration  is  a  four-dimensional  cube  which 
can  be  described  as  follows:  With  a  and  f3  fixed,  the  limits 
of  integration  on  7  and  6  are  (Fig.  2.4-3) 


for  a  0 
for  a  /  0 
for  6  >  0 
for  P  <  0 


a  -  L/2  <  7  <  L/2 

-  L/2  <  7  <  a  +  L/2 
6  -  w/2  <  6  <  w/2 

-  W/2  <  b  <  6  +  W/2 
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The  limits  of  integration  on 


a  and  p  are 


-L  /  a  /  L 
-W  /  p  <  W 


If  the  relative  average  noise  intensity  is  defined  by 


then 


^n(U'v)  > 


^I'n(u,v)> 

j/>2f2 


(2 


/I  (u,v)>  =  2/L  /W  (1  -  a/L)  (1  -  p/w )  f  (a  ,p,u,v  ’ 
n  0  0 

+  fL  J°  (1  -  a/L)  (1  +  p/w)  f  (a.  . B ,u, v ’ 

0  -W 

+  f°  fW  (1  -  a/L )  (1  -  P./W)  f  (a,  ,p  ,u,  v  ’ 

-L  0 

-I  f°  f°  (1  +  a/L)  (1  +  0/W)  f  (a  ,P  ,u ,  v 

-L  -W 


where 


f  (a.  ,p  ,u,v)  r  4>r(n  ,P  )  e 


-  j2;r  (au+pv) 


Using  '/■  r  ( a  ,  P  )  =  4>r*  (-a,-p)  and  tf>rVa,-p)  =0r*(-  a,p), 

v»rite  Eq.  (2.4-18)  as 


<1  (u,v))  =  2/L  /W  (1  -  a/L)  (1  -  p/W)  Re[  f  (a  ,p , 
0  0 


+  f (a  ,-p,u,v) 


]  dadP 


(2 
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4-17) 

dadB 

da  dp 
da  dp 

)  da  dp 

.4-18) 

we  can 

»v) 

.4-19) 


+L  ,+w  -U2ir(ua+v0  ) 

/  /  0r(a,fi)  e  dadp 

-L  -W 


2 fL  fW 
0  0 


(a/L  +  0/W  -  a0/LW)  Re[ f (a ,(3 ,u, v ) 


+  f(a,-0,u,v)l  dad0 


(2.4-20) 


Hence,  in  the  limit  as  L  and  W  become  arbitrarily  large,  we  obtain 


lim  I  iu,v)‘-  =  W  (u,v) 
L,W->« 


-  lim  (  i/L  f  fi/w  -  otp/LW) 

L,W-o  0  ‘.'i 

.  Re  [  f  ( a  ,  i  J< ,  u  ,  v  )  +  f(a.-0,u,v)l  da.d0 

(2.4-21) 

If  the  magnitude  of  the  covariance  0r(a,lO  tends  to  zero  at 

a  rate  greater  than  |  ctP  |  ^  1+'>'  \  where  A  )>  0,  in  the  region 

where  a  is  greater  than  some  and  p  is  greater  than  some 

0  .  the  second  term  of  Eq.  (2.4-21)  becomes  zero  and 
m’  ^  '  ' 


lim  I  (u,v>  =  W  (u,v) 
L  ,W^oo  * 


(2.4-22) 


When  the  covariance  function  is  zero  or  negligible  at  all  points 
outside  a  region  which  is  small  compared  with  the  area  of  the 
aperture,  an  excellent  approximation  to  an  estimation  of  Wr(u,v) 
may  be  obtained  with  a  finite  aperture.  The  effect  of  finite  aper¬ 
tures  on  the  measurement  of  spectral  densities  is  investigated  in 
Chapter  3 • 

The  astigmatic  system  shown  in  Fig.  2.4-2b  can  be 
used  for  obtaining  one-dimensional  spectral  densities.  If  the 
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length  of  the  aperture  is  D  ,  the  relative  average  light 
intensity  as  a  function  of  u  is  given  by 


2  2 

If(u)N  =  D  [  Rq  |  sine  uD  +  <'In(u)). 


(2.4-23) 


where  /In(u)x  is  the  relative  noise  intensity  given 


[  (u)N  =  1/D  f*D/ 2  r+D/2  0  (x-f)  e^2"  u(x-n 
n  '  -D/2  -D/2  r 


dxd? 


(2.4-24  ) 


With  a  =  x  -  ('  and  ft  =  x,  Eq.  (2.4-24)  can  be  written  as 
the  single  integral 

<'I_(u);  =  2  /°  (1  -  a/D)  Re(0  (a)  e~^2nUa]  da 

n  v  r 


r+D  \  -32 nun 


-  ?/  a/D  Re(0  (a ) 


rj27TUrt)  da 


(2.4-25) 


If  there  exists  an  interval  [-a,,  ,a„]  such  that  for  all  Ja|  V, 

M  MJ  M 

|  0(a)  |  tends  to  zero  at  least  as  fast  as  a.  '  where  A  /0, 

then 


lim  In(u)  =  / +0°  0  (of  )  e”^27T  m  dti  =  W(u). 


(2.4-26) 


2.5  SPECTRAL  MEASUREMENTS  OF  THE  REAL  PROCESSES  lt(x.vd 
AND  ] V (x,yj| 

When  an  object  with  complex  transmittance  R(x,y)  = 
t(x,y)  e^7f  X,YS>  is  inserted  into  a  liquid  gate  with  refractive 
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index  matched  to  the  refractive  index  of  the  object,  the  opti¬ 
cal  path  length  for  the  coherent  wave  is  constant  across 
the  aperture  and  the  relative  phase  ^(x, y)  =  0.  Hence, 

R(x,y)  -  t(x,y),  where  t(x,y)  is  a  real  random  variable 
with  range 


u  /  t(x,y)  /  1 
t(x,y)>  -  tQ 


(2.5-1) 


When  a  liquid  gate  is  used  in  the  system  shown  in  Fig.  (2.4-1), 
the  relative  light  intensity  in  the  uv  plane  is 


l(u,v)s  =  LW  t  2  sinc2uL  sine2  vW 
o 

,  1  /LW  fL/2  ?/J  / ^  ftj  *  (x-f  ,y-it) 

-L/2  -W/2  -L/2  -W/2 

.  e  i2',ruvx  $  )+v(y  ri )  dxdyd^dT] 

(2.5-2) 


where 


</>T(x-^y-Tl)  -  </(t(x,y)  -  tQ)  (t(f,r))  -  tQ)  > 


?  <  *r(x,y  )  r  (  f  , rj )  > 


(2.5-3) 


If  the  aperture  dimensions  are  very  large  compared  with 
the  largest  intervals  over  which  the  random  process  has  signi¬ 
ficant  correlation  then  to  a  very  good  approximation 

4l(u,v)  \  -  LW  tQp  sine2  uL  sine2  uW  +  W^(u,v) 

(2.5-4) 
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where  W  (u,v)  is  the  power  spectrum  of  the  process  j  r  ( x ,  y  )  |  . 

The  direct  measurement  of  the  power  spectrum  of  the 

process  |"/(x,y)}  is  not  so  easily  obtained.  The  reason  for 

* 

this  is  that  it  is  generally  not  possible  to  relate  the  power 
spectrum  W  (u,v)  directly  to  the  spectra  of  {x(x,y)|  and 
|V(x,y)|  even  if  |r(x,y)J.  and  |V(x.yl(  are  independent 
random  processes.  To  see  this,  observe  that  the  spectra 
W^(u,v)  and  W  (u,v),  of  {r(x,y)|  and  |,'(x,y){,  respectively, 
depend  only  on  the  second  moments  of  |r(x,y1}  and  {.(x,y)| 
while  W  (u,v)  normally  depends  on  all  the  moments  of  |’/(x>y)}- 
When  |'.’(x,y)}  and  |r(x,y){  are  independent,  it  follows  that 


R(x,y)  t(x.y)  .W(*.y),  H  ,  {x>y  (x.Y  1 


and 


2.‘  1 


R ( x * y )  N  -  R0  =  toM„  (l) 

xy 


(2.  -i.l 


where  M,,  (  ■  ) 

v  xy 

of  ',’(x,y).  Since 


j  V .  I  (  x  V  )  ‘ 

eJ  Y  is  the  characteristic  function 


r(x,y)  -  T(x,y)e^X’y)  4  t  fe^’(x,yl  -  m,  (1)1 

o  V  * 

Y  xy 

(2.1-7) 

the  covariance  ,y-T])  is  given  by 

0r(*-; .y-n)  0 x(x_ '  •  y-  t )  m ,  .  (i,-i) 

^  xy  ?  ( 


t  2[ 

o  1 


M 


*xy  *1 


(1,-1)  - 


7  (D 

f  fj  Y  xy 


(2. 1,-8) 
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where  M  ,  (  • 

•  xy '  t  r( 

of  f(x,y  ),V  (f  .  r,) 


/  )  is  the  joint  characteristic  function 

and  is  given  by 


M., 


=  <  e 


xyv  f  , 


j[>.V(x,y)  +  '  (2.5-9) 


Now  M  (1,-1)  can  be  expressed  as 

s  xy’  •  •, 

.k 

^[V'(x»y)  -  f(?iri)  ]k  >  — - 


M.(  ,,,  (1,-1),= 


v  ....  •. 


xy  ►  , 


k-0 


•(P.5-10) 

which  clearly  indicates  the  dependence  of  0r  (x- £  ,y-Tj )  on 
all  the  moments  of  y  ( x  ,  y  ) . 


In  the  special  case  where  the  phase  errors  are  small, 
a  possibility. for  example, when  the  matching  liquid  in  a  film 
gate  has  a  refractive  index  slightly  different  from  the  re¬ 
fractive  index  of  the  transparency,  and  |t(x,y)|  and  |0(x,y)} 
are  independent  processes,  there  results 


t  v  x ,  y )  t  ( ? 


n)  ej  r^(x,y )-V'(P ,T|)] 


> 


=  [to2  +  ]  .  [i-0^(o,o) 

+  y-h)  ] 

where  0  ,  (x-i‘  ,y-i( )  =<  f(x,y)  ip  (f,ri)  \  .  Hence  the  relative 
light  intensity  in  the  uv  plane  is  given  by 

<l(u,v)>~  tQ2  [l-0^  ( 0,  0 )  J  LW  sinc2uL  sinc2vW 

+  [1-0V,(O,O)]  Wt  (u,v)  +  Wt  * 

+  V  Vu>v) 
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where  the  asterisk  denotes  convolution.  Thus  in  addition  to 
the  attenuated  zero-order  light,  the  average  light  intensity 
contains  terms  that  are  proportional  to  the  spectra  of  |~(x,y)| 
and  •|V'(x*y)j  and  their  convolution.  The  presence  of  the 
convolution  term  prevents  the  direct  measurement  of  |w^(u,v)| 
even  if  W  u.v)  is  known  from  an  independent  measurement. 

More  will  be  said  about  the  random  processes  |^(x,y)|  and 
|t(x,y)|  when  discussing  theoretical  models  of  film  noise  in 
Chapter  . 
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3.  MEASUREMENTS  OF  POWER  SPECTRA  WITH 
FINITE  SYSTEMS 


In  the  previous  chapter  it  was  shown  that  a  diffraction 
system  with  an  infinite  object  plane  aperture  produces  in  the  uv 
plane  a  light-intensity  distribution  whose  average  provides 
an  unbiased  estimate  of  the  power  spectrum  of  a  random  process. 
Spectral  density  measurements  with  a  real  system  always  en¬ 
tail  the  use  of  finite  apertures  and  scanning  slits.  The  princi¬ 
pal  effects  of  these  components  are  to  reduce  the  sensitivity 
of  the  system  and  to  bias  the  estimates  of  the  true  spectral 
density . 

It  is  generally  desirable  to  make  spectral  density 
measurements  with  large -aperture  systems.  A  large  aperture 
reduces  the  bias  of  the  measurement,  allows  the  observation 
of  correlation  effects  over  large  intervals,  and  improves  the 
ratio  of  scattered  spectral  light  to  unwanted  zero-order  light. 

If  the  power  spectrum  were  available  in  the  uv  plane, 
then  ideally  the  scanning  slit  should  be  an  infinitely  small 
pinhole.  Such  a  pinhole  would  faithfully  transmit  the  power 
spectrum  to  a  photomultiplier  tube.  However,  in  practice,  the 
light  intensity  is  an  observation  upon  a  random  variable  whose 
expected  value  approximates  the  true  spectrum.  Hence  some 
smoothing  is  desirable.  The  ideal  slit  in  this  case  is  a  slit 
which  is  large  enough  to  smooth  the  random  fluctuations  in  the 
light  intensity  yet  not  so  large  as  to  smooth  significant  changes 
or  peaks  in  the  power  spectrum. 

3.1  EFFECT  OF  A  FINITE  APERTURE;  ASTIGMATIC  SYSTEMS 

The  relative  noise  intensity  given  in  Eq.  (2.4-25)  and 

Vi 

obtained  with  a  one-dimensional  diffraction  system  can  be  writ¬ 
ten  as 
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MOMMA  WKI&II 


Vu  x  = 


/+D  (1  - 


I  ft  I  x  ,  /  «  -j27rua  , 
-FT*  )  'Mu)  e  dr< 


'+°°  (1  -  L-^--)rect  a/2D  0(a)  e  ^,Ul  da. 


-  W(u) 


(3.1-1) 


where  0(a)  is  the  covariance  of  the  process,  rect  a./cD  is 
defined  by  Eq.  (2.4-5), and  W(u)  is  the  system-generated  power 
spectrum  or  apparent  power  spectrum  of  the  random  process. 

The  Fourier  transform  of  Eq.  (3.1-1)  is  the  convolution  of  the 
true  power  spectrum  and  the  Fourier  transform  of  the  triangu¬ 
lar  function 


=  (1  -  ~^“)  rect  ci/2D 


(3.1-2) 


Hence , 


=  D  f  W  (£)  sine5  (u-OD  dT 


(3.1-3) 


where 


sme*- 


(3.1-4) 


and  D  sine2  uD  is  the  Fourier  transform  of  h(n. ).  Equation 
(3.1-4)  clearly  describes  the  effect  of  the  finite  aperture  of 
the  spectrum.  W(u)  is  the  result  of  averaging  over  frequency 
the  spectrum  W(£ )  with  weights  proportional  to  H(u-f)  -  D  sine2 
(u-f)D.  In  the  terminology  of  communication  theory  (Blackman  & 
Tukey  1958),  H(u)  is  the  spectral  window  corresponding  to  the 
spatial  lag  window  h(m).  From  Eq.  (  3 . 1—  3 )  it  is  seen  that 
W(u)  is  a  biased  estimate  of  W(u).  The  degree  of  bias 
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depends  on  D  and  on  W(u).  For  large  apertures  the  signi¬ 
ficant  averaging  over  frequency  is  confined  to  frequencies 
near  u,  and  W(u)  would  be  expected  to  be  a  better  replica 
of  W(u)  than  if  small  apertures  were  used,  for  which  the 
averaging  would  extend  over  a  large  frequency  spread.  Since 
over  90'  of  the  area  of  sinc2(u-£)D  lies  under  the  main  lobe 
and  the  largest  sidelobe  (the  first)  is  1 3  db  below  the  peak 
of  the  main  lobe,  most  of  the  averaging  is  done  over  the  fre¬ 
quency  interval  -ID  /u-£  /  1  7D .  Nevertheless,  the  contribu¬ 
tions  of  the  sidelobes  are  undesirable  since  they  involve 
frequencies  "far"  removed  from  u  =  f  and  thus  tend  to  dis¬ 
tort  the  spectrum.  The  problem  of  confining  the  smoothing  to 
the  main  lobe  is  tantamount  to  finding  a  function  with  a  low 
sidelobe  response  and  a  reasonably  narrow  main  lobe.  This 
problem  has  been  extensively  treated  by  several  authors  (Dolph 
1946,  Jacquinot  and  Roizen-Dossier  I966).  In  optics,  side¬ 
lobe  reduction  can  be  achieved  by  apodisation.  Apodisation 
may  be  accomplished  by  using  either  light-absorbing  aperture 
functions  (also  called  pupil  functions)  or  contoured  apertures. 
Advantages  and  drawbacks  of  the  methods  are  discussed  in  Chap¬ 
ter  6. 

In  the  measurement  of  power  spectra  by  diffraction,  a 
principal  reason  for  apodising  is  to  reduce  the  zero-order  dif¬ 
fracted  light  of  the  aperture,  i.e.,  the  term  that  always  ap¬ 
pears  with  /In(u)x  or  <'In(u,v))>  as  in  Eqs .  (2.4-23)  or 
(2.4-13).  This  unwanted  light  must  be  suppressed  if  observa¬ 
tions  are  to  be  made  on  the  power  spectrum . 

When  it  may  be  presumed  that  the  least  distance  over  which 
there  are  significant  changes  in  W(u)  are  small  compared  to 
1/D,  the  form  of  W(u)  is  essentially  reproduced  in  W(u). 

When  D  is  made  arbitrarily  large, 

lim  D  sine2  uD  =  <*  u  =  0 

D-hoo  =  0  u  ^  0 


and 

lim  W(u)  =  J4°°  W ( f ’ )  (u-^)  -  W(u) 

D  *CV  —00 

(3.1-5) 

The  bias  introduced  by  a  finite  aperture  is  the  differ¬ 
ence  between  W(u)  and  W(u).  It  will  now  be  shown  that  the 
bias  depends  on  the  Hilbert  transform  of  W(u).  Since 
</!*(■-*)  ~  0(-<),  it  follows  that 

W(,u)  -  j  0(a)  e  J  da  -  f  4-  {^  )  e  da 

—  r*  —oo 

(  ’.1-6) 


The  truncated  spectrum  W°  (u)  is  defined  by 

W°(u)  -  /+D  0(a)  e-i27rm  da 

-D 

and  the  truncation  error  a(u,D)  is  defined  by 

n  (u , D )  -  W(u)  -  w°(u)  -■  2Re  /  .'"(a)  e+^2vTUft  d, 


(3.1-7) 


(3.1-8) 


where  o(u,or)  0.  The  Hilbert  transform  of  W  (u)  is  given  by 


WD(u 


f  0  *  (a  )  e 
0 


j  (2nua-  /r/2  )  d( 


,0  , 

!  I  <t> 

-D 


(a)  e  j  ( 2;tUTt  /i/2 )  dj 


2Re  /  (-j)  e 

0 


j2.ru» 


c/> "  (  a. )  da 


(  -M-9) 


where  A  denotes  the  Hilbert  transform. 

Hence  from  Eqs . (2 . 4-25 U 3 . 1-7  )  >and ( 3 . 1~9  )there  results 
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(3.1-10) 


«<»)  "“(u)  -  573  ^  «“(») 

and  the  bias  is  given  by 

W(u)  -  W(u )  -  o(u,D)  +  A  ^(u)  “  a(u>D-) 

(3.1-H) 

Equation  (3.1-11)  shows  that  the  error  in  measuring  the 

/N 

spectrum  results  from  two  effects.  First,  observe  that  W(u) 

is  independent  of  D  and  therefore  —  p  du  ^(u)  decreases 

inversely  with  increasing  D.  This  error  could  be  called  the 
resolution  error  since  it  is  independent  of  the  truncation  of 
V ( 0 .  The  terms  involving  o(u,D)  and  [c(u,D)]  are  trun¬ 

cation  errors  since  they  vanish  when  D  is  chosen  large  enough 
to  include  all  the  correlation  intervals  of  the  random  process. 

For  example,  suppose  that  0(a)  is  zero  (or  at  least  negligibly 
small)  outside  an  interval  A,  centered  at  the  origin  with 

2D .  Then  the  limits  in  Eq.  (3.1-7)  can  be  changed  to  (-00,00) 
and  o(u,D)  :  0  for  all  u.  In  this  case  there  is  no  con¬ 

tribution  to  the  bias  from  truncation  although  there  is  still 
a  bias  due  to  the  resolution  error.  In  the  spatial  domain 
the  bias  is  manifested  as  a  nonuniform  attenuation  of  0(a). 

This  may  be  seen  as  follows.  The  Fourier  transform  of  Eq.  (3.1-3) 
gives 

0(a)  =  0(a)  h (a  )  (3.1-12) 

where 

h(a  )  =  (1  -  L5"^-)rect  ~ 
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indicating  that  <p  (a  )  can  be  recovered  from  the  apparent  co- 
variance  0(a)  from  the  inverse  operation 

v  (a  )  --  0  (ot )  h  1  (a  )  (3. 1-13) 

This  situation  is  shown  in  Fig.  3. 1-1 a  . 

When  the  truncation  error  is  not  zero,  0(a)  cannot 
be  recovered  from  0(a)  since  there  is  no  information  about 
the  covariance  for  intervals  greater  than  D.  This  situation 
is  shown  in  Fig.  3.1-lb  . 

3.2  EFFECT  OF  A  FINITE  APERTURE;  ORDINARY  SYSTEMS 

The  situation  in  two  dimensions  is  essentially  a  direct 
extension  of  the  one-dimensional  case.  Equation  (2.4-18)  can  be 
written  as 


I  (u,v) 
nv  ’  ' 


.i"1"  /4“  (i  - 

—or  “oo 


a  I  .  a  ,  1st.  S 

rect  —  (1  -  rect— 


.  0  (a  ,h  )  e 
-  W(u, v ) 


2L  v  -  W  ' - 2W 

-j2 TT  (urHV|i)  did  ' 

(3.2-1) 


where  0  (a  ,ft  )  is  the  covariance  of  the  two-dimensional  random 
process,  W(u,vl  is  the  apparent  spectrum, and  all  other  quan¬ 
tities  have  been  defined  earlier.  Two-dimensional  Fourier 
transforms  have  properties  similar  to  one-dimensional  transforms 
(Sneddon  lq^l  )  .  Thus,  the  above  may  be  written  as  the  convolu¬ 
tion  of  the  true  spectrum  W(u,v)  and  H(u,v)  -  LW  sine2  uL 
.  sine2  vW.  This  gives 

W(u,v)  =  W(u,v)  *  H(u,v) 

~  LW  j‘+oc,j+0°  w(f,i|)  sine2  (u -£)  L  sinc2(v-ri)W  df)  di, 


1 
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a)  DISTORTION  OF  THE  COVARIANCE  IN  ABSENCE  OF  A 
TRUNCATION  ERROR 


b)  DISTORTION  OF  THE  COVARIANCE  IN  PRESENCE  OF  A 
TRUNCATION  ERROR 


A-006-1  -00-0061 


FIG. 3.1  -I  EFFECT  OF  FINITE  APERTURE  (WIDTH  =  D)  ON  COVARIANCE  ESTIMATE 
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The  apparent  covariance,  0(a,f3)  -  F-  1  |w(u,v)>  is ,  from 
Eq.  (3.2-1), 

0( f»  )  -  (1  -  (1  -  rect  rect  0(a,P  ) 

-  h(u  ,3  )  «/>(a,P  )  (  3.2-  -]) 


When  0(a,O)  is  non-negative,  0(a,P)  is  always  less  than 
0(a,P)  except  at  nulls  and  the  point  a  =  p  =  0,  where 

0(0,0)  =0(0,0).  Hence, 


j‘+on  j'4°  W(u,v)  dudv  J400  J4“’  W(u,v)  dudv  (  ^.2-^) 


The  function  h(a,P)  is  shown  in  Fig.  3.2-1  and  is  a  pyra¬ 
mid  of  height  unity  and  with  base  dimensions  2L  and  2W . 

If  0(a  ,3)  is  zero  for  all  points  outside  the  base  of  h(r/,p)} 
that  is,  for  points  |  a  |  )  L,  |  p  |  P  W,  then  by  analogy  with 
the  one-dimensional  results,  there  is  no  truncation  error  and 
the  true  spectrum  can  be  recovered  from  the  measured  spectrum 
by  digital-computer  filtering.  For  example,  W(u,v)  can  be 
obtained  from 


W(u,v  ) 


F 


J UjV 


h(a ,0  ) 


ll 


I  nt  |  P  L 

1 3 1  <r  w 

(3.2-5) 


Optical  spatial  filtering  of  W(u,v)  is  investigated  in  Sec. (3.5). 


From  Eq.  (3.2-2)  it  can  be  seen  that  L  and  W  must 
both  be  large  for  W(u,v)  to  be  a  reasonable  replica  of  W(u,v). 
Thus,  for  example  as  an  extreme  case,  if  only  L  is  large  while 
W  is  small,  W(u,v)  would  be  given  by 
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CUTTING  PLANE 
PARALLEL  TO  BASE 


<£ ( a,/3 )  =  h(a,/3)  <*>(a,£) 

WHERE 

£(a,/3)  =  EXPERIMENTALLY  GENERATED 
COVARIANCE 

<p  (a,/S)  =  TRUE  COVARIANCE  OF  PROCESS 
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FIG.3.2-1  ILLUSTRATING  THE  FUNCTION  h(o,/S)  OF  EQUATION  ( 3.2 -3  )  FOR 
A  RECTANGULAR  APERTURE  OF  WIDTH  W  AND  LENGTH  L 
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W(u,v)  ~  W 


/+°°  /+°°  W(£,r))  6(u-£)  sine2  (v-T))W  d£dr) 


=  W  /+0°  W(u,r,)  dn 

-oo 


(3.2-6) 


where  it  has  been  assumed  that  over  the  frequencies  for  which 
W(£,t])  is  not  negligible,  W  sine2  (v-r))W  is  essentially 
constant.  With  both  L  and  W  arbitrarily  large 

W(u,v)  -  /+°°  W(^,rj)  6(u-0  6  ( v-T) )  d^dq 

—oo 

=  w(u,v)  (3.2-7) 

The  effect  of  the  aperture  function  H(u,v)  is  to  aver¬ 
age  the  true  spectrum  over  frequencies  centered  about  u  =  £ 
and  v  =  t)  with  weights  proportional  to  H(u-£,  v-ti).  If  L-1 
and  W-1  are  small  compared  with  the  least  distances  by  which 
W(u,v)  changes  by  any  significant  amount,  w(u,v)  will  be  a 
reasonably  good  replica  of  W(u,v).  The  degree  of  bias  is  again 
a  function  of  the  truncation  and  resolution  errors.  Before 
computing  the  bias  error,  W(u,v)  -  W(u,v),  it  is  convenient 
to  introduce  some  transforms  which,  by  direct  analogy  with  the 
ordinary  Hilbert  transform,  will  be  referred  to  as  two-dimen¬ 
sional  Hilbert  transforms. 

3.3  TWO-DIMENSIONAL  HILBERT  TRANSFORMS 

This  section  is  a  brief  digression  from  the  main  work 
of  this  chapter.  In  what  follows,  certain  Hilbert  transform- 
type  integrals  are  derived  for  functions  of  two  variables. 

The  results  of  this  section  will  be  used  in  the  following  sec¬ 
tion  to  derive  the  bias  in  the  system-generated  spectrum  W(u,v). 
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The  notation  used  in  this  section  is  unrelated  to  the  notation 
used  elsewhere  in  this  chapter.  Thus,  for  example,  the  reader 
should  refrain  from  associating  (a,£)  with  the  spatial  "lag" 
domain  a.id  (u,v)  with  the  spatial  frequency  domain. 

Although  the  Hilbert  transform  has  been  used  in  Sec.  3.1  > 
it  is  best  to  begin  this  section  with  its  formal  definition. 

The  ordinary  Hilbert  transform  of  a  function  x(a)  is  defined 
as  the  Cauchy  principal  value  (C.P.V.)  of  the  integral 

x(a)  =(1/7 r)  /+0°  dA  (3.3-1) 

—  00 

/s 

where  x(a  )  denotes  the  Hilbert  transform  of  x(a)  (Titchmarsh 

A 

1948).  The  function  x(a)  can  be  recovered  from  x(cx)  by 
the  C.P.V.  of  the  integral 

x(a.)  =(-1/tt)  /+00  giAl  d*  (3.3-2) 

“00 

If  X(u)  denotes  the  Fourier  transform  of  x(a)  and  Y(u) 

A 

denotes  the  Fourier  transform  of  x(a),  then 


Y(u)  =  -jX(u) 
=  +jX(u) 


u  >  0 

u  <  0 


(3.3-3) 


The  substitution  of  Eq.  ( 3 . 3— 1 )  into  Eq.  (3.3-2)  gives 


x(a)  .(-1/n ')/-*“  d*x(>)  !*"  U-rflr- a)  (3.3-4) 

—00  —00  '  ^  ^  ' 

Hence , 

W>C  c  6<A-a)  ( 3 . 3" 5 ) 
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I 


3 


and 


of  the  fol- 


Consider  now  the  C.P.V.  about  a  =  f, 
lowing  integrals 


6  =  n 


X 


(«,.')  -(1/0  /*“  xfc.n) 


—  oo  — oc 


(a 


-eTffTn)  (3.3-6) 


g(a,p)  =(1/.A)  /+”  /+"  (a;^' dCd,, 


(3.3-7) 


—  oo  — oo 


Define  the  following  quantities  (F  |  |  r  Fourier  transform) 


L (u , v  )  =  F 

1  '  1 

(  7T2  m3  | 

X (u , v  )  =  F 

|x(a,p ) | 

Y  ( u  ,  v  )  =  F 

)x(a,3  )} 

G(u , v  )  =  F 

jg(«,o )  |  . 

For  L(u,v)  there  results 

L (u, v )  .(!//  )  /+“  /+“  e-32.i(ua+vp) 


—  on  —  oo 


/,  /.  i  ,+oo  sin  2tt  ua  ,  ,  ,+oo  sin  2tt  v(3 

(1/ Dir  )  /  - * -  da  ,  I/jtt  /  - g - u  dfi 


(3.3-8) 


Since 


(  1/j  7T)/ 


+oo  sin  tx 
x 


dx 


-  D 

+  j 


t  >  0 

t  ^  0 


(3.3-9) 
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it  follows  that 


L(u,v)  -  -1 

=  +1 


uv  3  0 

uv  <  0 


(3.3-10) 


If  the  Fourier  transform  is  taken  of  both  sides  of  Eqs  .  (3.3-6) 
and  (3.3-7),  then 

Y ( u , v )  =  X(u,v)  L(u,v)  (3.3-11) 

G(u,v)  =  Y(u,*)  L(u,v)  (3.3-12) 

Hence  G(u,v)  -  X(u,v)  L2(u,v)  =  X(u,v),  and  it  is  established 
that 

xK6)  g(a  ,6  )  (3.3-13) 

This  result  could  also  have  been  achieved  by  direct  substitu¬ 
tion  of  Eq.  ( 3.3-6)  into  Eq.  (3.3-7)  and  using  Eq.  (3.3-5). 

Equations  ( 3 . 3—  6 )  and  (3.3“7)  are  the  two-dimensional 

£ 

analogs  of  Eqs.  (3.3-1)  and  (3.3“2)  and  hence  x(a,P)  can 
be  viewed  as  the  two-dimensional  Hilbert  transform  of  x(or.,6). 
For  what  follows  it  is  also  useful  to  define  the  quantities 

xa  (a  ,6  )  and  x^(a,6  )  by 


*o<° 

,6) 

=  (1/tt)/+" 

—00 

x(a  .6  ) 
a-C 

d£ 

(3. 

3-14) 

,P) 

=(1  A)/+" 

—00 

x(a  ,e ) 
6-C 

d£ 

(3. 

3-15) 

3.4  THE  BIAS  ERROR  IN  THE  APPARENT  SPECTRUM  Wfu.v) 

The  results  of  the  previous  chapter  enable  the  writing 
of  the  two-dimensional  Hilbert  transform  of  the  true  power 
spectrum  W(u,v)  as 


-53- 


where 


A 

/v 


W(u,v) 


-2  f°  Ja  Re  [  f  (nr  ,3  ,u,v)  -  f  (a  , -3  ,u,  v  )  j  dadr 

u  0 

(3.4-1) 


0  (  t  ,1 ■  ) 


e~j2n  (un+v1/  ) 


Assuming  that  differentiating  inside  the  integral  sign  is  per¬ 
mitted,*  there  results 


1 

(  2  rr ) 2  LW 


?  WfUjV ) 
dUov 


Re  [  f  (a,('i  ,u,v) 


•if  (a  , -3  u,v)]  dadr 


(  ;>"2) 


In  deriving  Eq.  (  -1.4-1)  use  was  made  of  the  fact  that  W(u,v) 
could  be  written  as 


W(u,v) 


.-for  .-fr* 

I  i 

.■  J 

“l"  ™  Ou 


J2 


( ua  +\  3  ) 


dad: 


(3.4-3) 


Similarly, 


oW  ( u , v ) 


u 


2t  L 


iU 


2  r  f'  ~T~  Re  [ f ( r ,3 ,u , v ) 

0  0 


-I  f  (a ,-p  ,u , v )  ]  dad(3  (  Vl-4  ) 


*  This  requires  that  the  integrand  in  Eq.  3.4-1  has  section- 
ally  continuous  partial  derivatives  with  respect  to  u,v  and 
the  integrals  in  Eqs . ( 3 . 4-1  )and( 3 .4-2 )converge  uniformly 
(Courant,  Vol  III , "Differential  and  Integral  Calculus.") 


and 


1  s<V(u’v>  _  R  _  _ 

27w  ~ jw -  =  2/  /  w  Re  [£(a.P.u,v) 


0  0 


+  f(ra,-£,u,v)]  da  dp 


where 


wu(U,v)  =  (i/ 7r )  /+0° 


(3.4-5) 


(3.4-6) 


wv(u,v)  =  (1/tt  )  /+00 


(3.4-7) 


Let  the  truncated  spectrum  W^L,W^(u,v)  and  the  truncation 
error  o(u,v,L,W)  be  defined  by 


W(L-W>(u,v)  =  /L  /«  0(r,,B)  e-:2'T(uc.+v(j) 


da  dp. 


-L  -W 


(3.4-8) 


a(u,v,L,W)  r  W(u,v)  -  W(L’W)(u,v)  (3.4-9) 

Then  Eq.  (2.4-20)  can  be  written  with  <In(u,v)>  -  w(u,v),  as 
w(u,v)  =  w<L’w>(u,v)  -  ^  WU(L’W>(U,V) 


l  -iL-  M  (L»W),„ 

27T  W  c)V  V  (  ’  ) 

1 _ 72  *(L,W  )  , 

(2tt)2  IA*  SuSv  W  (u’v 


(3.4-10) 
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When  Eq.  (.'.4-  i)  is  used,  the  bias  W(u,v)  -  W(u,v)  can  be 
written  as 


W ( u , v )  -  W(u , v ) 


_1 _ 1 

2  7T  L  C'U 


4 


_1 _ 

2tt  L  c>v 


(U.v) 


1 

(2/r)2  LW 


eUc'v 


W(u,v) 


+  terms  involving  the  truncation  error, 
its  derivatives  and  transforms . 


C3.4-11) 


It  may  be  seen  from  Eq.  (3.4-8)  that  the  truncation  error  is 
identically  zero  for  all  u.v  if  0>(a,’3)  is  zero  outside  a 
region  bounded  by  |  a|  L ,  I  P|  )>  W.  In  this  situation,  the 

only  remaining  error  in  the  measurement  of  W(u,v)  is  the 
resolution  error  which,  in  principle  at  least,  can  be  corrected 
to  yield  w(u,v\  The  spatial  filtering  of  W(u,v)  or  W(u), 
as  the  case  may  be,  is  discussed  next. 

1  •  SPATIAL  FILTERING 

When  the  truncation  error  is  insignificant  compared  with 
the  resolution  error,  the  bias  in  the  measured  spectrum  may  be 
eliminated  by  further  processing.  This  is  implied  by  Eq.  (3.2~5) 
which  is,  perhaps,  the  way  a  digital  computer  might  proceed.  In 
this  section,  the  feasibility  of  spatial  filtering  is  considered, 
and  it  is  shown  that  direct  spatial  filtering  is  not  physically 
realizable . 

Consider  a  one-dimensional  system  first.  The  true  co- 
variance  is.  given  by 

<4(r')  :  g(m)  0(a)  (3.4-1) 


-^6- 


where 


g(a )  =  (1 


I  a|  <f  A  <  D 


0  otherwise 


(3.5-2) 


and  it  is  assumed  that  0(a)  =  0  or  is  insignificant  for 
|  / 1  >  A  <  D.  If  G(u)  denotes  the  Fourier  transform  of 
g(a).  then  the  Fourier  transform  of  Eq.  ( 3 . 5~ 1 )  gives 

W(u)  =  G(u)  *  W(u)  (3.5-3) 

This  operation  can  be  performed  optically  in  the  fre¬ 
quency  plane  by  collecting  W(u)  through  a  filter  with  inten¬ 
sity  transmission  function  G(u)  and  integrating  over  the  u 
domain  with  a  photomultiplier  tube.  In  order  for  the  operation 
expressed  in  Eq.  (3.5“3)  to  have  any  physical  significance, 

G(u)  should  satisfy  the  following  three  conditions: 


(i) 

lim  G(u)  =  0 

U  ->  oo 

(ii) 

G(u)  should  be  real  - 

(iii) 

G (u  )  >  0  -oo  <  u  ^ 

(3.5-4) 


Condition  (i)  requires  that  any  degree  of  accuracy  in  recon¬ 
structing  W(u)  can  be  achieved  by  using  a  filter  of  finite 
dimensions.  Conditions  ( i i )  and  (iii)  reflect  the  fact  that 
W(u)  and  W(u)  are  intensities  and  therefore  are  nonnegative, 
real  quantities. 

That  G(u)  satisfies  condition  (i)  is  easily  demon¬ 
strated  by  the  Riemann-Lebesgue  lemma  which  states  that  if  a 
function  x(a)  is  integrable  on  [cx,p],  then 

lim  x(a'  cos  27T  ua  da  =  0  (3. 5“ 5) 

u  -•*  oo  a 
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I  a  l  _i 

Since  (1  -  )  is  integrable  on  [0,A]  for  A  A  D, 

it  follows  that 


The  function  g(a  )  is  real  and  even.  This  requires  that 
G(u)  be  a  real  function  of  u.  Finally,  to  test  condition 
(iii),  observe  that  at  u  -  1/2A, 

G(1/2A)  =  2D)J 

=  2</2  cos » [-fe  -  --ir-.  a ]dli 

where  a  D/a.  In  the  interval  [0,/i/2j,  cos  3  )  0.  Now 
G(1/2A)  will  be  negative  if  l/frrf'v-l)  +  p  )  )>  l/(;:i-3)  for  all 
3  in  [0,tt/2];  however,  this  will  be  true  if  23/tt  A  1  for 
all  3  in  f0,7t/2].  Clearly  this  is  so,  proving  that 
G(1/2A)  A  0.  Hence,  G(u)  violates  condition  viii)  and  can¬ 
not  correspond  to  a  physically  realizable  filter. 

The  function  G(u)  can  be  written  in  terms  of  tabulated 
functions  (Jahnke  <$•  Emde  194 p)  as  follows: 

G(u)  =  D  cos  27T  uD  [Ci  (  2tt  uD )  -  Ci  (2rr  uDr  )  ] 

-f  D  sin  277  uD  [Si  (2n  uD )  -  Si  (27T  uDP  )  J 

(3.9-7) 

where 


-58- 


(3.5-8) 


Ci (x )  =  -  f °°  co-°  ^  d£ 
x  ^ 

Si(x)  -  f*  dC 

0  ^ 

and  P  =  1  -  A  /D  <  1 . 

G(u)  is  plotted  in  Fig.  3 . 5~1 . 

The  analog  of  Eg.  (3«5“1)  in  two  dimensions  is 

0(a,0)  =  g(a,3)  0(a,3) 


where 

<  L, 

<  w 

(3.5-9) 

The  function  g(a,3)  which  would  rectify  the  distortion  is 

g(a,3)  =  [h(a  ,3  )  J  _1  |a  I  <  A  x  <  L 

13  I  <  A2  <  W 


h(a  ,3 )  =  (1  -  4-L  )  (1  -  LE-L  )  for  |a  | 

131 

-  0  otherwise . 


(3.5-10) 


where  Ai  and  a2  are  intervals  beyond  which  0(a,3)  =  0. 
Clearly  if  G(u,v)  =  F  jg(a,3)|  ,  then 


G(u,v) 


J-+00  J+0°  g(a,3)  e~i27T  (Ua+Vp) 


G (u )  G(v). 


(3.5-H) 


•] 


\ 


i 


i 


.1 

* 


! 


i 
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DIMENSIONLESS  VARIABLE 


But  the  properties  of  G(u)  have  been  determined  in  the  one¬ 
dimensional  analysis,  where  it  was  found  that  G(u)  could  not 
correspond  to  the  space  transmission  function  of  a  physically 
realizable  filter.  Therefore  G(u,v)  is  not  physically  real¬ 
izable  . 

An  optical  filtering  technique  for  reducing  the  resolu- 
tion  inherent  in  W(u)  when  the  sari;  le  of  data  is  small  and 
the  effective  aperture  is  limited  by  the  sample  and  not  by  the 
lens  is  the  following.  Let  W(u)  be  recorded  on  low-noise 

film  so  as  to  let  the  average  amplitude  transmittance  of  the 

-w 

film  be  proportional  to  W(u),  i.e.,  t(u)  =  kw(u), where  t(u) 
denotes  the  amplitude  transmittance  and  k  is  some  constant. 

If  this  transparency  is  placed  in  the  aperture  of  the  optical 
system  shown  in  Fig.  3*5-2,  then  the  Fourier  transform  of  t(u) 
is  proportional  to  0(a),  and  a  filter  in  the  Fourier  transform 
plane (P2)  with  transmittance  g(a)  will  yield  a  light- ampli¬ 
tude  variation  proportional  to  0(a).  In  plane  (P3),  the 
Fourier  transform  W(u)  of  0(a)  is  obtained. 

The  method  described  above  has  several  drawbacks.  First, 
a  very  low- noise  film  such  as  a  Lippmann  emulsion  type  (Mees 
and  James  I966)  is  required  for  storing  W(u).  This  film  must 
also  be  linear  over  a  dynamic  range  sufficient  to  accommo¬ 
date  the  significant  variations  in  W(u).  Finally,  if  the 
sample  from  which  W(u)  was  obtained  is  very  small,  a  large 
aperture  may  be  required  for  good  results . 

3.6  EFFECT  OF  THE  SCANNING  SLIT 

Consider  a  scanning  aperture  (henceforth  called  a  slit 
to  avoid  confusion  with  the  object-plane  aperture)  which  to¬ 
gether  with  a  photomultiplier  tube,  amplifier,  and  recording 
device  enables  the  recording  of  the  light-intensity  distribution 
at  any  point  (x' ,y' )  in  the  frequency  plane.  Recall  that  any 
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FIG.  3.5-2  OPTICAL  SYSTEM  FOR  THE  SPATIAL  FILTERING  OF  THE  MEASURED 

SPECTRUM  W(u,v) 
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point  in  the  frequency  plane  corresponds  to  spatial  frequen¬ 
cies  u  =  x'/>f,  v  r-  y'/xf  [see  Eq.  (2.4-2)  and  Sec.  2.4]. 
If  the  slit  is  a  rectangular  hole  of  width  t  and  length  h, 
the  total  illumination  received  by  the  photomultiplier  tube 
through  the  slit  is  obtained  by  integration,  and  the  current 
in  the  photomultiplier  tube  when  the  center  of  the  slit  is  at 
(x1  ,y'  )  is  given  by 


i  =  k  /y'+b/2  ;*'+  V2  I(f;, n)  dfdi) 

y'-b/2  x'-h/2 


(3.6-1) 


where 


I  ( r,  H  )  =  light  intensity  at  (f,  ,l) 

k  -  gain  of  the  photomultiplier  tube 
i  =  current 

Equation  (3.6-1)  can  be  written  as 


i  ^  k  /+r°  /+°°  I  (x'  ,y*  -Tj )  g'U,T))  d^dr) 

—co  —00 


where  the  function  g’  (£,r|)  is  given  by 


(3.6-2) 


<3'  (f  ,'l)  “  rect  rect 


(3.6-3) 


and  is  the  ratio  of  the  transmitted  intensity  to  incident  in¬ 
tensity  at  the  point  (£,,tj).  The  coordinates  (£,tj)  are  de¬ 
fined  with  respect  to  an  origin  at  the  center  of  the  slit. 
Since  scanning  takes  place  in  the  frequency  plane,  it  is  con¬ 
venient  to  write  g' (£,tj)  in  terms  of  spatial  frequencies. 
Hence,  with  u'  =  t/~h£  and  v'  =  Tj/xf,  one  may  write 
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g(u'  ,v<  )  -  1, 


h 

2>f 


/  u'  / 


h 

2Af 


b 

2Af 


/  v' 


b 

21  f 


=  0  otherwise 


(3.6-4) 


For  any  given  light  intensity  l(u,v),  the  light  in¬ 
tensity  observed  through  the  slit  can  be  defined  by 


where 


Io<u-v> 


f+c°  J+Pc  l(u'  ,V  )  g(u-u'  ,v-v'  )  du’dv' 

S  —oc  —oo 

(3.6-5) 


IQ(u,v)  r 
I(u'  ,v«  )  - 
g(u-u'  ,  v-v'  )  = 


(  ‘  f  )P  ' 


observed  intensity  at  (u,v) 

intensity  at  (u*  ,V  ) 

the  filter  function  of  the  slit  cen 
tered  at  u  =  u'  ,  v-v' 

area  of  slit 

"area"  of  slit  in  units  of  spatial 
frequency 


The  slit  produces  an  averaging  effect  which  causes  blur¬ 
ring  by  masking  high -frequency  components  in  the  distribution 
l(u,v)  (see  Fig.  3.6-1).  Clearly,  the  best  way  to  minimize 
blurring  is  to  use  the  smallest  possible  slit  cone is tent  with 
the  sensitivity  of  the  photomultiplier  tube.  In  the  limit  when 
k  -  oo  and  h  and  b  become  zero,  blurring  will  vanish  and 
IQ(u,v)  -  l(u,v).  When, however, the  problem  is  to  estimate  the 
spectral  density  of  a  random  process,  there  are  considerations 
other  than  blurring  that  dictate  a  good  choice  of  scanning 
slit.  When  spectral  density  measurements  are  made  with  the 
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FIG.  3.6-1  ILLUSTRATING  THE  DISTORTING  EFFECT  OF  A  SCANNING  SLIT 

ON  TWO  DIFFERENT  SPECTRA 


-6*S- 


system  of  Fig.  2.4-2  ,  the  light  intensity  in  the  uv 
plane  is  an  estimate  of  W(u,v)  and  hence  subject  to  varia¬ 
bility.  The  stability  of  the  estimate  will  decrease  with  de¬ 
creasing  slit  size  since  the  averaging  takes  place  over  fewer 
spatial-frequency-resolution  cells.  Hence,  from  the  viewpoint 
of  stability,  too  small  a  slit  is  undesirable.  Observe  also 
that  a  small  slit  implies  a  narrow-bandpass  filter.  The  nar¬ 
rower  the  band,  the  greater  is  the  frequency-resolving  capabil- 
ity  of  the  scanner.  However,  the  apparent  spectrum  W(u,v) 
is  an  averaged  version  of  the  true  spectrum  W(u,v),  the  aver¬ 
aging  occurring  over  frequency  intervals  of  the  order  of  1/L 
for  the  "u"  frequencies  and  1  /W  for  the  "v"  frequencies 
[ 1/L  and  1/w  are  approximately  the  half-power  widths  of  the 
function  H(u,v)  defined  by  Eq.  (3.2-2)].  Hence,  it  would  be 
wasteful  to  demand  resolving  capabilities  far  in  excess  of 
l/L  and  1/W  since  this  would  furnish  unstable  estimates. 

The  above  considerations  illustrate  that  stability  may 
be  increased  at  the  expense  of  frequency  resolution  and  vice 
versa.  Since  stability  improves  with  increasing  slit  size, 
the  problem  of  what  constitutes  an  appropriate  upper  bound  on 
slit  size  must  be  considered.  Clearly,  it  is  undesirable  to 
let  the  bandwidth  of  the  slit  be  greater  than  the  smallest 
frequency  intervals  over  which  the  spectral  density  changes 
by  a  significant  amount.  For  a  one-dimensional  process,  con¬ 
sider  a  data  sample  large  enough  to  contain  information  about 

all  correlation  intervals  a  such  that  I  at  I  a  where 

'  '  max 

a  is  the  longest  correlation  distance  of  interest.  The 

max  ’ 

least  distance  over  which  W(u)  will  change  by  a  significant 

amount  is  then  of  the  order  of  a  .  Hence,  the  one-di¬ 

max 

mensional  slit  of  assumed  width  d,  should  satisfy 

d  / 

*f  ^  amax  (3.6-6) 
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By  way  of  illustration,  consider  the  measurement  of 
isotropic  grain  noise  with  a  system  for  which  7  =  6328  8 
(He-Ne  laser),  f  -  1  meter,  and  L  =  W  =  60  mm.  Then  assum¬ 
ing  that  uM  -  fJ>M  -  1000  microns  (which  is  very  large  and 
leaves  room  for  considerable  "clumping"  -  see  Sec.  4.1)  a  rea¬ 
sonable  guide  as  to  the  width  d  would  be  furnished  by 

X  '  d7  (3-6-7) 

max 

which  gives  approximately 

10  microns  <  d  <  600  microns. 

The  problem  of  obtaining  a  stable  estimate  is  best  determined 
by  experimentation.  It  may  turn  out  that  very  stable  estimates 
of  W(u,v)  are  obtained  with  even  a  5-micron  slit.  Or  else 
perhaps  frequency  resolution  in  excess  of  what  was  planned 
may  have  to  be  sacrificed  to  obtain  reasonable  stability. 

The  effect  of  slit  on  the  covariance  may  be  illustrated 
as  follows.  From  Eq.  (3*6-5)  it  follows  that  the  observed 
spectrum  Wq(u,v)  is  related  to  the  apparent  spectrum  W(u,v) 

Wq(u,v)  ~  f+0°  J’4'00  W(u>  ,v'  )  g(u-u'  ,v-v'  )  dudv 

S  —00  —00 

(3*6-8) 


The  Fourier  transform  of  the  above  gives 

0o(a,e)  =  sinc  ff  sinc  ff 

which  shows  how  the  blurring  affects  the  correlation  function. 
Clearly,  if  h  7  <  7f/a  and  b  <f  7  7f/f3  the  blurring 

is  negligible  since  the  sinc  functions  are  almost  constant  over 
regions  for  which  there  is  significant  correlation. 


by 
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4  .  FILM-NOISE  MODELS 


The  experimental  work  in  this  thesis  concerns  the  meas¬ 
urement  of  the  power  spectra  of  film  noise.  Results  are  included 
for  both  |t(x,y)}  and  |R(x,y)|  .  In  this  chapter,  the  pro¬ 

perties  of  existing  and  proposed  film-noise  models  are  discussed. 

4  .1  GENERALITIES 

As  stated  earlier,  in  many  optical  signal- processing  and 
wavefront  reconstruction  systems,  the  fundamental  limitation  to 
high  performance  is  film  noise.  Very  often,  predictions  of  achiev¬ 
able  signal- to-noise  ratios  of  such  systems  take  into  account  film 
grain  noise  only  since  the  phase  noise  can  be  significantly 
reduced  by  using  so-called  liquid  gates  (Thiry  1963 )  .  A  liquid 
gate  consists  cf  two  optically  flat  glass  sections  containing 
between  them  a  liquid  whose  refractive  index  is  chosen  to  be 
equal  to  that  of  the  object  which  is  to  be  immersed  in  the 
gate.  Ideally,  the  optical  path  length  across  any  point  in 
the  aperture  is  then  a  constant  for  all  ray  trajectories  and 
the  phase  fluctuations  resulting  from  random  optical  path  length 
differences  vanish  (Fig.  4.1-1).  When  optical  path  length  dif¬ 
ferences  are  due  to  variations  in  the  refractive  index  of  the 
film  emulsion  or  when  the  object  is  of  a  compound  nature,  that 
is,  made  up  of  several  materials,  the  liquid  gate  can  reduce 
but  not  eliminate  phase  errors.  The  relative  contributions 
of  film  thickness  irregularities  and  variations  of  refractive 
index  to  the  phase  deterioration  of  a  quasi-coherent  beam  has 
been  measured  for  photographic  film,  for  which  it  has  been 
shown  that  thickness  irregularities  are  the  chief  sources  of 
phase  noise  (Ingalls  i960). 
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OBJECT  WITH 


PHASEFRONT  OF  COLLIMATED  DETERIORATED  PHASEFRONT 

COHERENT  WAVE  DUE  TO  THICKNESS  VARIATIONS 


LIQUID  WITH 


GLASS  CONTAINER 


A-0Q6-I-00  -006S 

FIG.  4.1-1  ILLUSTRATING  THE  ACTION  OF  A  LIQUID  GATE  IN  REDUCING  NOISE 
DUE  TO  RANDOM  THICKNESS  VARIATIONS 


As  a  measure  <  what  constitutes  acceptable  variations 
in  the  optical  path  length,  that  is, variations  not  sufficiently 
large  to  deteriorate  the  phase  front  of  a  coherent  beam  sig¬ 
nificantly,  the  Rayleigh  limit  may  be  employed.  This  criterion 
specifies  that  optical  path  differences  no  more  than  one-quarter 
of  the  wavelength  of  light  are  permissible  over  the  active  aper¬ 
ture.  Hence,  for  film*  immersed  in  air  and  coherently  illumni- 
ated  in  an  optical  system  employing  a  He-Ne  laser  (bj28  a),  the 
maximum  allowable  variation  in  film  thickness  is 


ATmax  4(nf-l ' 


032 -mm 


which  is  a  very  small  variation. 

When  film  with  refractive  index  n2  is  immersed  in  a 
liquid  gate  with  refractive  index  n2 ,  the  Rayleigh  limit  is 
given  by 

At  -  - ...  ■■■  .. 

max  4  |  n2_  n J 


A  plot  of  the  Rayleigh  allowable  film-thickness  variations  vs 
immersion -liquid  index,  for  A  =  6328  A  and  film  with  refrac¬ 
tive  index  of  1.5,  is  given  in  Fig.  4.1-2. 

In  the  next  section  several  models  of  film  granularity 
are  discussed.  These  models  are  highly  simplified  and  their 
use  in  signal- to-noise  predictions  is  at  best  a  first-order 
approximation  to  the  problem.  Real  grain  structure  is  very 
complicated.  Photomicrographs  have  shown  that  developed  grains 
rarely  have  the  shape  of  the  undeveloped  grains, which  are  halide 
crystals  of  definite  and  regular  shape.  The  developed  silver 
grains  are  either  twisted  filaments  of  silver  or,  usually, 
masses  of  such  filaments.  Even  for  uniformly  exposed  film, 
the  density  of  grains  is  not  uniform  throughout  the  emulsion, 
the  number  of  developed  grains  being  most  numerous  near  the 

*  A  refractive  index  of  1.5  for  film  is  assumed. 
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ALLOWABLE  FILM  THICKNESS  VARIATION  (MMxlO 


10 


REFRACTIVE  INDEX  OF  FILM  =  1.5 

X  x  6328  8 


_ A -006- 1 -00  -  0066 

FIG.  4.1-2  MAXIMUM  ALLOWABLE  VARIATION  OF  FILM  THICKNESS  (RAYLEIGH 

LIMIT  CRITERION  ) 
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outer  surface  of  the  emulsion  and  decreasing  with  increasing 
displacement  from  the  surface.  The  distribution  in  depth  de¬ 
pends  on  such  features  as  the  concentration  of  the  silver 
halide,  the  size  of  the  grains,  the  thickness  of  the  emulsion 
and  the  composition  of  the  developer  (Mees  &  James  19 66). 


Formerly,  it  had  been  considered  that  the  grain  is  the 
unit  of  photographic  action  and  that  a  grain  develops  inde¬ 
pendently  depending  on  its  quantum  efficiency  (the  probability, 
e,  that  a  photon  incident  upon  the  grain  is  absorbed  by  that 
grain)  and  on  the  minimum  number  of  photons  m  required  to 
make  a  grain  developable.  Subject  to  the  assumption  of  Poisson- 
distributed  photon  arrivals,  the  probability  that  a  given  grain 
with  parameter  (f,m)  becomes  opaque  during  development  is  then 
given  by  (Silberstein  1922,  Goodman  1967 ) 


P(t  =  0)  =  £ 
k=m 


(e  /Ka 0 


■LVs 


(M.1-1) 


where 

E  =  total  exposure  (energy  density)  to 

T  which  the  grain  is  subjected 

-Tf  =  Planck' s  constant  divided  by  27t 
f  =  frequency  of  the  radiation 
O)  =  27Tf 

Ag  =  the  area  of  the  grain 

However,  as  early  as  1922  it  was  found  that  certain 
grains  develop  by  virtue  of  their  contact  with  developed  grains 
and  that  the  process  may  repeat  several  times,  ultimately  form¬ 
ing  a  clump  whose  size  may  be  many  times  the  size  of  a  single 
grain.  This  phenomenon  is  called  clumping  and  may  be  aggrav¬ 
ated  by  certain  developers,  notably  those  that  contain  hydrazine 
(Trivelli,  et.al.  1922). 


If  grains  were  to  develop  independently,  then  the  co- 
variance  function  would  be  expected  to  be  zero  for  correlation 
intervals  greater  than  the  length  of  the  longest  grains.  When 
clumping  is  present,  this  is  not  the  case  however;  the  clump¬ 
ing  may  be  viewed  as  a  random  process  with  a  Markov- type  compon¬ 
ent,  namely  the  state  of  a  grain  depending  in  some  degree  on 
the  state  (developed  or  undeveloped)  of  an  adjacent  grain.  When 
clumping  is  widespread,  the  length  of  significant  correlation 
intervals  increases  the  amplitude  of  the  low-frequency  components 
of  the  spectral  density. 

The  author  has  no  knowledge  of  stochastic  grain  models 
designed  to  include  the  statistics  of  clumping. 

4 .2  MODELS  OF  PHOTOGRAPHIC  GRANULARITY 

A .  The  Random  Checkerboard  Model 

The  random  checkerboard  model  is  shown  in  Fig.  4.2-1. 

The  transparency  is  assumed  to  be  composed  of  square  cells, 
each  of  length  t  and  arranged  in  parallel  rows  and  columns. 

A  given  cell  is  assumed  to  be  either  perfectly  opaque,  corres¬ 
ponding  to  a  grain  that  has  been  developed,  or  perfectly  trans¬ 
parent,  corresponding  to  a  grain  that  has  remained  undeveloped. 
Whether  or  not  a  grain  has  been  developed  depends  on  (1)  the 
sensitivity  of  the  grain,  (2)  the  energy  incident  on  the  grain 
during  exposure,  and  (3)  the  degree  to  which  the  photographic 
plate  has  been  processed. 

Let  the  transparency  be  uniformly  exposed.  In  the  ran¬ 
dom  checkerboard  model,  the  grains  are  assumed  to  develop  in¬ 
dependently  so  that  the  probability  that  a  given  grain  is  opaque 
is  independent  of  the  state  of  adjacent  grains.  The  amplitude 
transmittance  of  the  transparency  may  be  written  as 

t(x,y)  =  tQ  +  T(x,y)  tQ  =  < t (x,y )> 

(4.2-1) 
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FIG.  4.2-1  RANDOM  CHECKERBOARD  MODEL  OF  FILM  GRAINS 
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and  a  typical  sample  function  is  shown  in  Fig.  4.2-2.  The 
process  |t(x,-y)|  is  a  two-dimensional  analogue  of  the  random 
binary  signal  wave  occurring  in  communication  theory  in  which 
the  starting  time  of  a  pulse  is  a  random  variable,  uniformly 
distributed  over  the  duration  of  the  pulse*.  The  duration  of 
transmittance  pulses  must  be  multiples  of  f  .  The  correla¬ 
tion  function  of  the  process  jT(x,y)|  can  be  readily  shown 
to  be 

>,f-'  =-  te(i  -  to)  (i  _  laJ-  )  (i  .  JaL  )  |a|  <  , 

lei  7  f 

0T(a,f)  0  otherwise  (4.2-2) 


and  the  associated  power  spectrum  Wt(u,v) 


is  given  by 


WT(u,v) 


(sin  7tv^  \ 

m  )  (4.2-3) 


where  u  and  v  are  spatial  frequencies.  For  .u  and  v 
fixed,  W  (u,v)  has  a  maximum  at  t  =  .5.  At  frequencies 
given  by  u  n/-f  and  v  -  Y/l  where  n  =  1,  2,  .  .  . , 
and  k  1,  2,  .  .  Wt(u,v)  is  zero. 

From  Eq.  (4.2-2)  it  can  be  seen  that  0  (a,(3)  is  an 
even  function  of  both  its  arguments.  However,  since  the  co- 
variance  function  is  not  a  function  of  (a2  +  |32)?>  the  pro¬ 
cess  is  not  isotropic.  Although  the  random  checkerboard  model 
is  not  a  very  realistic  model  of  granularity,  it  has  been  used 
in  signal-to-noise  predictions  for  holographic  systems  (Good¬ 
man,  1967)  and  also  in  entropy  and  granularity  studies  (O'Neill, 
1958).  The  random  checkerboard  model  has  also  been  proposed 
as  a  standard  of  granularity,  against  which  the  granularity 
of  real  films  could  be  compared  (Fry,  1963). 


*  See  for  example  Schwarz  and  Friedland, I965. 
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'  FIG. 4.2*2  THEORETICAL  POINT  APERTURE  MICRODENSITOMETER  TRACE  OF  RANDOM 

CHECKERBOARD  MODEL  OF  FILM  GRAINS 


B. 


THE  OVERLAPPING  CIRCULAR  GRAIN  MODEL  (OCG) 


The  initial  treatment  of  the  transmittance  and  auto¬ 
correlation  function  of  this  model  was  given  by  Picinbuno 
(Picinbono  1955)  who,  in  deriving  the  autocorrelation  func¬ 
tion,  assumed  that  an  exposed  grain  is  completely  opaque  while 
the  unexposed  grain  and  the  emulsion  were  completely  trans¬ 
parent.  Savelli  (Savelli  1958)  extended  Picinbono' s  treatment 
to  allow  for  greater  than  zero  transmittances  for  exposed  grains 
and  less  than  unity  transmittances  for  the  unexposed  grains. 
However,  since  an  exposed  grain  is  essentially  a  particle  of 
silver  which  is  opaque  in  the  visible  light  region,  it  is  rea¬ 
sonable  to  assume  that  the  transmittance  of  an  exposed  grain 
is  always  zero . 


The  basic  OCG  model  is  based  on  the  following  assump¬ 
tions:  (l)  the  grains  are  circular  and  the  centers  of  these 

circles  fall  at  random  on  a  plane  as  independent  events; 

(?)  the  grain  centers  are  Poisson  distributed  in  the  sense 
that  the  probability  that  there  are  n  grains  in  an  area  A 
is  given  by 


P 


n 


(A) 


n! 


(4 .2-4) 


where  n  is  the  expected  number  of  grains  in  A?  and  (3)  the 
grain  sizes  are  random  variables  with  a  distribution  law 

P(* )  =  Prob  (b  <  j ) 

where  L  is  a  random  variable  denoting  the  grain  radius.  The 
factors  governing  grain  development  are  the  same  as  in  the  ran¬ 
dom  checkerboard  model . 

Subject  to  assumptions  (1)  through  (2),  it  can  be  shown 
(Picinbono  1955)  that 
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t  (x  ,y  V'  = 


(4.2-5) 


e"^d  f  [1  -  P(0]W 

0 


where  d  average  number  of  grains  per  unit  area.  The  correla¬ 
tion  function  for  this  model,  can  be  shown  to  equal 

(Picinbono  ly55) 

0t(a»P)  =  </t(x,y)t(x  +a,y  +p  )>  =  4>t(p) ,  where 

P  =  (a2  +  f2)r 

=  exp  Ud  (tJP/2  (1  -  P(0)C  dC 
L  0 

+  r  (7T  -  cos-1  p/2y)  (1  -  P (r))r  d")l 
P/2  J 

(4.2-6) 

By  Eq.  (2.2-6)  it  can  be  seen  that  the  OCG  model  has  isotropic 
statistics.  In  the  special  case  where  all  the  grain  are  assumed 
to  have  the  same  radius  t  ,  there  results  (Fig.  4-2-  ]) 


P(-f)  =  0 

p(*0  =  i 


*  <  fo 

*  1  ;o 


(4.2-7) 


so  that 


-Tiin.2 


t  =  e 
o 


o^d 


(4.2-8) 


The  autocorrelation  function  for  this  special  case  is 

-27 -J)  2d 


\(p)  =  e  ° 


=  t 


P  ">  2. 


„  ? -f  ed  F  ( p/2f  ) 

0t(p)  =  V  e  °  IP/ 


(4.2-9) 


where 


F(p/2^0)  =  2/tt  [(cos-1  (p/2f0)  -  p/2'o  -  (p/2»0)2  1 

(4  .2-10) 


P  < 
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FIG.  4.2-3  OVERLAPPING  CIRCULAR-GRAIN  MODEL  OF  PHOTOGRAPHIC  FILM. 


The  covariance,  0  (a,p),  is  given  by 


0T( a, 3)  =  <t(x,y)  t  (x-,A,y+fl)>=  0T(p) 


and 


♦T(P)  -  V  rt0'F(p/s<o)  -11 


=  0  otherwise 


(4  ,2-11) 


0  <  p  <  2^q  (4 .2-3.2  ) 


When  the  unexposed  portions  of  the  film  have  a  trans¬ 
mittance  h,  where  h  is  not  necessarily  unity,  then  the  pre¬ 
vious  results  are  modified  to 


-rr^  2d 

t.  =  <t(x,y)  >  =  h  e  c 


(4.2-13) 


and 


\ 


0-(p)  = 


■  v[ft) 


h\  p(P/2'o) 


‘] 


wl  re  F(p/2.f  )  is  as  defined  earlier 


(4  .2-14) 


The  power  spectra  w(u,v)  of  this  model  is  given  by 

w(u,v)  =  ^  r  Up)  e-i2™Pcos(n-p)p  dpd0 
0  0 


=  2v  f  0_(p)  j  (27iwp)p  dp 
0  t  ° 


(4.2-15) 


=  W  (w ) 

where  w  =  (u2  +  v2)^. 


-80- 


The  normalized  power  spectra,  defined  by 


W'n(w)  f - W-l-  V  (4.2-16) 

8''o  fco 

has  been  numerically  integrated  and  is  shown  in  Fig.  4.2-4 

as  a  function  of  the  dimensionless  variable  X  -  wj  with 

o 

tQ/h  as  a  parameter.  The  Picinbono  model  is  obtained  when 
h  =  1. 


This  model,  like  the  random  checkerboard  model,  has 
a  covariance  function  which  is  zero  for  intervals  greater  than 
the  grain  diameter.  The  relative  power  spectra,  at  dc,  of  the 
checkerboard  and  OCG  model  are  shown  in  Fig.  4.2-5,*  as  a 
function  of  amplitude  transmittance. 

4.3  MODELS  FOR  TOTAL  FILM  NOISE 

The  complex  transmittance  of  uniformly  exposed  film  is 
a  stochastic  quantity  that  describes  the  joint  effect  on  co¬ 
herent  light  of  film  granularity  and  variations  in  the  optical 
path  length  through  the  cross  section  of  the  film.  In  writing 
for  the  complex  transmittance  at  a  point  (x,y) 


R(x,y) 


t(x,y) 


ejV  (x,v) 


<t  (x,y  )^>  =  tQ 
<fR(x,y)>  =  Rq 

(4.3-1) 


the  amplitude  transmittance  t(x,y)  is  primarily  a  conse¬ 
quence  of  the  granularity  of  the  film  emulsion,  and  ip(x,y)  is 
due  to  non-uniformities  in  the  thickness  of  the  film  and  the 
refractive  index  of  the  emulsion.  Attempts  to  assess  the  rela¬ 
tive  contribution  of  these  two  factors  to  the  phase  noise  pro¬ 
cess  {V'(x,y)i  have  shown  that  the  major  cause  of  phase  deter¬ 
ioration  of  coherent  light  as  it  passes  through  film  is  the 
random  variation  in  film  thickness  (Ingalls  i960 ) . 


*  No  significance  should  be  attached  to  the  difference  in  heights 
between  the  two  curves. 
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NORMALIZED  POWER  SPECTRUM 


l0  >  GRAIN  RADIUS 
10  *  AVERAGE  TRANSMITTANCE 
h  >  TRANSMITTANCE  OF  GELATIN 
X  >  DIMENTIONLESS  VARIABLE 


A-006-  I  -  00*  0070 


FIG.  4.2-4 THE  NORMALIZED  POWER  SPECTRUM  OF  GRAIN  NOISE  AS  COMPUTED 
FOR  THE  OVERLAPPING  CIRCULAR -GRAIN  MODEL.  (  ADAPTED  FROM 

SAVELLI  1958  ) 
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FIG. 4.2-5 VARIATION  OF  POWER  SPECTRUM  AT  ZERO  SPATIAL  FREQUENCY  WITH 
AVERAGE  AMPLITUDE  TRANSMITTANCE .  (ADAPTED  FROM  GOODMAN  1967) 
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In  attempting  to  construct  stochastic  models  for  the 
complex  transmittance  of  uniformly  exposed  film, . it  is  of 
primary  importance  to  determine  whether  |t(x,y)f  and 
|V/(x>y)|  can  be  considered  as  independent  processes.  Al¬ 
though  it  would  appear  that  surface  corrugations  of  the  emul¬ 
sion  have  little  to  do  with  film  amplitude  transmittance,  there 
is  in  fact  evidence  to  suggest  that  these  are  not  independent 
phenomena.  For  example,  it  has  been  found  that  when  developed 
film  has  been  immersed  in  an  aqueous  solution,  the  more  highly 
dense  (low  transmittance)  areas  dry  more  quickly  than  regions 
of  low  density  (high  transmittance)  and  thus,  while  drying, 
pull  the  moist  and  flexible  low-density  regions  with  them. 

The  resulting  shift  in  the  emulsion  which  is  both  parallel  and 
normal  to  the  plane  of  the  emulsion  is  called  the  gelatin  effect . 
The  changes  in  the  thickness  of  the  emulsion  layer  are  visible  in 
relief  on  the  surface  (James  and  Higgins  I960).  Thus,  the  re- 
liefing  is  a  function  of  transmittance  and  in  dealing  with  de¬ 
terministic  transmittance  variations  such  as  images,  the  phase 
and  amplitude  transmittance  cannot  be  adjudged  to  be  independ¬ 
ent  quantities.  It  has  been  suggested  that  the  phase  modula¬ 
tion  that  results  from  reliefing  can  be  described  by  a  complex 
transmittance  with  a  phase  dependence  that  is  linearly  propor¬ 
tional  to  the  relative  transmittance  at  that  point  (Leith  1962). 

When  considering  the  random  properties  of  uniformly  ex¬ 
posed  film,  however ,  the  average  deterministic  reliefing  is  zero 
and  inhomogeneities  in  film  thickness  are  probably  more  likely 
a  reflection  of  the  film  manufacturing  process.  In  any  case, 
the  reliefing  at  any  point  on  the  film  is  the  sum  of  the  pull¬ 
ing  forces  of  many  clumps  of  grains,  randomly  situated  and  ex¬ 
erting  their  forces  in  random  directions.  For  these  reasons, 
it  may  be  reasonable  to  assume  a  simplified  model  of  the  com¬ 
plex  transmittance  of  uniformly  exposed  film  consisting  of  the 
independent  processes  |t(x,y)[  and  |V'(x*y)|  with  |^(x,y)| 
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a  wide-sense  stationary,  Gaussian  process.  The  assumption  of 
a  Gaussian  process  allows,  in  principle  at  least,  the  covari¬ 
ance  of  the  process  |$(x,y)}  =  to  be  computed  in 

terms  of  the  moments  of  |^(x,y)|.  This  follows  from  the  fact 
that  the  density  functions  of  all  orders  are  determined  by  the 
covariances  of  all  pairs  of  random  variables  and  no  higher-order 
moments  are  needed.  Under  these  assumptions,  the  correlation 
function  of  the  process  |R(x,y)J  is  given  by 

*R(a,f)  =  r.‘>T(a,e)  +  to?l  e"(°2  "  (4.3-2) 

and  -  1 ^  ( 0,0).  The  covariance  function  of  |r(x,y)J  is 
given  by 

.\(a,f)  =  4R(a,P)  -  |R0I  ?  (4.3-3) 

-  (  1  /P  )  n  2 

where  Rq  =  <^R(x,y)v  -  tQe  v  ’  ,  a  result  that  follows 

from  the  formula  for  the  characteristic  function  of  a  Gaussian 
variate*.  When  the  phase  errors  are  small,  such  as  they  might 
be  when  the  liquid  gate  is  slightly  mismatched  to  the  film,  the 
exponential  can  be  expanded  into  a  Taylor  series: 

-  (o?  -  b  (a,P  ) ) 

e  ~  1  -  os  +  2y(a,P  )  o2  1 

(4.3-4) 

and  the  diffraction  system  of  Fig.  2.4-1  will  produce  an  aver¬ 
age  relative  light  intensity  in  the  (u,v)  plane  given  by 

^l(u,v)  S  =  tQ2(l  ”  °2>)  LW  sine2  uL  sine2  vW  +  Wr(u,v) 

(4.3-5) 

where  it  has  been  assumed  that  all  significant  correlation  in¬ 
tervals  are  very  small  compared  with  the  dimensions  L,  W. 

*  See  for  example  Schwarz  and  Friedland , I965. 
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Since  g2  is  small,  the  amount  of  light  energy  scattered  away 
from  the  zero-order  light  through  phase  errors  is  small  and 
the  attenuation  is  essentially  caused  by  the  factor  t  .  When 
the  phase  angles  are  large,  g2  sx  1,  the  average  light  in¬ 
tensity  is  given  by 

_n2 

/l(u,v)N  T  t  2  e~  LW  sine2  uL  sinc^  vW  +  W  (u.v'i 
o  r 

(4.3-6) 

Clearly  in  this  case,  a  great  deal  of  energy  has  been  scat¬ 
tered  away  from  the  zero  order  into  the  power  spectrum  of  the 
complex  process  jr(x,y)|. 

Rowe  has  dealt  with  the  problem  of  phasefront  deter¬ 
ioration  of  coherent  light  by  lenses  with  large  random  imper¬ 
fections  (Rowe  1966).  The  imperfections  are  assumed  to  create 
phase  errors  with  stationary,  isotropic  Gaussian  statistics. 
This  problem  is  similar  to  the  problem  of  random  aperture  er¬ 
rors  on  antenna  radiation  patterns  first  treated  by  Ruze 
(Ruze  1952).  One  would  expect  the  average  back- focal 
plane  light  intensity  distribution  in  a  coherent  optical  sys¬ 
tem  to  be  very  similar  in  the  two  cases  of  phasefront  deter¬ 
ioration  by  film  followed  by  a  very  good  lens  or  phasefront 
deterioration  by  a  poor  lens  alone  if  the  statistics  in  the 
two  cases  were  the  same.  Rowe  considered  typical  examples  of 
both  differentiable  and  non-dif ferentiable  processes.  If  the 
covariance  is  assumed  to  be  of  a  type 

*„,(p)  -  p  =  (a2*  Pz)l/s 

*  (4.3-7) 

corresponding  to  a  differentiable  random  process,  the  asso¬ 
ciated  power  spectrum,  W^(w),  is 

W  (w)  =  g2  (tt/c2)  g*^7™/0)  w  =  (u2  +  v2)1/2 

(4.3-8) 
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and  if  the  covariance  is  of  the  type 


0^(i  )  =  e  c  IN 


(4.3-9) 


corresponding  to  a  non-dif ferentiable  process,  the  associated 


power  spectrum  is 


W  (w)  -  - 2%C" 

j-]_  +  /2™ 


(4.3-10) 


The  covariance  of  the  zero-mean  process  0(x,y)  associated 
with  the  process  6(x,y)  is 


^  (a ,(3  )  -  /  [^(x+a,y^p  )  -  *(x,y)]  >  _ 


(4.3-H) 


If  the  process  |^(x,y)}  has  a  covariance  given  by 
Eq.  (4.3-7),  the  power  spectrum  of  |o(x,y)|  is  approximately 
given  by  ( ,-2  '  '  1  ) 


Wt.(w)  -  > 


-  (2;iw /a  '  )2 


(4.3-12) 


while  if  the  process*  |^(x,y)f  has  a  covariance  given  by 
tq.  (4.3-9),  the  power  spectrum  of  {0(x,y)}  is  approxi¬ 
mately  given  by  - 

Vw)  =  l\  ./2^\213  p  (4.3-13) 


[i 


In  Eqs.  (4.3-12)  and  (4 . 3-13 ) : 


•2  =  7 


d0(x 


I  2  +  I 


<30(x , 


I  2  > 


m  =  the  initial  ^slope  of  0^  (p) 
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c  =  ^  for  Eq.  (4.3-7 ) 
c  m  for  Eq .  (4 .3-8 ) 


The  power  spectrum  of  the  random  process  |r)(x,y)}  gather 

case,  subject  only  to  independence  between  jt(x)J  and  {  (x)f 
is 

2 

W  (w)  =  e“°  W  (w)  +  t 2  Wfl(w)  +  W  (w)*Wfl(w) 
r  T  °  9  T  0  (4.3-14) 

It  is  to  be  emphasized  that  the  spectra  W^(w)  given  by 
Eqs .  (4.3-12)  and  (4.3-13)  are  approximate.  Truncated  ver¬ 
sions  of  0  (p)  (p  small)  were  used  in  their  computation. 
Nevertheless  these  spectra  serve  as  a  first  approximation  to 
the  problem. 
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VARIANCE  OF  POWER  SPECTRUM  ESTIMATE 


r 


In  this  chapter,  the  variance  of  the  power  spectral  den¬ 
sity  estimate  is  derived.  A  real,  wide-sense  stationary,  Gaussian 
process  is  assumed.  The  estimate  occurs  in  the  form  of  a  light 
intensity  in  the  (u,v)  plane  of  the  two-dimensional  optical  sys¬ 
tem  shown  in  Fig.  2.4-2  .  The  assumption  of  a  Gaussian  law 
when  the  process  is  a  random  transmittance,  as  it  must  be  when 
spectral  density  estimates  are  made  by  diffraction,  is  somewhat 
unrealistic  since  by  definition, transmittances  must  be  con¬ 
fined  to  values  greater  than  zero  and  less  than  unity.  However, 
in  the  case  of  film  grain  noise,  microdensitometer  traces  bear 
out  the  fact  that  the  transmittance  deviations  about  the  mean 
are  in  good  agreement  with  the  Gaussian  probability  density 
[Mees  and  James  (5^  1966]. 

The  results  derived  in  Sec.  5.2  are  exact  for  rectangu¬ 
lar  or  diamond- shaped  apertures.  The  form  of  the  solution  when 
the  aperture  is  large  suggests  that  the  results  may  be  used  as 
guides  for  other  apertures  with  comparable  resolution  and  area. 

5 . 1  FORMULA  FOR  THE  VARIANCE 

Consider  the  Gaussian  process  |t(x,y)|  with  mean  value 
t  for  which  the  estimated  spectrum  is  •  (u,v).  Then  for  a 

rectangular  aperture  of  area  A  -  LW,  the  estimated  apparent 
spectrum  is 


•  (u,v)  =  I ( u , V  )  -  to2  H (u, V  ) 


and  the  mean  and  variance  are  given  by 


(5.1-1) 
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H(u,v) 


<  ‘Vf  ( u ,  v  )  >  =  W(u,v  )  =  w(u,  v  )  * 

(5.1-2) 

var  (^f  ( u  ,  v  )  )  =  W2 ( u, v )  +  F2(u,v)  +  2tQ2  H(u,v)  [F(u,v) 

-t  W(u,v  )  ]  (5.1-3) 


where 

total  relative  light  intensity 

system  function  for  a  rectangular 
aperture 

LW  sinc2uL  sinc2vW 

+00  +00 

LW  /  /  w(f,Tj)  sine  (u+f  )L  sine  (u-f)L 

—00  —00 

.  sine  (v-tj)W  sine  (v+rj)w  df  dT} 
w(UjV')  apparent  or  biased  spectrum 

-  f+n  j'4"  Wf?,’i)  sine  (u-“)l  sine  (v-])w  dfd| 

—  C  Y  —  CX 

5.2  DERIVATION  OF  EQUATION  (5.1-3) 

The  relative  light  intensity  l(u,v)  is  defined  by 


I (u, v  )  = 
H(u,v)  = 


F  ( u ,  v  )  - 


I(u,v) 


_L  , L/2  rL/2 

LW  ■'  /  J  / 

-L/2  -L/2 


rW/2  , Vi/2  t  t  e-  j27r(u(x1-x2  )-fv(yi-y2  )  ) 
-W/2  -W/2  1 

(5.2-1) 


.  dx:dx2dy1dy2 


where 

t.  -  t(x.,y.) 

i  =  1,2, 

The  square  of  l(u,v)  is 
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e-  j27r(u(x1-x2+x3-x4  )+v(yi-y2+y3-y4  ) 


I2(U,V)  -  (Lw)2  /•••/  tit2t3t4 


8  times 


(5.2-2) 


II  dx .  dy . 
1  1  i 

i=l 


To  evaluate  <I2(u,v)  )>  it  is  required  to  obtain  Ot1t2t3t4)>. 
Observe  that 


^1^2^3^4 


4  3  4 

t  4  +  t  3  F.  T.  +  t  2  V  £  T  T 

°  °  i=l  1  °  i=l  j>i  1  3 


+  to  ill  j>i  k>  j  Ti^Tk  +  TlTeT3T4 


(5.2-3) 


The  Tj‘s  are  Gaussian  variates  with  zero  mean.  Hence 


<'TlT2T3T4>^^34  +^3024  + 

<TiTiTk  >  =  0 


2  3 


(5.2-4) 


(See  for  example.  Miller  1964.) 

In  Eq.  (5.2-4)  0i;j  £  ^x(xi,yi)  t  (x..,y..)  >.  It  follows  that 


<tlt2t3t4  >  =  t4  +  V  (012  +  013  +  014  +  023 


+  02  4  +  0  3  4  )  +  012034  ■*"  013024  01402  3 


(5.2-5) 
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It  is  convenient  to  write 


Ir (u.v  ) 


V 


u 

7  I. 
i-  1  ' 


ir:l 


?-r,) 


where  the  1^0  s  are  t^ie  contribution3  of  integrals  involving 
the  covariances  and  I:  is  the  contribution  of  the  tQ4 

term.  Consider  the  contribution  of  I31,  I32,  and  I33.  For 
I31  there  results 


ls,  =  (£)*  /.../  e-i2;r<u<xi-x2+x3-x4)  +  v(y1-ye+y.)-y4) 


8  times 


(5-2-7) 


.  n  dx.dy. 
i=l  1  1 


4  times 


/I  %  r  ,  ,  -j2/r[u(x3-x  4)  +  v(y  -y  )  ] 

!nj )  /---J  e  '  4  dx3dx4dy ady4 


4  times 


By  Eq.  (2.4-14),  each  integral  is  I  1  (u ,  v )  /[(j/>.2  f20~  W(u,v). 

o  II 

Hence  I31  =  W  (u,v).  The  contribution  of  I33  which  involves 
$ 14023  is  also  W2(u,v).  For  I32,  there  results 


T  /  1  v  f  ,  ,,  -  j2Tr[u(xa+x3 )  +  v(yj4  y3 )  ] 

132  =  ( LW '  '13  e 

4  times 


dx1dx3dy1dy. 


/-i_\  r  r  ,  Qj27r[u(x  +x4)  -!-v(y  +y  )] 

•  'uj)  J---J  't‘s- »  e  dXsdx4dy2dy4 


4  times 


(5.2-8) 
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which  may  be  written  as 


I  —  T I  t" 
■*■32  A  32 


(5.2-9) 


where  the  primed  quantities  are  the  two  integrals  involved  in 
the  product.  Consider  I'32,  given  by 


!•*  -  £  /V2  /V2  /"(2  i"?  y,-y.) 


-L/2  -L/2  -W/2  -W/2 


-  j27r[u(x1+x3  )  +  v(yj+y3)] 
e  dx jdx3dy jdy3 

(5.2-10) 


This  integral  can  be  simplified  by  considering  the  transformation 


a 

P 

7 

6 


X1  "  X3 


yx  -  y3 

+  x3 
Vi  +  y3 


(5.2.11) 


Q  ^  ^  6 

which  has  a  Jacobian  J(~  )  =  4.  The  field  of  inte- 

Xl'X3'"ljy3 

gration  is  shown  in  Fig.  5.2-1.  It  is  bounded  by 


(L+a  )  <  7  <  L  +  a 

(L-a  )  <  7  <  L  -  a 

■  (W+p  )  <  6  <  W  +  p 

(W-P)  <  6  <  W  -  3 


-L  a  <_  0 

L  >  a  >  0 

-W  <_  p  <  0 

w  >  p  >  0 


(5-2.12) 


The  result  of  using  this  transformation  is 
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FIG. 5.2- 1  DIAGRAM  SHOWING  THE  FIELD  OF  INTEGRATION  FOR  EO.  (5.2-10) 

WITH  THE  TRANSFORMATION  OP  EO.  (5.2-11 ) 
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32  u  ' LW 


z  iffe'  fL  fw  0(“.6) 


sin  2ttu(L  -  |a|  ) 


-L  -W 


7TU 


(5.2-13) 


s  in 


27;v(W  -let) 


7FV 


da  dp 


Since  I r  ,j  is  an  even  function  of  u  and  v,  I "  3p  =  I '  ^ 

F  ( u  ,  v  )  and  I  32  =  F2(  u,v).  The  function  F(u,v)  can  be 
written  as  follows 


F(u,v)  -  i  (ij) 


—  oc  —  oc 


/  dfdr,W((  .r,)/ 

-L 


L  sin  2ttu(l  -la|)  ej2/T^ada 

7TU 


/W 

-w 


sm 


2  v(W  -  iBl  )  j2TP,p 


/TV 


dP 


(5.2-14) 


where  use  was  made  of  the  relation 

H-x, i0  =■-  r  f°°  w(f,i)  ej27r(fa+Tifi)  dfdn.  (5.2-15) 

—  00  —  oc 

The  integrals 


f  _  ( U  ,  f  )  -  -  -i— 

/L 

sin  2?ru  (L  - 

lot  1  ) 

J27Tta  da 

-L 

(5-2-16) 

fTI(v,rl)  - 
W  V  1  '  7TV 

/W 

-W 

sin  2ttv  (w  - 

IPI  ) 

e  j2TT11&  jp 

can  be  written,  after  some  tedious  but  straightforward  algebra,  as 
fT(u,f)  -  2L2  sine  (u+f  )L  sine  (u-£)L 

L  (5.2-17) 

f^(v»rl)  -  2W2  sine  (v+i|)w  sine  (v-rj)W. 
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Hence 


Fixu,v^  -  LW  J'00  W(f  ,^,)  sine  ( utf  )L  sine  (u-t  )L 

—  co  —no 


.  sin  (v-r,  )W  sine  (v-4  ,)w  cH  cl  ( 


( l . 2-18 } 


Observe  that,  when  L  is  large,  the  product  sinc(u,  £'L  sinc(u-t;'L 
is  small  for  any  value  of  u  except  those  near  the  origin  'see 
Fig.  5.2-2  )  ,  because  the  overlap  occurs  in  the  "tails"  of 
the  sine  functions.  When  L  becomes  very  large,  (a  similar 
argument  holds  for  W)  the  amount  of  overlap  becomes  a  vanish¬ 
ingly  small  and 

lim  F(u,v)  -  0  u,v  not  zero  (5.2-19) 

WjL-.; 

The  exception  occurs  at  the  or  .gin;  there 

F(0,0)  -  LWj  °°  j"  W(f,’|)  sine2  fL  sine2  ijW  d^dij 

(5.2-20) 

for  which 

lim  F ( 0, 0 )  =  W(0,0)  (5.2-21) 

W  t  L-«oo 

The  total  contribution  of  the  last  three  terms  in  Eq.  (5.2-2) 
is  2W2(u,v)  -f  F2(u,v). 

Consider  now  the  second  summation  in  Eq.  (5.2-6). 

There  are  two  types  of  terms  -  those  where  the  arguments  of 
the  covariance  and  the  exponents  contain  the  difference  of  the 
variables  of  integration  and  those  where  the  exponent  contains 
the  sum  while  the  covariance  contains  the  difference.  There 
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are  four  terms  of  the  first  type  and  two  of  the  second.  Each 
of  the  four  terms  where  the  arguments  agree,  contributes  an 
amount 


*«i  = 


(LW) 


+■  2 

O  r  r  v  “i27Tru(x,-X„  )  +  V(V.-V„)1 

J7T2-  /.../  0(x1-x2,y1-y  )  e  J  L  v  i  2  ;  x  y2  n 


4  times 


dx1dx2dy1dy2 .1/ 


L/2  e~  j27rux  dx 

L/2 


fW/2  e-K-HVV  dy 
'-w/2 


(5.2-22) 


*2  3  *2  4  *2  6 


to2  w(u,v)  H(u,v). 


The  two  remaining  terms  involve  the  covariances  <£±3  and 
0  .  whose  arguments  do  not  agree  with  their  respective  expon- 
ents .  Thus, 

I22  =  I25  <Hx1-x„Vl-V,)  e-j2x[u(x1+x,)+v(y1+y3)l 

4  times 


.  dx  dx  dy  dy  /L/2  e":’27TUX 

-L/2 


tQ2F(u,v)  H(u, v ) ]  / 


dx 


2  ■  ,W  2  ■  ^  *  2 

/  -  32  nvy 

-W  2 


dy 


( b .2-23 ) 
(5.2-24) 


I  y  which  involves  only  the  constant  tQ4  ,  is  given  by 

Zl  =  to4  H2(u>v)-  (5-2-25) 


/ 


/ 
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Summing  all  the  contributions  gives 


I2 (u,v ) 


H (u, v )  +  W  (u,  v  )  )2  4  W2 ( u , v )  +  F2(u,v) 


4  2t  2  H  (u  ,  v  )  [  F (u , v  )  +  W  (u ,  v  )  ] 


(5.2-2  r>) 


The  variance  of  XiV^is  given  by 


var  (1&*  =  /'UJyj  -  '  TJm  ? 

=  ^(Ifu.v)  -  tQ2  H (u, v ) )2>  -  W2 (u , v  ) 
W2(u,v)  4  F2(u,v)  4  2tQ?  H(u,v) 

.  [  F ( U , v  ^  4  W(u, v ) ]  . 

(5.2-27) 

For  large  apertures,  var  (\J>)  is  given  to  a  good  approximation 
by  (except  at  the  origin  1 


var  (UVfi)  ~  W2(u,v) 


(5.2-28) 


a  result  which  follows  from  the  argument  leading  to  Eq.  (5.2-19) 
and  from  the  fact  that  H(u,v)  is  very  small  away  from  the 
origin  for  large  apertures.  The  variability  of  ■UV)  is  quite 
pronounced  since 


var  ( U*  )  m  W2  (u.v)  ___  1 

.  <vr>2  w2  (u,v) 
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6.  EXPERIMENTAL  TECHNIQUES  AND  APPARATUS 


The  feasibility  of  reliable  spectral  density  measure¬ 
ments  by  diffraction  is  borne  out  by  the  theoretical  considera¬ 
tions  of  the  earlier  chapters.  The  results  of  chap.  4  provide 
a  basis  against  which  experimental  results  can  be  compared. 

Such  experimental  results  are  given  in  chap.  7.  The  present 
chapter  is  a  discussion  of  the  laboratory  apparatus  and  pro¬ 
cedures  used  to  obtain  the  experimental  data. 

The  first  section  gives  a  brief  description  of  the  ex¬ 
perimental  apparatus  used  for  obtaining  spectral  density  meas¬ 
urements.  Section  6.2  describes  the  basic  measuring  procedure  and 
its  advantage  over  more  direct  methods.  Section  6.3  summarizes 
the  techniques  by  which  various  film  noise  parameters  can  be 
measured.  Section  6.4  describes  the  use  of  contoured  aper¬ 
tures  and  aperture  functions  to  achieve  "desirable"  zero-order 
diffraction  patterns  since  intense,  spatially  distributed,  zero- 
order  diffraction  levels  preclude  the  measurement  of  power  spec¬ 
tra.  Finally,  Sec.  6.5  contains  a  description  of  a  magnification 
technique  by  which  low-noise  spectra  can  be  enhanced. 

6.1  DESCRIPTION  OF  MEASUREMENT  EQUIPMENT 

The  spectral  density  measurements  were  performed  on  the 
coherent  optical  system  shown  schematically  in  Fig.  6.1-1  . 

The  same  basic  system  has  been  extensively  described  in  the 
thesis  of  J.  Minkoff  (Minkoff  1967 ) .  A  brief  description  of 
the  system,  together  with  some  peripheral  equipment  is  provided 
here . 

The  basic  system  consisted  of  a  granite  bench  upon  which 
was  mounted  a  He-Ne  (6328  A)  laser,  spatial  filters,  collimating 
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and  integrating  lenses ,  a  scanning  assembly  and  a  photomulti¬ 
plier  tube.  The  output  of  the  photomultiplier  tube  fed  an 
(x,y)  recorder  through  an  amplifier  which  could  be  operated  in 
a  linear  or  logarithmic  mode.  All  the  optical  components  and 
light  beams  were  shielded  against  extraneous  light,  dust,  and 
air  currents  by  cylindrical-bellow  light  shields  and  light¬ 
tight  boxes.  As  an  added  precaution,  all  unnecessary  lights 
were  turned  off  during  the  actual  experiments  . 

The  optical  bench  was  16  ft  x  18  in.  x  18  in.  with  a 
steel  base  mounted  in  tripod  fashion  on  air-operated,  self- 
regulated  damping  mounts.  The  top  surface  of  the  bench  was 
ground  to  a  flatness  of  ±  0.0001  in.  A  T-shaped  slot  was 
machined  into  the  bench  to  accommodate  the  optical  components . 

One  surface  of  this  slot  provided  a  reference  surface  for  the 
components  and  was  machined  to  a  flatness  of  .001  in.  The  ac¬ 
tual  system  showing  the  laser,  bellows,  etc.,  is  shown  in 
Fig.  6.1-2  . 

The  laser  provided  a  1.8-mm  diameter  vertically  polar¬ 
ized  beam.  This  beam  was  focused  by  a  spreading  lens  (micro¬ 
scope  objective)  onto  a  20u.  pinhole  serving  as  a  spatial  low- 
pass  filter  (Fig.  6.1-3  )•  Beyond  the  pinhole,  the  beam  spread 
sufficiently  to  cover  the  entire  collimating  lens.  The  loca¬ 
tion  of  the  pinhole  in  transverse  directions  was  controlled  by 
a  set  of  positioning  screws  and  the  axial  location  was  con¬ 
trolled  by  positioning  a  barrel-type  pinhole  mount  along  the 
optic  axis  (Fig.  6.1-4  ).  The  pinhole  was  positioned  at  the 
focal  point  of  the  objective  in  a  way  that  would  allow  uniform 
illumination  to  appear  across  the  entire  collimating  lens  aperture. 

The  data  sample  was  located  in  the  region  between  the 
collimating  and  integrating  lenses  -  the  so-called  collimated 
region.  The  sample  was  rigidly  held  between  two  optically  flat 
glass  pieces  and  mounted  on  a  rotary  table  which  provided  six 
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1 .  Lens  Mount 

2.  Microscope-Barrel  Assembly 

3.  Coarse  Focusing  Adjustment 

4.  Fine  Focusing  Adjustment 

5.  Pinhole  Assembly 

6.  Pinhole  Positioning  Screws 

7.  Microscope  Objective  Lens 

8.  Objective-Lens  Holder 

9.  Spreader-Lens-Assembly  Mounting  Plate 

10.  Lens-Mount  Locking  Knob 

11.  Mounting-Plate  Adjustment  Screws 


I 

\ 

J 


Fig.  6.1-4  Spreader-Lens  Assembly 
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degrees  of  freedom  for  precise  positioning.  The  aperture  stop 
was  located  immediately  after  the  collimating  lens  and  was  in¬ 
terchanged  according  to  the  type  of  apodisation  used  (Sec. 

6.3^.  The  phase-eliminating  liquid  gate  [see  Fiq.  4.1-1]  was 
located  in  the  collimated  region  and  is  shown  in  Fig.  6.1- 3 
with  the  film  sample  in  place.  The  light  transmitted  by  the 
sample  was  then  focused  on  the  x'yf  plane  by  the  integrating 
lens  and  detected  with  the  scanning  mechanism  shown  in  Figs. 
6.1-7  and  6.1-7  .  Scanning  was  possible  in  two  orthogonal 
directions  and  it  was  possible  to  allow  any  size  scanning  slit 
to  be  mounted  on  the  scanning  assembly.  Between  the  scanning 
slit  and  the  photomultiplier  tube  were  located  adjustable  neu¬ 
tral-density  filters  to  limit  the  peak  intensity  of  light  en¬ 
tering  the  photomultiplier.  By  adjusting  the  neutral  density 
filters,  the  light  level  at  the  photomultiplier  was  held  nearly 
constant  despite  the  varying  optical  densities  of  the  various 
samples.  The  entire  scanning  mechanism  was  housed  in  a  light¬ 
tight  box.  System  alignment  was  accomplished  using  an  auto¬ 
collimator  telescope  shown  in  Fig.  6.1-6  . 

The  transmittance  of  the  sample  was  measured  by  an  auto¬ 
matic  recording  microdensitometer  (Joyce  Ltd.,  Model  MK  111C) 
shown  in  Fig.  6.1-8  .  The  microdensitometer  provided  a  vol¬ 
tage  proportional  to  optical  density,  which  is  related  to  inten¬ 
sity  transmittance  by 

Density  =  -log10T 

The  square  root  of  T  is  denoted  as  the  average  amplitude 
transmittance  T^  but  for  reasons  given  in  Sec.  6-3,  the 
amplitude  transmittance  obtained  by  incoherent  illumination 
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1.  Light  Shield 

2.  Slit-Traversing  Knob 

3.  Traversing-Mechanism  Connecting  Shaft 

4.  Pellicle-Assembly  Support  Bar 

5.  Mirror-Position  Lever 

6.  Window  for  Pellicle  Assembly 

7.  Photomultiplier  Voltage  Divider  Unit 


Fig.  6.1-6  Scanning-Mechanism  Assembly 
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FIG  6.1-7  SCHEMATIC  REPRESENTATION  OF  SYSTEM  FOR  MEASUREMENT 
BACK-FOCAL  PLANE  LIGHT  INTENSITY 
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was  not  assigned  the  same  symbol  as  the  transmittance  measured 
with  the  laser.  Microdensitometer  traces  were  used  to  determine 
the  uniformity  of  film  exposure  (Fig.  6-1-9). 

In  some  cases,  grain  structure  of  film  samples  were 
altered  by  magnification.  Magnification  was  done  by  two 
separate  means,  both  involving  incoherent  illumination.  Low- 
resolution  magnification  was  obtained  with  a  Wollensak  "Raptar" 

F  A  .  5 ,  50-rom  focal  length  lens.  High  resolution  was  obtained 
with  an  American  Optical  "Fluorestar"  microscope,  Cat.  No. 
LlOTU-FDW.  The  microscope,  power  supply,  and  exciter  are  shown 
in  Fig.  6.1-10  . 

6.2  MEASUREMENT  TECHNIQUE  AND  CALIBRATION 

In  this  section,  the  basic  measurement  technique  is 
described.  A  rectangular  aperture  of  area  LW  is  assumed. 

If  L=W  the  results  can  be  applied  to  a  diamond  aperture  of 
the  same  area  provided  that  scanning  is  then  along  the  line 
uv.  For  the  analysis  that  follows,  it  is  sufficient  to  deal 
only  with  the  relative  light  intensity  l(u,v),  so  that  the 
factor  J2/(^f)?  of  Eq .  (2.4-13)  will  not  appear  in  any  of 
the  equations . 

The  zero-order  light  diffraction  pattern  of  a  rectan¬ 
gular  aperture  of  length  L  and  width  W  is  given  by 

I^(u,v)  =  H(u,v)  (6.2-1) 


where  H(u,v)  =  LW  sinc2uL • sinc2vW .  When  the  same  system  is 
loaded  with  an  object  of  average  complex  transmittance  Rq  the 
resulting  average  diffraction  intensity  is 


<l(u,v)> 

W(u,v) 


I  Rq  I  2  H  (u,  v )  +  W  (u ,  v  ) 

W(u,v)*H(u,v) 
apparent  spectrum. 


(6.2-2) 
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KODAK  2479  FILM  ; 


AVERAGE  INTENSITY  TRANSMITTANCE,  ^=0.28 
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FIG.  6.1-9  INTENSITY  TRANSMITTANCE  MEASUREMENTS  WITH  A 
MICRODENSITOMETER  .  (A)  10  MICRON  SLIT  >  (B)  50 

MICRON  SLIT. 


Power  Supply 
Camera 

Camera  Adapter 
Mercury-Xenon  Light  Source 
Trinocular  Microscope 
Film  Sample 


Microscope  and  Camera  Assembly 


If  the  system  is  unloaded  and  a  point  scanning  slit  is  tempor¬ 
arily  assumed,  the  output  of  the  log  amplifier,  when  the  slit 
is  at  (u,v),  will  produce  a  displacement  on  the  recorder 
given  by 


C1(u,v)  =  10  log  H(u,v) 

where  the  logarithm  is  to  base  10.  When  the  sample  is  inserted 
into  the  system,  the  displacement  will  be 

C2(u,v)  =  10  log  [h(u,v)  |Rq|2  +  W (u, v )] 

If  the  recorder  is  now  offset  by  a  displacement  Adb  so  that 
C2(0,0)  +  Adb  =  C1(0,0)  then  it  follows  that  (Adb  =  10  log  a) 


10  log  [A  |Rq|  2  LW  +  AW(0,0)  ]  =  10  log  LW 

=  C2' (0,0) 


(6.2-3) 


If  LW  |Rj  yy  w(0,0),  a  condition  which  is  easily  met  in 
most  cases,  then  to  a  high  degree  of  accuracy. 


A 


(6.2-4) 


To  get  an  idea  of  the  relative  magnitude  of  these  quantities, 
assume  that  W(u,v)  ~  W(u,v)  and  that  the  random  checkerboard 
representation  of  film  grain  noise  is  used  in  an  estimate  of 
W(0,0).  Then  from  Sec.  4.2 


W(0,0)  =  Is  tQ(l  -  tQ) 

where  i 2  =  the  area  of  a  grain  and  tQ  is  the  average  trans¬ 
mittance  of  the  sample.  Assume  the  following  typical  values 
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l  =10  microns 
W  =  L  =  50  mm 

Then  w(0.0)  is  74  db  below  the  peak  of  the  zero-order  light 
and  the  approximation  in  Eq .  (6.2-4)  is  justified. 

The  two  traces,  C  (u,v)  and  C2'  (u,v),  thus  obtained 
have  the  following  properties.  The  addition  of  Adt)  effec¬ 
tively  compensates  for  the  average  density  of  the  film  without 
affecting  the  power  spectrum  which  results  from  transmittance 
fluctuations  about  the  average.  If  the  sample  were  noiseless, 
the  traces  would  coincide  at  all  spatial  frequencies.  Since 
the  sample  is  not  noiseless,  the  curve  C2'  (u,v)  will  be  above 
C1(u,v).  Suppose  that  the  scanning  is  done  along  the  u-axis 

(v  0)  and  there  exists  a  frequency  u  such  that  for  u  >  u  , 

o  o* 

C„' (uQ,0)  Ci(uo,0)  +  10  db;  then  the  error  in  interpreting 
C '  (u,0)  as  10  log  w(u,0)  +  10  log  [  Rq  |-^  is  11$  or  less. 

This  technique  has  the  very  strong  advantage  of  avoiding  any 
kind  of  operation  on  the  data  except  compensating  for  |Rq|  . 

This  means  that  the  presence  of  unwanted  noise  such  as  might 
be  contributed  by  random  lens  errors,  scattering  from  mirrors, 
etc.,  which  add  significantly  to  the  zero-order  light  away  from 
the  origin,  (Fig.  6.2-1  ),  does  not  upset  the  validity  of  the 
measurement, since  the  power  spectral  light  component  must  be 
a  factor  of  9  or  more  above  the  diffracted  light  of  the  un¬ 
loaded  system  before  it  is  interpreted  as  such.  Fig.  6.2-2 
summarizes  the  measurement  procedure.  Smoothed  light  patterns 
are  shown  as  they  would  appear  through  a  slit  of  nonzero  area. 

A  disadvantage  of  this  procedure  is  that  the  spectrum 
W(u,v)  may  not  become  observable  in  the  sense  of  the  10-db 
criteria  of  the  previous  paragraph  until  the  zero-order  light 
is  relatively  small,  which  may  occur  at  fairly  high  frequencies. 
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a)  TYPICAL  RECORDER  TRACES  OF  LOADED  AND  UNLOADED  OPTICAL  SYSTEM 

BEFORE  COMPENSATION 


b)  SAME  AS  (a)  EXCEPT  FOR  COMPENSATION  FOR  |R0 1 "2 


c) 


RESULTING  SPECTRUM 
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FIG.  6. 2-2  ILLUSTRATING  POWER  SPECTRUM  MEASUREMENT  TECHNIQUE 


For  example,  an  estimate  may  be  obtained  by  taking  the  random 
checkerboard  configuration  as  a  model  of  film  grain  noise  and 
assuming  diffraction  by  a  rectangular  aperture  of  area  LW  for 
which  the  DC  value  of  the  power  spectrum  will  be  10  db  above  the 
zero-order  intensity  at  spatial  frequencies  satisfying 


t  (1  - 

o ' 


\ 


9LW  t  2 
o 

(7IUL)  P 


which  for  the  numbers  previously  chosen  gives 


u  y  100  cycles/mm  for  l  =  10ul 
u  N  10  cycles/mm  for  l  =  lOOu, 


The  discussion  so  far  has  assumed  a  point  scanning  slit 
and  that  W(u,v)  and  not  its  estimate  was  available  in  the 
(u,v)  plane.  For  calibration  purposes,  the  effect  of  slit  area 
must  be  taken  into  account.  The  slit  is  chosen  large  enough 
to  reduce  the  variability  of  the  estimate  so  that  repeated  ex¬ 
perimental  runs  yield  the  same  result,  and  not  so  large  as  to 
smooth  variations  in  the  power  spectrum.  In  measuring  the  grain 
noise  of  film,  the  slit  may  be  several  times  the  width  of  the 
main  lobe  of  the  zero-order  diffraction  pattern.  If,  simultan¬ 
eously,  the  variation  in  the  power  spectrum  over  the  frequency 
span  of  the  slit  is  sufficiently  small  to  enable  the  spectrum 
to  be  approximated  by  its  value  at  the  center  frequency  of  the 
slit  filter,  then  the  calibration  of  the  spectrum  is  straight¬ 
forward  as  will  now  be  shown. 

From  Eq.  (3*6-5)*  the  light  observed  through  a  slit  of 
area  Ag,  described  by  the  filter  function  '-r(ul  ,v'  )  is 

I  (u,v)  =  j0'  l(u'  ,v'  )  g(u-u'  ,v-v'  )  du' dv' 

°  S  -or  -or 

(6.2-5) 
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The  observed  zero-order  light  at  the  origin  is  then 


I 


o 


(0) 


(0,0) 


fa  f°°  H(u'  ,v'  )  g(-u'  ,-v'  )  du'  dv* 

00  —00 

(6.2-6) 


With  the  assumptions  of  the  preceding  paragraph,  there  results 


(0) 


(0,0)  = 


(\f): 


(6.2-7) 


The  observed  apparent  spectrum  at  the  origin  is  given  by 

(u'  ,v'  )  g(-u'  , v '  )  du' dv' 

f+  f+  g(u'  ,V  )  du'  dv' 

—  Or  —Or 

(6.2-8) 

If  the  measurement  is  properly  designed,  then  w(0,0)  ~  w(0,0). 
Following  the  measurement  procedure  outlined  in  this  section, 
assume  that  the  extrapolated  value  of  the  curve  C2' (u,0), 

(Fig.  6.2-2)  is  -r  db  below  the  peak  zero-order  light. 

Then  it  follows  that  the  suitable  calibration  in  the  case  of 
grain  noise  is 


wo(0,°)  -  /+7+” 

S  —or  —  or 

~  i-llll  W(0,0^ 
As 

~  W(0,0) 


or 


10  log 


t  W<0,0) 
o  ' 

)2/As 


-r 


(6.2-9) 


W(0,0) 


10 


-r/io 


(6.2-10) 
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Eq.  (6.2-10)  is  then  the  proper  calibration  for  the 
spectrum.  Observe  that  the  spectrum  has  units  of  meter 
From  Eq .  (6.2-9) >  it  may  be  seen  that  the  ratio  of  power  spec¬ 
tral  to  zero-order  light  at  the  origin  depends  on  the  param¬ 
eters  X,  f.  and  A  ;  since  A  is  usually  determined  by 

S  S 

other  considerations  (Sec.  3*6)  and  A  is  the  wavelength  of 
the  coherent  source  and  is  therefore  fixed,  improvement  in 
this  ratio  can  best  be  achieved  by  using  lenses  with  short 
focal  lengths . 

From  the  discussion  so  far,  it  can  be  seen  that  a  central 
problem  in  spectral  measurements  by  diffraction  is  the  observation 
of  the  spectral  component  of  light  in  the  presence  of  the  strong 
zero-order  diffraction  of  the  aperture.  Measurements  of  the 
spectrum  at  the  origin,  where  the  zero-order  is  a  maximum,  cannot 
be  made  unless  phase-grating  type  aperture  functions  are  used 
(Sec.  6.5).  In  general,  with  the  origin  excluded,  it  is  desir¬ 
able  to  use  a  large  aperture  which,  in  addition  to  adding  to  the 
resolution  capability  of  the  system,  would  concentrate  the  zero- 
order  light  in  the  vicinity  of  the  origin  while  reducing  the  en¬ 
velope  of  the  sidelobes .  The  envelope  of  the  relative  diffraction 
pattern  of  a  square  aperture  is  the  curve  connecting  the  local 
maxima  of  the  function  H(u,v)  and  is  given  by 


Z(u,v) 


u,v  4  0 


which  decreases  as  (lw)  •  Apodisation  by  aperture  shaping 
can  be  used  to  reduce  the  envelope  even  further.  This  is  dis¬ 
cussed  in  Sec.  6.4. 

6 .3  DETERMINATION  OF  FILM  NOISE  PARAMETERS 

The  purpose  of  this  section  is  to  show  how  various  film- 
noise  parameters  can  be  determined  experimentally. 
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(a )  Measurement  of  Intensity  Transmittance 
of  Film  ~ 

The  intensity  transmittance  of  a  film  sample,  denoted 
by  T^Xjy)  is  the  ratio  of  the  transmitted  to  incident  light 
intensity  at  the  point  (x,y).  Sample  traces  of  |T(x,y)| 
can  be  obtained  with  the  microdensitometer  described  in 
Sec.  6.1.  When  the  microdensitometer  is  used  with  a  large 
scanning  aperture,  the  traces  will  have  few  fluctuations  and 
essentially  furnish  the  average  value  of  |T(x,y)|  -  T  .  When 
a  small  scanning  aperture  is  used,  the  resulting  trace  furnishes 
qualitative  information  regarding  the  grain  noise  of  the  film. 


It  should  be  noted  that  the  square  root  of 


is  not 


the  average  amplitude  transmittance  t  ,  defined  as  the  aver¬ 
age  magnitude  of  the  ratio  of  transmitted  to  incident  amplitude 
despite  the  tendency  for  some  authors  to  ignore  the  distinction 
(Kozma  I966).  The  inequality  can  be  demonstrated  as  follows.  Since 


T  (x,  y )  - 


iT(x,y) 


ET(x,y) 


-  t  (x,y) 


it  follows  that 


Tq  =  <T(x,y)>  =  tQp  +  a/ 


where  a  r  =  S (t(x,y)  -  t  )?  )> .  Hence  t  /  T  £.  If  t  and 
t  r  can  be  measured  with  sufficient  accuracy,  the  above  relation 
provides  a  scheme  by  which  o  p  could  be  determined.  Exper¬ 
iments  have  shown,  however,  that  the  measured  values  of  TQ 
obtained  with  a  microdensitomefer  are  not  reliable  enough  to 
provide  an  accurate  value  of  o  p . 

Sample  functions  of  intensity  transmittance  can  be  used 
to  test  the  uniformity  of  film  transmittance.  Uneven  exposures 
result  in  systematic  transmittance  variations  £i.e.,  Tq  =  TQ(x,y)J 
which  can  easily  be  observed  in  the  sample  function. 


-121- 


(b )  Measurement  of  the  Magnitude  of  the  Average 


Complex  Transmittance,  |r 


The  complex  transmittance  is  defined  by  Eq.  (2.1-5); 
the  magnitude  of  the  average  complex  transmittance  can  be 
measured  with  the  coherent  optical  system  shown  in  Fig.  6.1-2 
by  measuring  the  attenuation  in  the  peak  zero-order  light  when 
the  sample  is  inserted  in  the  collimated  region.  In  this  meas¬ 
urement  no  liquid  gate  is  used.  It  is  seen  from  Eqs.  (6.2-1) 
and  (6.2-2)  that  if  W(0,0)  is  neglected,  the  ratio  of 
<1(0,0))  to  I(0)(0,0)  is  lRo|a. 


( c )  Measurement  of  the  Average  (Coherent) 

Amplitude  Transmittance,  t 

The  amplitude  transmittance  t  is  measured  by  the 
same  technique  used  to  measure  |  Rq  |  except  that  the  film 
sample  is  immersed  in  a  matched  liquid  gate  which  has  the  ef¬ 
fect  of  negating  phase  noise.  In  this  case  Rq  =  tQ  and  the 
attenuation  in  the  peak  zero-order  light  when  the  film  sample 
is  in  the  system  is  tQ2 . 

(d)  Determination  of  the  Variance  of  the  Random 

Phase  Process  {#(x,y)f  ' 

When  the  processes  \  |t(x,y)f  and  |V'(x,y)}  are  con¬ 
sidered  independent  and  j is  a  wide-sense,  zero-mean 
Gaussian  process,  then 


<R(x,y)>  =  Rq 


tQ  e"(1/2)a/  [Sec.  4.3,  Eq.  (4.3-3)]. 


where 


a  ^  =  the  variance  of  |^(x,y)} 


Under  these  conditions,  a  graph  of  R  vs  t  should  be  a 

-(1/2)S/ 


straight  line  with  slope  m  =  e 


Hence 


obtained  from 


V' 


can  be 


if 


=  2  In  1/m 
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(e )  Determination  of  the  Isotropic  Statistics 
of  Film  Noise 

If  film  noise  is  isotropic,  the  power  spectrum  is 
a  function  of  the  frequency  w  where  w  =  (u2  +  v2)^.  In 
this  case,  the  power  spectrum  is  an  even  function  of  spatial 
frequencies  and  depends  only  on  the  distance  from  the  origin 
of  the  spatial  frequency  plane.  In  general,  when  the  noise 
is  anisotropic,  the  two-dimensional  power  spectrum  will  not  be 
an  even  function  of  frequency  as  discussed  in  Sec.  2.3  and  will 
not  be  a  symmetric  function  of  its  arguments.  Hence,  anisotropic 
noise  can  be  detected  experimentally  by  observing  asymmetric  spectra 
about  the  origin  or  unequal  spectra  along  different  orthogonal  axes. 

( f )  Variations  of  the  Power  Spectrum  with 
Transmittance 

Both  the  random  checkerboard  and  overlapping  circular 
grain  model  predict  a  maximum  in  the  power  spectrum  of  grain 
noise  in  the  vicinity  of  tQ  -  0.3.  Following  the  procedure 
outlined  in  Sec.  6.2,  this  can  be  verified  experimentally  by 
obtaining  the  spectra  at  fixed  spatial  frequencies  for  a  large 
number  of  film  samples  with  varying  transparencies.  The  varia¬ 
tion  of  the  spectrum  of  the  total  noise  with  Rq  can  be  ob¬ 
tained  with  the  same  procedure,  provided  that  the  liquid  gate 
is  not  used. 

(g )  Measurements  by  Magnification 

Under  certain  conditions,  magnification  can  be  used  to 
obtain  the  spectra  of  low  noise  films  and  the  frequency  response 
of  optical  systems.  This  technique  is  discussed  further  in 
Sec.  6.5. 

6.4  REDUCTION  OF  ZERO-ORDER  SIDELOBES  BY  APERTURE  SHAPING 

The  reduction  of  the  zero  order  light,  a  problem  of 
central  importance  in  spectral  measurements  by  diffraction. 
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can  be  accomplished,  in  principle  at  least, in  two  ways.  The 
most  direct  way  is  to  decrease  the  dimensions  of  the  aperture 
diffraction  pattern  by  increasing  the  transverse  dimensions 
of  the  aperture.  However,  in  practice,  the  aperture  is  usually 
already  as  large  as  can  be,  being  limited  by  the  size  of  the 
lens,  the  uniformity  of  the  collimated  beam,  and  the  off-axis 
capability  of  the  optical  system. 

Much  study  and  literature  have  been  devoted  to  the  other 
method,  which  suppresses  the  zero-order  light  by  "apodisation" 
(Jacquinot  and  Roizen-Doz zier  I966),  by  which  is  meant  the  ap¬ 
plication  of  pupil  or  aperture  functions  with  sidelobe-suppress- 
ing  properties.  Sidelobe  reduction  is  usually  achieved  at  the 
expense  of  main-lobe  broadening,  which  has  the  effect  of  di¬ 
minishing  the  resolving  power  of  the  system. 

The  pupil  function  can  be  generated  by  an  aperture  with 
uniform  (unity)  transmission  together  with  a  suitable  contour 
or  by  an  absorbing  filter  of  appropriate  transmission  with  a 
"normal"  contour.  The  latter  technique  offers  the  advantage 
of  not  confining  the  desired  effect  to  any  one  or,  at  most,  a 
few  particular  directions  in  the  frequency  plane  as  is  charac¬ 
teristic  of  apodising  by  contour  shaping.  However,  the  construc¬ 
tion  of  accurate  abosrbing  filters  is  in  itself  a  difficult  task 
since  great  care  must  be  taken  to  avoid  filter-generated  noise. 
The  other  possibility  -  apodisation  by  aperture  shaping  -  avoids 
this  problem,  providing  the  contour  is  reasonably  constructed, 
as  is  not  difficult  in  practice. 

When  contoured  apertures  are  used,  the  power  spectrum 
should  be  measured  along  a  direction  of  apodisation.  If  the 
random  process  has  isotropic  statistics,  the  measurement  of  the 
spectrum  in  any  direction  is  sufficient  to  obtain  second  moment 
information  about  the  process;  if  the  process  is  not  isotropic, 
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then  rotating  the  sample  with  respect  to  the  aperture  will 
effectively  allow  the  measurement  of  the  power  spectrum  along 
different  directions. 

It  will  now  be  shown  how  apodisation  by  aperture  shap¬ 
ing  is  possible.  The  diffraction  amplitude  of  an  aperture  of 
area  A  with  uniform,  unity  transmission,  is  given  by 

Ef(u,v)  -  K  f  f  e-32;‘(ux+vy)dxdy  (6.4-1) 

A 

where  the  constant  K  required  for  energy  conservation  is 
of  no  interest  in  this  discussion  and  is  therefore  set  equal 
to  one . 

When  the  diffraction  intensity  of  a  square  aperture  of 
side  D  is  scanned  by  a  point  aperture  along  v  ^  0,  the 
intensity  is 

If(u.v)  d"  (siB. j£SbV  (6.4-2) 

_  2 

The  envelope  is  proportional  to  ( -ud)  and  decreases  at  the 
rate  of  -  6  db/octave.  A  diamond  aperture  with  the  same  area 
produces  a  diffraction  intensity  with  an  envelope  function  pro¬ 
portional  to  (ttud)-1  which  decreases  at  12  db/octave.  The 
diffraction  pattern  of  a  diamond  aperture  can  be  obtained  by 
scanning  the  diffraction  pattern  of  a  square  aperture  along  the 
line  u  v.  The  diamond  aperture  wa«  the  first  and  simplest 
example  of  apodisation  along  a  special  direction  fcouder  and 
Jacquinot  1939).  An  adjustable  diamond  aperture  constructed 
for  the  experiment  is  shown  in  Fig.  6.1-5  . 

The  diffraction  pattern  resulting  from  an  aperture 
with  contour 

y  +  1/2  C(x)  -a/2  '  x  <  a/2 

0  otherwisa 


(6.4-3) 
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is  given  by 


E.(u,v)  =  Ja/2  dx  /c(x)/2 

f  -a/2  -C(x(/2 


e~  j2?rvy  fly 


which  gives,  along  the  u-axis 


E.fUjO)  =  /a/2  C (x )  e"^2  ux  dx 
f  -a/2 


(6.4-4) 


This  is  the  one-dimensional  Fourier  transform  of  C(x)  over 
the  interval  f-a/2,a/21.  Complete  sidelobe  elimination  c  n 
be  achieved  with  the  Gaussian  aperture  defined  by 

y  -  D  e_x2//2a"  -a-  <  x  <  a  (6.4-5 ) 


which  gives  a  diffraction  amplitude  of  the  form  e-2 ^  u  ' 

The  diffraction  patterns  of  square,  diamond,  and  Gaussian  aper- 
turesof  equal  area  are  shown  in  Fig.  6.4-1  .  In  practice,  the 
aperture  must  always  be  of  finite  extent  and  the  diffraction 
pattern  is  then: 


Efr(u,0)  =  /+Pr  e  x  /2a  rect  x/a  e-12"Tux  dx 

£  —  or 

=  S (u)  *  a  sine  ua 


(6.4-6) 


where 

S  (u )  =  >/5T  crD  e*2(,Tuo)* 

Gaussian  apertures  used  in  the  experiment  are  shown  in  Fig. 
6.4-2  .  The  widths  are  approximately  3cr. 

In  practice,  aperture  shaping  does  not  provide  a  com¬ 
plete  solution  to  the  zero-order  light  suppression  problem. 
First  of  all,  the  source  light  at  the  origin  remains  unaffected 
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FIG. 6. 4-1  ENVELOPE  OF  DIFFRACTION  PATTERNS  FOR  THREE  APERTURES 

(THEORETICAL) 
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and  may  be  many  orders  of  magnitude  above  the  spectral  light 
component,  thus  prohibiting  its  measurement  at  or  near  zero  fre¬ 
quency.  ■*■  Also,  apodisation  cannot  compensate  for  light  scattered 
by  random  imperfections  in  the  lenses.  Regardless  of  the  care 
and  precision  with  which  contoured  apertures  are  constructed, 
the  light  levels  in  the  diffraction  pattern  cannot  be  reduced 
further  than  the  basic  scattering  level  of  the  optical  system. 

For  the  system  described  in  Sec.  6.1,  this  level  was  found  to 
be  of  the  order  of  70-80  db  below  the  peak  incident  light  for 
the  largest  off-axiu  scan  displacements  of  interest. 

6.5  USE  OF  MAGNIFICATION  IN  POWER  SPECTRA  MEASUREMENTS 

There  are  certain  situations  where  it  is  desirable  to 
magnify  the  film  grains  before  attempting  to  measure  the  power 
spectrum  by  diffraction.  For  example,  it  may  be  difficult,  if 
not  impossible  to  measure  the  spectrum  of  low-noise  film  in 
which  the  grains  are  too  small  to  scatter  sufficient  light  that 
can  be  distinguished  from  the  zero  order.  Even  if  the  film  is 
such  that  the  spectrum  of  its  noise  can  be  measured,  the  off- 
axis  capability  of  the  optical  system  may  preclude  the  observa¬ 
tion  of  the  spectrum  at  high  spatial  frequencies.  From  Eq . 

(2.4-2)  it  can  be  seen  that  high  spatial  frequencies  corres¬ 
pond  to  large  transveise  displacements  in  the  frequency  plane 
which  result  from  large  diffraction  angles.  The  response  of 
most  lenses  becomes  poorer  for  large  diffraction  angles  and 
hence  this  mode  of  operation  should  be  avoided. 

It  is  shown  in  this  section  that  if  the  film  grains  of 
a  sample  are  magnified  and  transferred  to  low-noise  film,  the 
spectrum  of  the  magnified  grains  can  often  be  easily  obtained 
by  diffraction  and  after  suitable  calibration  can  be  interpreted 
as  the  spectrum  which  would  have  been  produced  by  unmagnified 
grains.  Magnification  also  effectively  "compresses"  the  fre¬ 
quency  axis  so  that  much  higher  spatial  frequency  components  or 
the  original  spectrum  can  be  observed. 


+  See  Appendix  B 
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The  above  advantages  accrue  provided  that  the  magni¬ 
fication  is  done  without  distortion.  In  practice,  the  elimin¬ 
ation  of  all  distortion  is  not  possible.  However,  the  use  of 
precision  lenses  with  large  numerical  apertures  can  minimize 
the  distortion  as  will  now  be  shown. 


Consider  first  "ideal"  magnification.  This  involves 
changing  all  dimensions  by  the  magnification  factor  K  so  that 
a  grain  of  radius  rQ  now  has  radius  KrQ .  The  amplitude  trans 
mittance  t(x,y)  of  the  original  process  is  assumed  to  have 
mean  value  t  which  does  not  change  through  magnification 
provided  that  the  magnified  sample  contains  a  very  large  num¬ 
ber  of  grains .  That  the  mean  value  remains  the  same  follows 
rrom  the  observation  that  the  ratio  of  opaque  to  total  area  is 
the  same  before  and  after  magnification.  If  the  covariance  of 
the  original  process  |r(x,y)}  -  |t(x,y)  -  tQ|  is  P  (a,f  ), 
then  the  covariance  of  the  new  process  ■Jt'  (x,y)|  (Fig.  6.5-1 
is 

Tt  (n,P  )  =  4T(a/K,p/K)  (6.5-1) 


and  the  power  spectrum  of 


4*00  --f-Or 


wT,  (u,v) »  r  i 


-ft'  —cr 


=  K-7  *  /  “ 


—  Of  —or 


s 

’  T1 


'  T 


(x,y) } 

(a,p)  e 

(a,p)  e 


is  given  by 

-jS7T(”+vp)  da, dp 


-  j2?r(ua+vp,  )k 

da,  dp 


*=  K"  W  ( Ku ,  Kv  )  . 


(6.5-2) 
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and  ma'’  be  many  orders  of  magnitude  above  the  spectral  light 
component,  thus  prohibiting  its  measurement  at  or  near  zero  fre¬ 
quency.*  Also,  apodisation  cannot  compensate  for  light  scattered 
by  random  imperfections  in  the  lenses.  Regardless  of  the  care 
and  precision  with  which  contoured  apertures  are  constructed, 
the  light  levels  in  the  diffraction  pattern  cannot  be  reduced 
further  than  the  basic  scattering  level  of  the  optical  system. 

For  the  system  described  in  Sec.  6.1,  this  level  was  found  to 
be  of  the  order  of  J0-80  db  below  the  peak  incident  light  for 
the  largest  off-axis  scan  displacements  of  interest. 

6.5  USE  OF  MAGNIFICATION  IN  POWER  SPECTRA  MEASUREMENTS 

There  are  certain  situations  where  it  is  desirable  to 
magnify  the  film  grains  before  attempting  to  measure  the  power 
spectrum  by  diffraction.  For  example,  it  may  be  difficult,  if 
not  impossible  to  measure  the  spectrum  of  low-noise  film  in 
which  the  grains  are  too  small  to  scatter  sufficient  light  that 
can  be  distinguished  from  the  zero  order.  Even  if  the  film  is 
such  that  the  spectrum  of  its  noise  can  be  measured,  the  off- 
axis  capability  of  the  optical  system  may  preclude  the  observa¬ 
tion  of  the  spectrum  at  high  spatial  frequencies.  From  Eq . 

(2.4-2)  it  can  be  seen  that  high  spatial  frequencies  corres¬ 
pond  to  large  transverse  displacements  in  the  frequency  plane 
which  result  from  large  diffraction  angles.  The  response  of 
most  lenses  becomes  poorer  for  large  diffraction  angles  and 
hence  this  mode  of  operation  should  be  avoided. 

It  is  shown  in  this  section  that  if  the  film  grains  of 
a  sample  are  magnified  and  transferred  to  low-noise  film,  the 
spectrum  of  the  magnified  grains  can  often  be  easily  obtained 
by  diffraction  and  after  suitable  calibration  can  be  interpreted 
as  the  spectrum  which  would  have  been  produced  by  unmagnified 
grains.  Magnification  also  effectively  "compresses"  the  fre¬ 
quency  axis  so  that  much  higher  spatial  frequency  components  or 
the  original  spectrum  can  be  observed. 


+  See  Appendix  B 


-129- 


The  above  advantages  accrue  provided  that  the  magni¬ 
fication  is  done  without  distortion.  In  practice,  the  elimin¬ 
ation  of  all  distortion  is  not  possible.  However,  the  use  of 
precision  lenses  with  large  numerical  apertures  can  minimize 
the  distortion  as  will  now  be  shown. 

Consider  first  "ideal"  magnification.  This  involves 
changing  all  dimensions  by  the  magnification  factor  K  so  that 
a  grain  of  radius  rQ  now  has  radius  KrQ.  The  amplitude  trans¬ 
mittance  t(x,y)  of  the  original  process  is  assumed  to  have 
mean  value  t  which  does  not  change  through  magnification 
provided  that  the  magnified  sample  contains  a  very  large  num¬ 
ber  of  grains.  That  the  mean  value  remains  the  same  follows 
from  the  observation  that  the  ratio  of  opaque  to  total  area  is 
the  same  before  and  after  magnification.  If  the  covariance  of 
the  original  process  jr(x,y)|  -  |t(x,y)  -  tQ }  is  ( a , (  ) , 
then  the  covariance  of  the  new  process  |t'  (x,y)t  (Fig.  6.5-1  ) 
is 

(a,p)  -  (a/Ktfi/K)  (6.5-1) 

t  l 

and  the  power  spectrum  of  |t'  (x,y)  }  is  given  by 

WT,  (u,v)  „  r  I+"  fT,  (CP)  e-j27’(ua+VB)  d„,dp 


-  *>/  °  / 


—  or  —or 


'/•T(a,p ) 


-j27i(ua+vp  )K 
fc  da ,  dp 


-  Kp  W(Ku,Kv) . 


(6.5-2) 
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t(x)  =  !0  +  T(x)  <t(x)>  =  t0 


-1 — I  n  i — i  r  r. 


a)  FILM  SAMPLE  WITH  TYPICAL  TRANSMITTANCE  TRACE 


a  t'(x)  =  t0+r'(x)  <  t'(x)>  =  t0 


1 


o 


r~n 


b)  SAMPLE  ANO  TRANSMITTANCE  TRACE  AFTER  MAGNIFICATION  BY  K 

A  COVARIANCE 


A-0Q6-I-00-00C 


FIG. 6.  5-1  THE  EFFECT  OF  LINEAR  MAGNIFICATION 
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Hence,  the  magnification  produces  a  power  spectrum  that  is  K? 
times  as  large  at  the  origin  as  the  original  spectrum  and  K” 
times  as  narrow.  The  total  power  in  the  two  processes  are  the 
same  since 


r  r  wTt  (u>v) dudv  =  r  r  wT(u>v) dudv 

—  oo  —or  —On  —on 

(6.5-3) 


which  follows  from  Eq.  (6.5-1). 

When  magnification  is  done  by  a  real  optical  system, 
the  frequency  response  characteristics  of  the  lenses  must  be 
taken  into  account.  It  has  been  shown  (Hopkins  1955>  Born 
and  Wolf  1964  )  that  the  maximum  frequency  that  an  optical 
imaging  system  can  transmit  is  a  function  of  the  width  of  the 
aperture,  the  wavelength  and  type  (coherent  or  incoherent)  of 
illumination  and  the  location  of  the  image  plane.  For  incoherent 
illumination,  the  frequency  response  of  an  optical  system  is 
the  Fourier  transform  of  the  transmission  function  Bj (xi-xQ ,y 1-yo  ) 
defined  as  the  distribution  of  light  intensity  in  the  image  plane 
at  (x1,y1)  due  to  a  point  source  of  unit  strength  in  the  ob¬ 
ject  plane  at  (xQ,yo).  By  superposition,  the  imaging  can  be 
described  by  the  convolution 

I1(x1>y1)  =  r  r  Bi(xrx0,yi-y0)  dxodyo 

—00  —00 

(6.5-4 ) 

where 

I1(x1>y1)  =  light  intensity  at  (x1,y1) 
in  the  image  plane . 

Io(xo'yo)  =  1;L9ht  intensity  at  (xQ,yQ) 
in  the  object  plane. 
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A  similar  relation  holds  for  coherent  systems,  namely. 


,  (x, ,y, )  =  r  f0'  B  (x  -x  ,y  -y  )  E  (x  ,y  )  dx  dy 
iv  i”i'  J  J  cv  i  o,-ri  Jo  o  o  o  o  c 


—  00  —on 


where 


(6.5-5) 


El(xl,y1)  =  complex  light  amplitude  in  the 
image  plane  at  (x^yj 

E0(xn,y0)  =  complex  light  amplitude  in  the 
object  plane  at  (xQ,yo) 

B  (x^Xg  ,y  1-yQ  )  complex  amplitude  distribution 
c  ''  u  in  the  image  plane  due  to  a  point 

source  of  unit  amplitude  and  zero 
phase  in  the  object  plane. 


In  this  case,  the  frequency  response  is  the  Fourier  transform 
of  Bc(x,y).  Denoting  by  b  (u,v)  and  bc(u,v)  the  fre¬ 
quency  response  characteristics  of  incoherent  and  coherent 
systems,  Hopkins  has  shown  that  for  aberration- free  systems, 
bj(u,v)  is  proportional  to  the  self-correlation  of  the  aper¬ 
ture  function  and  bc(u,v)  is  equal  to  the  aperture  function. 
For  a  circular  aperture  of  diameter  d,  the  aperture  function 
G(x,y)  is  given  by 


G(x,y)  =  1 


=  0 


x'  +  y‘ 


d?/4 


otherwise 


For  coherent  illumination,  all  frequency  pairs 
that  (\  =  wavelength  of  radiation,  R  =  radius 
reference  sphere,  i.e.  ~  the  distance  from  the 
image  plane ) 


p  n  \ 

ur  +  v‘  > 


(u,v)  such 
of  the  Gaussian 
aperture  to  the 


will  not  be  transmitted.  When  the  illumination  is  incoherent 
all  frequency  pairs  such  that 
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will  not  be  transmit  ted, and  the  frequency  response  is  not  uni¬ 
form  as  shown  in  Fig.  6.5-2  . 

In  order  to  derive  the  conditions  under  which  power 
spectra  can  be  magnified  with  high  fidelity  by  imaging,  con¬ 
sider  the  incoherent  optxcai  system  shown  schematically  in 
Fig.  6.5-3  .  The  optical  system  denoted  by  L  can  be  a 
microscope  or  an  ordinary  lens  such  as  those  used  in  the  ex¬ 
periments  of  this  thesis.  Assume  the  coordinates  in  image 
and  object  planes  are  normalized  so  that  the  magnification  is 
unity.  If  the  intensity  of  the  source  is  unity,  it  follows 
from  Eq.  (6.5-4 ),  that  if  T  ( P  ,  *i )  denotes  the  intensity  trans¬ 
mittance  of  the  object  0,  and  ( a ,  P  )  /T(<s,t))  T^+n 

then  the  light  in  plane  P  is  given  by 


Ip(x,y)  =  f"  /"  T(P,,n)  Bj  (x-^  ,y-Tj )  d£d-|  (6.5-6) 

—  r.  —  or 


and 


Wr 


(P) 


(u,v)  |bI(u,v)r  WT(u,v) 


(6.5-7) 


where 


W_(P)(u,v)  = 


the  power  spectrum  of  the  imaged 


process  in  plane  P 
WT(u,v)  ap{0T(a,p.)l 

f|  }  -  Fourier  transform  of  |  |. 


Eq.  (6.5-1)  involves  intensity  spectra.  The  intensity  spectrum 
will  be  equal  to  the  amplitude  spectra  if  the  sample  functions 
of  the  intensity  and  amplitude  transmittances  are  equal  such 

v  A  similar  relation  holds  for  coherent  systems  namely  W  ^P^(u,v) 
|  bc (u,v)l r  Wt(u,v),  where  Wfc(u,v)  is  the  power  spectrum  of  the 

amplitude  transmittance.  bj  is  the  self-correlation  function  of 
t>c  (Born  and  Wolf  1964). 
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bc  (u  to)  *  FREQUENCY  RESPONSE  FOR 
COHERENT  ILLUMINATION 


bc  (u,o) 


bc  (o  ,o) 


b  x<u  ,o)  *  FREQUENCY  RESPONSE  FOR 
INCOHERENT  ILLUMINATION 


A-006- 1 -00-0003 


NORMALIZED  FREQUENCY  RESPONSES  FOR 
SYSTEMS  FREE  OF  GEOMETRICAL  ABERRATIONS 
(HOPKINS,  1952  •  BORN  AND  WOLF,  1964) 


OBJECT  PLANE 


IMAGE  PLANE 


A-  006-1-00-0004 


FIG.  6. 5*3  IMAGING  OF  GRANULARITY  THROUGH  AN 
INCOHERENT  OPTICAL  SYSTEM 
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as  in  the  random  checkerboard  and  overlapping  circular-grain 
models  . 

It  can  be  seen  from  Eq .  (6.6-7),  that  high-fidelity  imag¬ 
ing  of  spectra  requires  that  I  b^(u,v)  I  be  nearly  constant  over 
the  frequencies  that  WT(u,v)  has  significant  amplitude.  If 
Wn(u,vj  is  fiat,  or  nearly  so,  over  all  the  frequencies  for 
which  bjfujv)  is  not  zero,  the  transferred  spectrum  will  be 
proportional  to  I  bI(u,v)r>  |  so  that  the  magnitude  of  the 
frequency  response  can  be  measured  by  measuring  the  magnified 
spectrum . 

For  the  experiments,  two  magnification  schemes  were  em¬ 
ployed.  High  resolution,  wide  band  magnification  was  obtained 
with  the  "Fluorestar"  microscope  shown  in  Fig.  6.1-10  .  An 
indication  of  the  frequency  performance  of  a  high-quality  micro¬ 
scope  is  its  resolving  power  which  is  often  specified  in  terms 
of  the  numerical  aperture  (N.A.);  the  N.A.  can  be  expressed  as 

M  a  objective  diameter  _  _d_ 

"  2x  focus  ~  2f 

The  theoretical  cut-off  frequency  for  an  aberration- free 
objective  is 

w  d/\f  -  2  N.A.  A  w  =*J  u'  -!*  v ” 
c 

Depending  on  the  magnification  used,  the  microscope  objectives 
used  in  these  experiments  had  values  of  w  that  ranged  from 
wc  •  1850  cycles /mm  (N.A.  =  .5)  to  wc  =  4600  cycles/mm 

(N.A.  -  1.25).  These  values  were  adequate  for  observing  the 
power  spectrum  of  film  noise  at  frequencies  not  above  several 
hundred  cycles/mm. 

Low  resolution  magnification  was  obtained  with  an  in¬ 
expensive  commercial  enlarging  lens,  operated  in  white  light 
(see  Sec.  6.1).  When  this  lens  was  used  for  magnification,  the 
resulting  spectra  were  significantly  distorted  by  the  frequency 
response  of  the  lens. 
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7.  EXPERIMENTAL  RESULTS 


This  chapter  contains  the  results  of  the  experimental 
phase  of  the  research.  The  experiments  that  were  performed  can 
be  broadly  classified  into  three  categories.  In  the  first 
(Sec.  7. I-?. 2)  are  those  experiments  designed  to  verify  the 
theoretical  conclusions  regarding  the  influence  of  system 
parameters  on  spectral  density  measurements .  Included  in  this 
group  are  the  determinations  of  the  effects  of  the  liquid  gate, 
aperture  sizes  and  shapes,  and  scanning  slit  widths.  The  sec¬ 
ond  category  (Sec.  ( .  $- 7.4)  contains  experiments  concerned 
with  the  measurement  of  film  noise  and  related  subjects .  Com¬ 
parisons  between  measured  spectra  and  the  spectra  of  several 
mathematical  models  of  granularity  are  furnished.  A  posteriori 
discussions  are  provided  in  the  several  instances  where  the 
results  were  not  anticipated.  The  last  (Sec.  7.5'*  category 
contains  those  experiments  related  to  magnification  and  the 
measurement  of  the  frequency  response  of  optical  systems . 

7.1  REDUCTION  OF  FILM  NOISE  WITH  A  LIQUID  GATE 

The  contention  made  in  earlier  chapters  regarding  the 
reduction  of  film-scattered  light  by  the  use  of  a  liquid  gate 
was  experimentally  verified  by  several  workers  (Thiry  I963, 
Leith  I962).  In  practice,  the  choice  of  liquid  which  furnishes 
the  optimum  refractive  index  can  be  made  by  using  an  interfer¬ 
ometer  in  the  collimated  region  and  observing  the  interference 
fringes  or  by  examining  the  structure  of  the  diffraction  pat¬ 
tern  with  the  film  sample  in  place.  The  second  technique  is 
decidedly  preferable  since  the  diffraction  pattern  is  extremely 
sensitive  to  phase  errors  in  the  coherent  beam. 
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Under  the  condition  of  optimum  match,  the  chief  con¬ 
tributor  to  the  scattered  light  is  the  film  grain  noise  which, 
being  much  smaller  than  the  phase  noise,  would  not  be  expected 
to  distort  the  diffraction  pattern  by  a  significant  degree. 
Hence  little  difference  should  be  discerned  between  the  dif¬ 
fraction  patterns  of  the  loaded  and  unloaded  systems. 

In  Fig  7.1-1  the  diffraction  pattern  of  a  carefully 
aligned,  coherent  optical  system  is  shown.  The  parameters  of 
the  optical  system  in  this  case  were:  f  =  focal  length  of 

lens  :  1  meter,  '  -  6328  A,  d  =  length  of  rectangular  aper¬ 
ture  -  34  mm.  The  null-to-null  width  of  the  main  lobe  is 
~  37  ;!  which  is  2yf/d.  The  distortion  in  the  zero-order 
when  film  is  placed  in  the  collimated  region  (Fig.  2.4-2  )  is 
clearly  seen  in  Fig.  7.1-2  for  which  no  liquid  gate  was  used. 
When  a  liquid  gate  is  used  with  a  mismatched  liquid,  the  dif¬ 
fraction  pattern  still  suffers  serious  deterioration  as  can  be 
seen  from  Fig.  7.1-3  •  Here  the  refractive  index,  n,  of 
the  liquid  was  1  5  while  the  best  average  value  of  n  for  the 
film  was  1.610.  When  this  value  of  n  was  used,  the  diffrac¬ 
tion  pattern  was  very  nearly  the  symmetrical  pattern  of  Fig. 

7,1-1  except  for  the  slight  attenuation  of  the  peak  of  the 
zero-order  component,  a  consequence  of  the  optical  density  of 
the  film.  The  diffraction  pattern  is  shown  in  Fig.  7.1-4  . 

7  .2  EXPERIMENTAL  DETERMINATION  OF  THE  EFFECTS  OF  THE 

APERTURE  AND  SCANNING  FILTER 

It  was  shown  in  Sec.  3-2  that  the  apparent  spectrum  is 
the  convolution  of  the  true  spectrum  and  the  aperture  function 
H(u,v'),  defined  for  a  rectangular  aperture  of  dimensions 
L,  W,  by  Eq .  (3.2-2).  The  apparent  spectrum  is  a  biased  esti¬ 
mate  of  the  true  spectrum  and  the  bias  [Eq .  (3.4-11)]  can  be 
seen  to  depend  on  the  dimensions  L  and  W  of  the  aperture. 
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FIG.  7.1-1  FOCAL-PLANE  LIGHT  INTENSITY  DISTRIBUTION 
ALONG  X*  AXIS  FOR  CLEAR  (UNLOADED)  SYSTEM 
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FIG.  7.1-2 


FOCAL-PLANE  LIGHT  INTENSITY  DISTRIBUTION  ALONG  x' 
AXIS  WHEN  SYSTEM  IS  LOADED  WITH  FILM 
WITHOUT  PHASE  COMPENSATION 
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RELATIVE  LIGHT  INTENSITY  (4h) 


B -006 -I  -00-0050 


FIG.  7.1-3  FOCAL- PLANE  LIGHT  INTENSITY  DISTRIBUTION 
ALONG  X*  AXIS  WHEN  SYSTEM  IS  LOADED  WITH 
FILM  WITH  PARTIAL  PHASE  COMPENSATION 
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RELATIVE  LIGHT  INTENSI TY  (dfe) 


FIG.  7.1-4  FOCAL- PLANE  LIGHT  INTENSITY  DISTRIBUTION 
ALONG  X’  AXIS  WHEN  SYSTEM  IS  LOADED  WITH  FILM 
WITH  OPTIMUM  PHASE  COMPENSATION 
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The  bias  error  varies  inversely  with  the  size  of  the  aperture, 
becoming  zero  when  the  aperture  is  infinitely  large  [Eq.  (3*2-7^)* 
Figure  7*2-1  shows  the  effect  of  different  aperture  sizes  and 
shapes  on  the  measured  spectrum.,  It  can  be  seen  that  the  four 
apertures  listed  in  Table  7.2-1  furnished,  despite  minor  differ¬ 
ences,  essentially  similar  results  at  most  frequencies  despite 
radical  differences  in  size  and  shape.  The  results  demonstrate 
the  earlier  contention  that  when  the  aperture  dimensions  are 
large  compared  with  the  largest  significant  correlation  intervals, 
the  measured  spectral  estimates  are  relatively  independent  of 
the  aperture  and  reflect  a  fundamental  property  of  the  film. 


TABLE  7.2-1 


SHAPE 

WIDTH  OF  SIDE  OR  BASE 
(mm ) 

AREA 

(mm° 

Gaussian 

30 

400 

Diamond 

30 

900 

Gaussian 

50 

1100 

Diamond 

50 

2500 

It  will  be  remembered  from  Sec.  3*2  that  when  the  aver¬ 
aging  over  frequency  is  confined  to  a  single  point  as  in  Eq. 
(3*2-7)>  the  bias  is  zero.  The  best  approximation  to  Eq.  (3*2-7) 

is  achieved  with  the  50-nun  diamond  aperture  which  furnished  a 

—  2 

maximum  half-power  bandwidth  of  1.27  *  10  cycles/mm. 

The  slight  variations  that  exist  among  the  spectra  in 
Fig.  7*2-1  relate  to  two  different  phenomena  which  will  pres¬ 
ently  be  discussed.  The  decrease  at  high  spatial  frequencies 
in  the  apparent  spectra  obtained  with  the  three  smaller  aper¬ 
tures  is  the  attenuation  due  to  the  phenomenon  called  vignetting 
(Born  and  Wolf  1964  ) .  Vignetting  occurs  when  the  diffraction 
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FILM  TYPE  *  ROYAL -X  PAN 
DEVELOPER  *  OK  90 

AVERAGE  AMPLITUDE  TRANSMITTANCE,  »0  =0  53 


FREQUENCY  (CYCLES/  mm) 


FIG.  7.2-1 


EFFECT  OF  APERTURE  SHAPES  ON  THE  MEASUREMENT 
OF  THE  POWER  SPECTRUM  OF  GRAIN  NOISE 


angles  are  sufficiently  large  to  force  some  rays  to  miss  the 
lens.  When  the  object  which  generates  the  angular  diffraction 
modes  is  placed  near  the  lens,  the  vignetting  is  reduced  or 
eliminated.  Slight  vignetting  is  also  present  in  the  apparent 
spectrum  furnished  by  the  largest  aperture.  In  this  case  the 
vignetting  appears  at  lower  frequencies  .  In  the  extreme  case 
when  the  aperture  is  as  large  as  the  lens,  vignetting  will  be 
present  for  all  spatial  frequencies  except  DC. 

The  second  discernible  difference  in  Fig.  7  2-1  is  the 
decrease  in  variability  that  results  when  the  largest  aperture  is 
used.  The  relation  between  the  aperture  size  and  the  variance 
of  the  estimated  power  spectrum  was  determined  for  rectangular 
and  diamond  apertures,  in  Chap.  5  under  the  assumptions  of  wide- 
sense  stationary,  Gaussian  statistics.  The  key  result  was  Eq . 

(5. 1-3)  which  can  be  rewritten  as 


var(  <r'uJ1(u,v ) )  t=W2(u,v)  +  F2(u,v)  +  tQ2H(u  ~)g(u,v) 

(7.2-1) 


where 

W(u,v) 

H(u, v ) 
F(u,v) 


apparent  spectrum 
W (u,v ) *  H(u,v) 

LW  sinc2uL  sinc°vW 

LW  f°°  fK  W(£,tj)  sinc(u  +  £  )L  sinc(u  -  %  )L 

—  00  —  or 


•  sine  (v  -  r|  )w  sinc(v  +  t)  )w  d^dr) 
=  average  amplitude  transmittance 


o 

G  (u,  v )  =  LW  /°°  /“  W  ( I* ,  T| )  •  |  sine  (u  +  5  )L  sinc(v  +  )w 

—00  —00 


£ 


+  sinc(u  -  £  )L  sinc(v  -  r^)W 


■] 


d^dr] 
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The  functions  H,  G  and  F  are  positive  for  all  u,  v, 

L,  W  and  have  envelopes  that  decrease  uniformly  with  in¬ 
creasing  aperture  size.  For  large  apertures,  all  the  terms 
except  for  Ws(u,v)  become  insignificant  away  from  the  origin. 
Hence,  on  the  basis  of  Eq.  (7.2-1),  less  variability  is  ex¬ 
pected  in  the  estimate  when  a  large  aperture  is  used  as  opposed 
to  a  smaller  one.  Figure  7.2-2  shows  more  explicitly  the  in¬ 
fluence  of  aperture  size  on  the  variability  of  the  measured 
spectrum.  The  extrapolated  spectrum  (dashed  lines)  corresponds 
to  the  region  where  the  zero-order  light  was  too  intense  for  an 
accurate  measurement  of  the  spectrum.  Fig.  7.2-2  also  illus¬ 
trates,  as  predicted  by  Eq.  (7.2-1)  that  the  envelope  of  the 
variance  decreases  with  increasing  frequency. 

The  result  of  Chap.  5  shows  that  for  a  rectangular  aper¬ 
ture,  the  ratio  of  the  variance  of  the  estimated  spectrum  to  its 
square  mean  cannot  be  made  less  than  unity.  However,  additional 
smoothing  is  provided  by  the  scanning  filter. 

The  action  of  the  scanning  filter  was  discussed  in 
Sec.  3-6  and  illustrated  in  Fig.  3.6-1  .  The  smoothing  of 
the  high  frequency  components  in  the  spectrum  by  the  scanning 
filter  is  shown  in  Fig.  7.2-3  .  In  this  experiment,  a  4"  diam¬ 
eter  lens  with  a  focal  length  of  f  =  .4  meters  was  used.  Both 
filters  shown  in  Fig.  7.2-3  had  diameters  well  within  the 
range  of  reasonable  values  established  by  the  criterion  of  Eq. 
(3-6-7)  which  for  this  case  was  5u  ^  d  <(  250u.  As  can  be  seen 
the  120|t  diameter  slit  furnished  a  uniformly  smoother  spectrum 
over  the  entire  frequency  band  of  interest  than  the  50' i.  slit. 
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SPECTRUM  (mm*)  »  O  -  SPECTRUM  (mm 


FIG  7.2-2  INFLUENCE  OF  APERTURE  ON  POWER  SPECTRUM  MEASURE¬ 
MENT  FILM  TYPE  IS  ROYAL-X  PAN  DEVELOPED  IN 
-148-  DK-50  SCANNING  FILTER  IS  50  /x  DIAMETER  PINHOLE 


SPECTRUM  (mm)  a>  SPECTRUM  (mm 


0  | - 1 - 1 - 1 - 1 - 1 

0  32  64  96  128  160 


A -006- 1-00-0022 


FIG.  7.2-3  INFLUENCE  OF  SCANNING  FILTER  ON  POWER  SPECTRUM 

MEASUREMENT.  FILM  TYPE  IS  ROYAL-X  PAN  DEVELOPED 

2 


IN  DK-50.  APERTURE  IS  60  x60  mnT  DIAMOND 
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7.3 


MEASUREMENTS  OF  THE  SPECTRA  OF  FILM  GRAIN  NOISE 


The  procedure  described  in  Sec.  6.2  was  used  to  obtain 
the  spectra  of  several  common,  commert  ially  available  films. 
The  important  characteristics  of  these  films  are  summarized  in 
Table  7.3-1.* 


TABLE  7.3-1 


ASA  Speed 


Resolving  Power 


Granularity  Class 


am 

800  (fast) 
Low 

<  55  cpm** 
Coarse 


2479 

400  (medium) 


Medium 
69-95  cpm 


Coarse-medium 


ROYAL  X-PAN 

1000-1200 
( fast ) 

Low 
55  cpm 

Coarse 


Special  Properties  Extended  sensitivity  Panchromatic 

to  red  light 

Thre«*  developers***  were  usedi  Kodak  DK-50,  D-19  and  Dektol . 


*  Data  obtained  from  Kodak  Tech  Bits,  1966,  No.  2.  Also 
Kodak  Reference  Handbook,  No.  1. 

**  Cycles  per  millimeter. 

***  The  developing  agent  is  a  chemical  compound  which  reduces 
silver  ions  to  metallic  silver.  The  basic  reaction  equa¬ 
tion  is  Ag+  +  e  — ►  Ag.  A  condition  that  any  developer 
must  fulfill  is  that  it  must  reduce  exposed  silver  halide 
grains  at  a  higher  rate  than  the  unexposed  grains.  Other 
conditions  that  a  developer  must  satisfy  are  (1)  it  must 
be  soluble  in  water  or  in  an  alkaline  solution,  (2)  be 
reasonably  stable  and  resistant  to  aerial  oxidation, 

(3)  yield  soluble  oxidation  products,  and  (4)  be  non- toxic. 
The  developer  makes  visible  the  latent  image  on  the  film 
after  exposure .  The  hypo  on  the  other  hand,  stops  the  de¬ 
velopment  process  and  fixes  permanently  the  visible  image. 

The  mechanism  of  the  formation  of  the  latent  image  was  ex¬ 
plained  by  the  theory  of  R.  w.  Gurney  and  N.  F.  Mott,  (1938) 
who  formulated  what  is  now  known  as  the  Gurney-Mott  principle 
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The  latter  is  what  is  known  as  a  paper  developer  and  is  nor¬ 
mally  not  used  with  the  processing  of  negatives.  D-19  is  a 
high  contrast  developer  usually  applied  for  10  minutes  at 
68°F .  DK-50  is  a  medium  grain  producing,  moderately  fast, 

very  popular  developer  used  in  the  same  way  at  D-19.  The 
Dektol  developer  was  applied  for  3  minutes  at  72°.  According 
to  Kodak  authorities,  the  use  of  Dektol  with  negatives  would 
result  in  grainy  transparencies.  It  was  of  interest  to  deter¬ 
mine  whether  the  diffraction  method  could  discern  the  difference 
between  film  treated  with  Dektol  and  a  recommended  developer. 

It  is  for  this  reason  that  Dektol  was  used. 

Figure  7*3-1  are  photomicrographs  of  the  grain  structure  of 
Kodak  2ky*  film  developed  in  D-19  at  three  different  levels  of  mag¬ 
nification.  To  obtain  these  photographs,  the  equipment  shown  in 
Fig.  6.1-10  was  used.  The  choice  of  film  is  incidental.  The  clump¬ 
ing  of  grains  is  evident  by  the  chain-like  aggregates  of  silver 
which  extend  in  some  cases  to  40  or  more  microns.  These  aggregates 
can  be  much  larger  than  the  actual  exposed  silver  halide  crystals. 

The  spectra  discussed  in  the  subsequent  paragraphs  were 
calibrated  by  using  Eq.  (6.2-10),  and  extrapolating  the  spec¬ 
trum  back  to  the  origin.  It  will  be  remembered  that  the  meas¬ 
urement  of  the  spectrum  near  the  origin  is  difficult  to  make 
because  of  the  presence  of  the  strong  zero-order  light.  Thus, 
the  extrapolated  spectrum  corresponds  to  the  region  where  the 
difference  between  the  zero-order  and  the  scattered  light  was 
less  than  10  db. 

Figures  7*3-2  to  7*3*5  are  actual  retraces  of  exper¬ 
imentally  obtained  spectra  of  the  films  listed  in  Table  (7. 3-1). 
Superimposed  on  each  graph  are  the  theoretical  spectra  predicted 
by  the  overlapping  circular  grain  (0CG),  the  random  checkerboard 
(RC)  and  the  modified  overlapping  circular  grain  (MOCG)  models. 

Each  film  sample  was  carefully  tested  with  a  microdensitometer 
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Fig.  7.3-1 
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Grain  Structure  of  Film  (Kodak  2496) 
Viewed  at  Three  Levels  of  Magnification 
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FIG.  7.3-2  POWER  SPECTRUM  OF  GRAIN  NOISE  OF  ROYAL-X  PAN 
COMPARED  WITH  THEORETICAL  MODELS 
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FIG.  7.3-3  POWER  SPECTRUM  OF  GRAIN  NOISE  OF  KODAK  2479  FILM 
COMPARED  WITH  THEORETICAL  MODELS 
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FIG  7  3-4  POWER  SPECTRUM  OF  GRAIN  NOISE  OF  KODAK  2475 
DEVELOPED  IN  D- 19,  COMPARED  WITH  THEORETICAL  MODELS 
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FIG. 7.3-5  POWER  SPECTRUM  OF  GRAIN  NOISE  OF  KODAK  2475 
(DEKTOL)  COMPARED  WITH  THEORETICAL  MODELS 
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to  insure  that  the  noise  was  isotropic.  The  theoretical  curves 
do  not  represent  "best  fit"  curves  in  any  sense.  The  parameters 
used  for  the  theoretical  models  were  obtained  by  independent 
measurements.  For  example,  in  order  to  specify  the  spectrum 
of  the  RC  model  as  given  by  Eq .  (4.2-3),  it  is  necessary  to 
specify  W  (0,0)  and  the  grain  size  Jf  .  Since  W^(0,0)  is 
estimated  by  extrapolating  the  experimental  spectrum  and  t 
is  measured  independently,  >  can  be  calculated  from 


f 


(7.3-1) 


By  this  procedure,  the  theoretical  spectra  and  the  extrapolated 
spectra  were  made  to  coincide  at  the  origin.  For  the  OCG  model, 
numerical  integration  of  Eq .  (4.2-15)  was  required.  From  Eqs . 
(4.2-14)  and  (4.2-15),  it  can  be  seen  that  the  power  spectrum 
depends  on  the  parameters  tQ,  h  and  j(q  where  JfQ  is  the 
mean  grain  radius  and  h  is  the  transmittance  of  the  gelatin 
(Sec.  4.2).  For  h  =  1  and  tQ  measured  independently  by  the 
procedure  in  Sec.  6.3  (part  c),  t  can  be  determined  by  nu¬ 
merically  computing  the  normalized  spectrum  Wn(0),  which 
does  not  depend  on  lQ,  and  using  Eq.  (4.2-16): 


o 


1 

87Tt  2 

O 


(7.3-2) 


With  f Q  computed  as  above,  the  power  spectrum  for  the  OCG 
model  can  be  computed  for  all  other  frequencies.  The  parame¬ 
ter  d,  which  is  the  average  number  of  grains  per  unit  area, 
can  be  determined  from  Eq.  (4.2-9). 

The  modified  overlapping  circular-grain  model  (MOCG  )  results 
when  the  effect  of  the  transmittance  of  the  gelatin,  h,  is  considered 
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(Savelli  1958).  A  very  good  estimation  of  h  can  be  made  by 
measuring  the  transmittance  of  unexposed  film,  i.e.,  film  freshly 
removed  from  the  pack  without  prior  exposure  to  light.  However, 
despite  precautions  to  the  contrary,  there  will  always  be  some 
unavoidable  exposure  partly  due  to  photon  leakage  into  the  emul¬ 
sion.  For  the  MOCG  model,  t  is  computed  from  Eq.  (7.3-2)  but 
Eq.  (4.2-13)  with  h  as  measured  by  the  procedure  of  Sec.  6.3 
(part  c)  is  used  for  computing  d. 

of  the  three  models ,  the  MOCG  gave  the  best  fit  to  the 
measured  spectrum,  particularly  at  lower  frequencies.  The  RC 
and  simplified  OCG  models  furnished  curves  which  were  not  sig¬ 
nificantly  different  even  though  the  latter  is  the  more  realistic. 
The  simplified  OCG  model  is  less  convenient  to  use  since  numeri¬ 
cal  integration  is  required  for  computing  its  power  spectrum. 

The  RC  spectrum  was  used  with  v  =  0  (Eq.  4.2-3)  which  corres¬ 
ponds  to  a  specific  grain  "alignment"  in  the  film  plane.  This 
rather  artificial  situation  is  a  consequence  of  trying  to  fit 
isotropic  models,  and  hence  the  problem  of  grain  alignment  is 
not  relevant.  Another  advantage  of  these  models  is  that,  un¬ 
like  the  RC  model,  they  do  not  predict  nulls  at  frequencies 
which  are  multiples  of  the  inverse  grain  sizes  . 

The  RC  and  OCG  models  yielded  power  spectra  which  were 
larger  than  the  MOCG  and  the  measured  spectra  at  most  frequencies  . 
This  behavior  is  a  consequence  of  the  larger  variance  of  these 
models  which,  because  the  variance  is  the  integral  of  the  power 
spectrum,  tended  to  produce  larger  spectra . 

Table  7*3-2  summarizes  some  of  the  principal  data  ex¬ 
tracted  from  Figs.  7*3-2  through  7*3“5  •  The  data  listed  in 
Table  7*3“2  are  in  good  agreement  with  the  film  properties  listed 
in  Table  7*3“1  •  The  known  coarsest  film  of  the  set  (Royal-X  Pan, 
when  the  other  films  are  used  with  their  normal  developers)  can 
be  seen  to  have  the  largest  grains  while  the  known  finest  film 
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of  the  set  (Kodak  2479)  had  the  smallest  grains  and  the  largest 
grain  density.  The  data  reveal  a  well-known  property  of  fast 
films  (Mees  and  James  1966),  namely, that  such  films  have  large 
grains  and  low  grain  densities. 

It  was  of  interest  to  determine  why  the  various  models 
furnished  different  results.  For  example,  for  every  film  and 
developer  combination,  the  estimated  grain  size  computed  from 
the  RC  model  is  smaller  than  the  corresponding  grain  sizes  com¬ 
puted  from  the  OCG  models.  Equally  apparent  is  the  discrepancy 
between  the  two  models  with  respect  to  the  computed  grain  density. 
In  every  case,  this  figure  is  largest  when  computed  by  the  OCG 
model  and  smallest  for  the  MOCG  model. 

The  theoretical  justification  for  these  differences  can 
be  found  by  considering  the  properties  of  the  RC  and  OCG  models. 
First,  in  all  three  models,  it  is  assumed  that  there  is  no  vari¬ 
ability  in  the  opacity  of  exposed  grains  and  that  an  exposed 
grain  has  zero  transmittance.  Thus,  in  the  OCG  model,  the  fact 
that  two  exposed  grains  may  overlap  does  not  further  reduce  the 
transmittance  of  a  point  in  the  overlapped  region.  Second,  since 
the  random  checkerboard  model  does  not  provide  for  overlapping, 
it  follows  that  the  total  opaque  area  is  simply  the  number  of  ex¬ 
posed  grains  multiplied  by  the  grain  size  -  a  result  not  true  for 
the  OCG  models.  For  the  latter,  the  developed  grains  do  not  re¬ 
move  clear  area  equal  to  their  size  because  some  grains  overlap. 
Hence,  for  a  given  transmittance,  the  nature  of  the  OCG  models 
requires  that  a  greater  amount  of  silver  be  deposited  than  in 
the  RC  model.  This  greater  silver  deposit  can  result  from  a 
greater  number  of  developed  grains,  larger  grains,  or  both. 

The  large  grain  areas  and  low  grain  densities  supplied  by 
the  MOCG  model  are  a  consequence  of  the  additional  parameter  h, 
which  must  be  determined  before  the  model  is  completely  specified. 
Comparing  Eq.  (4.2-1)  with  Eq.  (4.2-8)  it  is  seen  that  since  h 
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is  less  than  unity,  less  silver  (fewer  grains)  must  be  ex¬ 
posed  for  the  same  value  of  t  than  in  the  OCG  model.  This 
accounts  for  the  low  densities  recorded  for  the  MOCG  model. 

Experimental  evidence  is  given  later  in  this  section  to 

show  that  grain  sizes  decrease  . 'ith  decreasing  transmittance 

over  a  wide  transmittance  range.  Since  the  MOCG  model  can  be 

regarded  (except  for  a  multiplicative  constant)  as  the  ordinary 
.  •  •  t0 

OCG  model  with  equivalent  transmittance  t^  =  — and  since 

t^  tQ,  the  larger  grain  sizes  furnished  by  the  MOCG  model 

are  not  surprising  in  view  of  the  above. 

The  last  two  columns  of  Table  7»3~2  refer  to  data  obtained 
on  Kodak  2475  film  samples  treated  with  different  developers. 

The  recommended  developer  in  this  case  is  D-19  while  the  Dektol 
(deemed  unsuitable  for  2475)  was  expected  to  yield  a  grainy 
transparency.  The  results  corroborated  this.  Using  the  MOCG 
model  as  a  basis,  the  Dektol  developer  produced  grains  40$ 
larger  than  the  D-19.  With  the  OCG  and  RC  models,  the  increases 
were  equally  significant  at  46$. 

When  the  spectrum  of  grain  noise  is  plotted  vs  tQ,  the  RC 
model  predicts  a  maximum  in  the  vicinity  of  tQ  R,0.5  while  from 
Fig.  4.2-5  t  the  peak  for  the  OCG  is  shifted  to  the  left  near 
tQ  =0.57.  Figure  7.3-6  shows  the  variation  of  the  actual  spectrum 
as  a  function  of  tQ  for  several  different  frequencies.  In  each 
case,  the  peak  occurred  at  approximately  the  value  of  t  =  0.6, 
although  the  peak  value  decreased  with  increasing  frequency. 

This  is  predicted  by  both  the  OCG  and  RC  models  -  for  the  latter, 
at  least  until  the  first  null.  All  the  film  samples  used  in  this 
experiment  were  developed  in  Dektol  for  3  minutes  at  72°F.  Figure 
7.3-7  shows  the  variation  of  the  extrapolated  spectrum  with  t 
at  the  origin  |see  Fig.  4.2-5  •  J 

The  variation  of  grain  noise  with  transmittance  was  deter¬ 
mined  for  still  another  sequence  of  transparencies  which  were 
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FIG.  7.  3-6  VARIATION  OF  GRAIN  NOISE  WITH  AMPLITUDE 
TRANSMITTANCE  FOR  KODAK  2475 
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AVERAGE  AMPLITUDE  TRANSMITTANCE,  t0 


_ 

FIG.  7.3-7  VARIATION  OF  SPECTRUM  AT  DC  WITH  t 
(EXTRAPOLATED  FROM  MEASUREMENTS) 


obtained  by  imaging  Royal  X-Pan  grain  on  low  noise  Kodak 
2^96  film  with  magnification  of  l8x .  Although  these  samples 
were  generated  for  another  purpose,  their  use  in  this  case 
proved  advantageous  because  the  magnified  grains  furnished 
a  very  strong  scattered  light  level,  thus  enabling  the  observa¬ 
tion  of  the  spectrum  at  much  lower  frequencies  without  the 
need  for  extrapolation.  The  results  are  shown  in  Fig.  7*3-8  , 
where  it  can  be  seen  that  the  film  noise  again  peaked  near 
t  Z  0.6, although  at  16  cycles/mm  the  peak  was  found  to  be  quite 
broad. 

When  the  grain  sizes  were  computed  from  the  theoretical 
models  and  the  extrapolated  spectra,  it  was  found  that  the  aver¬ 
age  grain  area  decreased  with  decreasing  transmittance  until  a 
minimum  was  reached  beyond  which,  for  still  lower  transmittances , 
the  grain  size  increased  (Fig.  7.7-9  ).  The  decrease  in  grain 
size  with  decreased  transmittance,  that  is,  with  increased  ex¬ 
posure,  was  first  observed, but  not  explained, by  M.  Savelli 
(Savelli  1958). 

The  two  regions  on  either  side  of  the  minimum  reflect  separ¬ 
ate  phenomena.  The  negative-slope  region  can  be  explained  by  con¬ 
sidering  Silberstein' s  quantum  theory  of  exposure.  According  to 
Eq.  (4.1-1),  for  an  emulsion  of  equally  sensitive  grains,  the 
large  grains  will  have  a  higher  probability  of  absorbing  the  nec¬ 
essary  number  of  photons  required  for  development  than  the  small 
grains.  Hence,  for  low  exposure  levels,  the  large  grains  will 
develop  first  and  the  average  grain  size  will  be  high.  For  ex¬ 
ample,  if  m  and  e  are  unity  in  Eq.  (4.1-1),  the  probability 
of  exposure  is 

p(t=0)  =  1  -  exp 
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FILM  GRAIN  NOISE  —  ( db  BELOW  PEAK  INCIOENT  LIGHT  ) 


FIG.  7.  3-8  VARIATION  OF  GRAIN  NOISE  WITH  AMPLITUDE  TRANSMITTANCE ,  t0 
FOR  ROYAL  X-PAN  GRAIN  EXPOSED  ON  KODAK  2496  FILM 
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GRAIN  AREA  (MICRONS2) 


where  Ag  is  grain  area  and  E,p  is  the  incident  energy  density 
(the  other  terms  are  defined  in  Sec.  4.1).  This  equation  shows 
that  for  low  exposure  levels,  the  probability  of  exposure  varies 
approximately  linearly  with  area. 

As  the  exposure  becomes  more  intense,  the  probability  that 

smaller  grains  will  develop  increases  and  the  development  of 

large  grains  becomes  exhausted.  Thus,  the  average  grain  size 

is  reduced.  This  process  continues  until  significant  grain 

clumping  begins  to  form.  As  mentioned  in  Sec.  4,  these  clumps 

are  principally  due  to  the  phenomenon  of  grain  exposure  by 

contact  (Mees  and  James  1966)  and  are  not  accounted  for  in 

Silberstein' s  theory.  It  should  also  be  remembered  that  with 

increased  exposure,  development  spreads  to  the  grains  in  the 

deeper  layers  beneath  the  surface  so  that  clumping  occurs 

throughout  the  gelatin.  The  projection  of  these  clumps  on  a 

plane  produces  large,  contiguous  opaque  shadings  which  increase 

with  increasing  exposure.  These  shadings  are  the  equivalent 

grains  sensed  by  the  models .  It  is  of  interest  to  note  that  all 

the  theoretical  models  furnish  the  two  regions  reflecting  these 

separate  phenomena  and  that  the  minimum  for  all  three  occurs  at 

the  same  t  . 

o 

Small  grains  are  associated  with  high  resolving  power. 

For  Kodak  2475  developed  in  Dektol,  the  minimum  grain  area 
occurs  near  t  =0.25.  That  a  given  film  could  have  its  re¬ 
solving  power  depend  on  exposure  was  discovered  by  M.  Levy  (Mees 
and  James  1966).  Levy  found  that  for  certain  aerial  films,  the 
maximum  resolution  occurs  at  exposure  levels  corresponding  to  an 
average  amplitude  transmittance  of  t  ~0.34. 

One  advantage  in  making  spectral  density  measurements  by 
diffraction  is  the  possibility  of  obtaining  two-sided  spectra 
by  simply  scanning  in  the  back  focal  plane  on  either  side  of  the 
origin.  Such  two-sided  spectra  are  shown  in  Fig.  7.3-IO  ;the 
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film  in  this  case  was  magnified  Royal-X  Pan  imaged  on  2496 . 

The  unusual  structure  is  partly  a  reflection  of  the  imaging 
optics.  The  evenness  about  the  origin  and  the  repeatability 
along  orthogonal  axes  is  a  characteristic  of  isotropic  noise, 
which  requires  that  the  spectrum  be  an  even  and  symmetric  func¬ 
tion  of  the  orthogonal  spatial  frequencies  (u,v)  .  An  example 
of  nonisotropic  noise  is  the  condition  of  granularity  ’’streak¬ 
ing",  which  is  characterized  by  pronounced,  parallel  transmittance 
variations  in  some  specific  direction.  Streaking  can  be  de¬ 
tected  from  micrcdensitometer  traces.  Such  traces  and  the 
spectrum  of  the  resulting  noise  are  shown  in  Fig.  7.3-11  • 

7 .4  MEASUREMENTS  OF  THE  SPECTRA  OF  TOTAL  FILM  NOISE 

The  scattered  light  that  results  in  the  Fourier  plane 
when  film  is  not  immersed  in  a  liquid  gate  was  shown  earlier  to 
be  an  estimate  of  the  power  spectrum  of  the  process  |R(x,y)|  . 

The  random  variable  describing  the  complex  transmittance 
R (x ,y )  defined  in  terms  of  t(x,y)  and  V'(x,y)  by  Eq.  (4.3-1) 
is  reduced  to  the  amplitude  transmittance  t(x,y)  when  phase 
correction  is  provided  by  way  of  a  liquid  gate  matched  to  the 
refractive  index  of  the  film. 

When  the  magnitude  of  Rq  -  <(  R(x,y)  ">  was  plotted  vs 

the  average  amplitude  transmittance  tQ,  a  linear  relationship 
was  observed  over  wide  transmittance  variations,  excluding  very 
low  transmittances .  If  the  model  of  Sec.  4.3  is  applied,  namely, 
the  assumptions  of  independence  between  |V'(x,y)}  and  |t(x,y)| 
and  Gaussian  statistics  for  the  former,  then  the  theoretical  re¬ 
lationship  between  Rq  and  tQ  is  given  by  Rq  =  tQ  exp (-a^/2), 
which  is  a  line  with  constant  slope  that  furnishes  the  variance 
of  the  phase  noise  Sec.  6.3,  part  (d)J. 

i 

To  see  if  the  RQ  vs  tQ  relation  given  above  was  valid  for 
commercial  film,  two  sets  of  transparencies  were  tested:  Kodak 
2496  developed  in  D-19  and  Kodak  2475  developed  in  Dektol.  The 
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FIG.  7.3-11  POWER  SPECTRA  OF  NON-ISOTROPIC  GRAIN  NOISE. 
■170- 
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two  films  had  radically  different  properties  .  Kodak  2496  is 
a  high  resolution,  low  noise  film  with  a  "fine"  granularity 
rating  while  Kodak  2475*  particularly  when  treated  with  Dektol, 
is  a  fast,  coarse  film  (see  Tables  7.3-1  and  7*3“2).  When 
optical  path  lengths  were  considered, however ,  the  two  films 
were  more  similar  in  that  both  had  the  same  clear,  fast  drying 
backing  and  same  4-mil  gray  Estar  base.  The  results  are  shown 
in  Figs.  7.4-1  and  7.4-2  .  It  can  be  seen  that  the  two 
curves  have  nearly  the  same  slope  over  the  central  values  of 
tQ,  strongly  pointing  to  the  fact  that, in  this  case,  is 

essentially  a  measure  of  the  base  and  backing  and  not  of  de¬ 
veloper  and  granularity.  Table  7.4-1  summarizes  the  results 
of  this  experiment. 


TABLE  7.4-1 


FILM  AND 
DEVELOPER 

MEASURED 

a* 

RMS  THICKNESS 
DEVIATION 

RMS  OPTICAL  PATH 
LENGTH  DEVIATION 

2479 

D-19 

.836  rad 

I.38  x  10~^mm 

0.84  x 

-4 

10  mm 

2475 

Dektol 

.844  rad 

1.40  x  10~^mm 

0.85  x 

-4 

10  mm 

Measurements  indicated  that  the  total  noise  spectrum  was 
considerably  larger  at  low  frequencies  than  the  spectrum  of  grain 
noise  alone.  This  was  expected  after  obtaining  the  diffraction 
patterns  shown  in  Figs.  7.1-2  and  7.1-4  , where  it  can  be  seen 
that  the  effect  of  the  grain  noise  on  the  zero-order  light  was 
hardly  discernible  Fig.  (7.1-4  )J  while  the  effect  of  the  total 
noise  on  the  zero-order  jVig.  (7  . 1—1  )J  was  clearly  evident. 

*  For  film  in  a  vacuum  or  air. 
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AVERAGE  TRANSMITTANCE  |RC| 


FIG. 7.4 -I  EXPERIMENTAL  RELATIONSHIP  BETWEEN  THE  REAL  AND  COMPLEX 

TRANSMITTANCE  OF  FILM 
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AVERAGE  AMPLITUDE  TRANSMITTANCE,  t0 


A -006  -1-00-0053 


FIG.  7.4-2  EXPERIMENTAL  RELATIONSHIP  BETWEEN  REAL  AND 
COMPLEX  TRANSMITTANCE  OF  FILM 
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Figures  7  .4-3  through  7  4-5  are  experimentally  obtained 

spectra  of  total  film  noise.  Unlike  the  behavior  of  grain  noise. 

the  spectrum  of  the  total  noise  is  peaked  near  the  origin  and 

falls  off  very  rapidly  for  higher  frequencies.  In  Fig.  7.4-3 

spectral  densities  of  the  grain  and  to+-al  noise  are  compared. 

An  enlarged  view  of  the  grain  noise  is  provided  Calibration 

of  the  total  noise  spectrum  was  effected  by  using  Eq.  (6  2-10) 

with  I  R  |r  substituted  for  t  r.  It  can  be  seen  that  the 
o  o 

spectrum  of  the  total  noise  is  nearly  an  order  of  magnitude 
above  the  grain  noise.  The  large  total  noise  spectrum  enabled 
its  measurement  near  the  origin. and  little  interpolation  was 
required.  The  low-frequency  nature  of  the  total  noise  spec¬ 
trum  can  be  illustrated  by  observing  that  the  spectrum  decreases 
to  half  of  its  maximum  value,  which  is  at  the  origin,  at  a  frequency 
~  16  cycles/mm,  whereas  the  half-maximum  of  the  grain  noise 
does  not  appear  until  80  cycles /mm.  These  results  imply  that 
the  phase-related  process  J°(x,y)  (x»y^|  is  a  low-fre¬ 

quency  phenomenon . 

Equation  (4.3-14),  which  was  derived  on  the  basis  of 
the  previous  assumotions  of  independence  between  j  /,  (x,y)J 
and  jt(x,y)|  ,  and  the  normality  of  {^(x,y)(,  predicts  that 
the  spectrum  of  the  total  noise  will  vary  approximately  with 
tQfc  provided  that  tQ  is  not  too  small  and  that  the  re¬ 
lated  phase  spectrum  W0(w)  is  the  dominant  term.  The  spec¬ 
trum  in  Fig.  7.4-4  was  obtained  for  a  sample  of  Kodak  2475 
(Dektol)  with  transmittance  tQ  -  O.58.  Comparing  the  spectra 
of  Fig.  7*4-3  and  7.4-4  ,  it  is  seen  that  the  ratio  of  the 
square  of  the  transmi ttances  (1.24)  closely  follows  the  ratio 
of  the  peaks  of  total  noise  spectra  (1.29). 

Examination  of  Eq .  (4.3-14)  shows  that  the  spectrum  of 
the  total  noise  is  the  sum  of  attenuated  versions  of  the  spec¬ 
tra  of  grain  noise  and  related  phase  noise  and  their  convolu¬ 
tion  product.  Under  the  conditions  of  small  t  ,  large  a,,  and 

-174- 


i 


POWER  SPECTRUM  (mm 


80  X  10 


AVERAGE  AMPLITUDE  TRANSMITTANCE,  »0  *  .58 
AVERAGE  COMPLEX  TRANSMITTANCE  ,  |Rq|  *  .  34 

DEV  =  DEKTOL,  3  MIN.,  72  0  F 


\ 

\ 

60  -\ 


TOTAL  FILM  NOISE 


GRAIN  NOISE 


FREQUENCY  (CYCLES/mm) 


A- 006-  I  -  00-0076 


FIG.  7.4-4  POWER  SPECTRUM  OF  TOTAL  FILM  NOISE  AT 
PEAK  GRAIN  NOISE  (KODAK  2475) 
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sufficiently  high  frequency,  the  spectrum  of  the  total  noise 
may  have  a  smaller  amplitude  than  the  grain-noise  spectrum. 

The  situation  where  Wr(w)  <(  W^(w)  was  experimentally 
observed  (Fig.  7.4-5  )  for  a  Kodak  2475  sample  with  low  trans¬ 
mittance.  The  prediction  and  interpretation  of  the  results  in 
Fig.  7 .4-^  in  the  light  of  Eq.  (4.3-14)  is  however  not  com¬ 
pletely  valid  since  the  |  R  I  vs  t  curve  for  Kodak  2475  de- 

o  o 

veloped  in  Dektol  undergoes  a  departure  from  linearity  at  low 
transmittances ,  thus  violating  the  assumption  of  complete  in¬ 
dependence  between  |^(x,y)}  and  {t(x,y)f  upon  which 
Eq .  (4.3-14)  was  derived. 

7  .5  SPECTRA  OF  MAGNIFIED  GRAIN  NOISE  AND  FREQUENCY 

RESPONSE  OF  OPTICAL  SYSTEMS 

The  discussion  in  Sec.  6.5  indicates  that,  based  on  the 
theoretical  considerations  leading  up  to  Eqs .  (6.5-2)  and  (6.5~7)> 
magnification  and  imaging  techniques  can  be  used  to  obtain  the 
extended  spectra  of  grain  noise  and  the  frequency  response  of 
optical  systems.  Because  of  a  limitation  in  the  availability 
of  laboratory  apparatus,  the  results  in  this  section  were  ob¬ 
tained  using  procedures  which  were  not  ideal.  Hence,  these  re¬ 
sults  are  intended  more  to  illustrate  feasibility  than  to  pro¬ 
vide  refined  data. 

Ideally,  imaging,  magnification , and  spectral  analysis 
should  be  performed  simultaneously  as  spatially  sequential  oper¬ 
ations  on  the  optical  bench  comprising  the  coherent  optical  sys¬ 
tem.  The  available  laboratory  equipment  did  not  permit  this, 
and  it  was  necessary  to  adopt  a  procedure  which  involved  doing 
the  imaging  on  a  separate,  incoherent  system,  storing  the  image 
on  low-noise  film, and  then  transferring  the  film. to  the  coher¬ 
ent  optical  system  to  measure  the  spectrum.  A  major  disadvantage 
of  storing  the  image  on  film  results  from  the  interaction  of  the 
fi'i.m' s  own  granularity  with  the  imaged  granularity.  In  view  of 
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FIG.  7.4-5  ANOMALOUS  BEHAVIOR  OF  DENSE  FILM  (KODAK  2475) 
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the  analysis  given  in  Appendix  A,  it  is  unclear  whether  or  not  the 
spectrum  of  the  granularity  simply  adds  to  the  spectrum  of  the 
granularity  of  the  storage  medium  as  some  authors  propose 
(Doerner  1962).  However,  in  the  present  case,  the  image  was 
stored  on  Kodak  2496  film  which , insured  that  the  magnified 
granularity  would  exercise  a  dominant  role.  Kodak  2496  is  a 
high-resolution  (136-225  cycles/mm),  low-noise  film  suitable 
for  applications  where  grain  noise  has  to  be  kept  small .  A 
second  serious  disadvantage  results  from  the  use  of  incoherent 
illumination  systems.  This  is  discussed  in  greater  detail  below. 
The  reader  is  reminded,  however,  that  these  problems  are  a  con¬ 
sequence  of  the  procedure  that  was  used  and  could  be  eliminated 
by  using  an  all-coherent  optical  system. 

In  order  to  determine  experimentally  the  frequency  re¬ 
sponse  of  the  "Raptar"  commercial  enlarger  lens  described  in 
Sec.  6.1,  the  following  procedure  was  used.  The  log  spectrum 
of  the  film  storing  the  imaged  magnified  granularity  of  Royal 
X-Pan  was  superimposed  on  the  original  log  spectrum  of  the 
Royal  X-Pan  after  calibrating  the  frequency  axis  according  to 
Eq.  (6.5-2),  that  is,  the  frequencies  associated  with  the  mag¬ 
nified  spectrum  were  ir  iltiplied  by  the  magnification  factor. 

The  log  spectra  were  made  to  coincide  at  the  origin  ,  there¬ 
by  normalizing  the  frequency  response  with  respect  to  its 
value  at  the  origin.  Away  from  the  origin,  a  separation 
that  increased  with  frequency  was  observed  between  the  log 
spectra.  The  interpretation  of  this  separation,  in  the  light 
of  the  imaging  theory  summarized  in  Sec.  6.5,  is  as  follows: 

The  imaging  of  granularity  in  incoherent  light  produced 
an  intensity  in  the  image  plane  given  by  (Eq.  6.5-6) 

Ip(x,y)  -  f  r  T(S,ri)  B!  (x-$,  y-rl)  <*£dr| 

*  -or-  —00 
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When  film  is  placed  in  the  image  plane,  it  is  possible  by  proper 
control  of  exposure  and  development  to  generate  an  amplitude 
transmittance  variation  proportional  to  the  light  intensity. 
Hence , 


t(p)(x,y)  =  K  Ip(x,y)  (7.5-1) 

where  K  is  the  constant  of  proportionality.  It  follows  from 
Eqs .  (7. 5”1)  and  (6,5-6)  that 

^P)(a,p)  =  Kp  .\p(a  ,p  j  *  G  (a  )  (7.5-2) 

where  G(a,i  )  is  the  self-correlation  function  of  the  trans¬ 
mission  function  B.j.(x,y),  and  T  and  are  the  autocor¬ 

relation  functions  of  the  intensity  and  amplitude  transmissions. 
The  Fourier  transform  of  Eq .  (7.5-2)  yields 

Wt'P\u,v)  K  WT (u, v )  I  b;[(u,v'!  I  r'  (7.5-3) 

which  establishes  the  proper  relationship  between  the  spectrum 
of  the  print  and  the  spectrum  of  the  original  sample.  However, 
knowledge  of  the  granularity  structure  of  the  original  sample 
was  in  the  form  of  Wt(u,v)  and  not  WT(u,v).  Hence,  the  ob¬ 
served  difference  between  the  two  log  spectra  can  be  interpreted 
as  log  |  bI(u,v)  |  provided  that  any  point  on  the  film  is 
in  a  dichotomic  state  -  either  opaque  or  transparent.  This  con¬ 
dition  implies  that  the  amplitude  and  intensity  spectra  are  the 
same.  Figure  7.5-1  shows  the  normalized  average  frequency 
response  of  the  lens  obtained  by  this  procedure.  The  term  aver¬ 
age  is  used  because  (a)  white  light  was  used  for  illumination 
during  imaging  and  (b)  most  of  the  lens  surface,  as  opposed  to 
the  axial  region  alone,  was  used  for  image  formation. 
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FIG  7.5-1  MEASUREO  AVERAGE  FREQUENCY  RESPONSE  OF  WOLLENSAK 
"RAPTAR"  ENLARGING  LENS  USED  IN  WHITE  LIGHT 


-1P1 


When  it  is  desired  to  obtain  the  extended-frequency  spec¬ 
trum  accurately  by  magnification,  the  imaging  optics  must  have 
a  frequency  response  which  is  flat  over  the  frequency  interval 
of  interest.  The  "Fluorestar"  microscope,  described  in  Sec. 

6.1  was  used  for  high-resolution  imaging.  In  this  case,  the 
attenuation  of  the  microscope  optics  was  negligible  as  can  be 
seen  in  Fig.  7.5-2  .  The  cut-off  frequency  (the  maximum 
spatial  frequency  transmitted  by  the  lens )  was  I85O  cycles  per 
millimeter  corresponding  to  an  N.A.  (for  definition  see  Sec. 

6.5)  of  0.5  and  a  wavelength  of  5400  A.  The  magnification  was 
measured  by  a  precision  stage  micrometer  and  was  found  to  be 
130X.  The  film  type  and  developer  was  Kodak  2479  developed  in 
D-19.  When  the  original  spectrum  was  superimposed  on  the  spec¬ 
trum  of  the  magnified  grains  so  that  the  two  coincided  at  the 
origin,  excellent  coincidence  resulted  at  low  frequencies  al¬ 
though  a  slight  drooping-off  occurred  in  the  spectrum  of  the 
unmagnified  grains  at  higher  frequencies  due  to  off-axis 
vignetting.  The  spectrum  of  the  magnified  grains  is  not  af¬ 
fected  by  off-axis  attenuation  because  the  actual  diffraction 
angles  are  very  small.  For  example,  at  a  spatial  frequency  of 
I.58  cycles  per  mm.,  corresponding  to  a  calibrated  frequency 
of  204  cycles  per  mm.,  the  diffraction  angle  is  I.58  x  .6328  x 
10~  -  10  radians  =  .057  degrees.  Hence,  magnification  offers 

the  advantage  of  observing  the  spectrum  at  high  spatial  fre¬ 
quencies  without  exceeding  the  off-axis  capability  of  most  lenses. 
The  extended  spectrum  shown  in  Fig.  7-5-2  shows  that  at  fre¬ 
quencies  above  -  150  cycles/mm,  the  spectrum  of  grain  noise  can 
be  approximated  with  high  accuracy  by  a  negative  exponential. 

In  this  particular  case,  the  high-frequency  spectrum  is  propor¬ 
tional  to  exp(-l6.2  x  10~3  ] w |  ). 
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FIG.  7. 5-2  SPECTRUM  OF  GRANULARITY  MAGNIFIED  I30X 

COMPARED  WITH  SPECTRUM  OF  ORIGINAL  GRANULARITY 
(KODAK  2479  DEVELOPED  IN  D 19) 

. 


-I83- 


APPENDIX  A 


FILM  GRAIN  NOISE  IN  THE  PRESENCE 
OF  SIGNAL 


This  appendix  represents  an  attempt  to  determine  the 
nature  of  film  grain  noise  in  the  presence  of  signal.  Using 
the  random  checkerboard  model  (chap.  4,  Sec.  2)  as  a  basis, 
it  is  shown  that  grain  noise  in  the  presence  of  signal  cannot 
be  considered  as  additive  or  multiplicative  noise. 

The  transmittance  fluctuation  about  the  mean  transmit¬ 
tance  t  is  defined  by 
o 

T(x,y)  -  t(x.y)  -  tQ 

where  t(x,y)  is  a  random  variable  denoting  the  transmittance 
at  (x,y) .  It  is  assumed  that  a  developed  grain  is  opaque  to 
light  [t(x,y)  =  0]  while  an  undeveloped  grain  is  completely 
transparent  [t(x,y)  =  1].  The  relative  phase  of  the  trans¬ 
mitted  light  is  undisturbed  in  passing  through  the  medium. 

The  correlation  function  for  this  model  is  given  by  Eq .  (4.2-2) 


4 


(x  ,x  ,y  ,y  )  =  0  (x  -x  ,y  -y  ) 
tt  i  2  1?  tt  2  12 


X  -X 


I  y  -y  | 

-  V1  -  ‘o'  u  -  f1  -  L~ViL> 


x2-xi!  >  ly2-Yi  I  <  * 


=  0  otherwise 


(A-l) 
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Assume  now  that  a  bipolar  transmittance  representing 
the  signal  s(x,y)  is  present.  If  all  correlation  intervals 
and  observations  are  restricted  to  the  x-axis  (y  =  0),  then 
the  analysis  can  be  made  conveniently  one-dimensional  with 
0  h, 0)  0  (a),  s(x,0)  r  s(x)  and  t(x,0)  e  t(x).  If  the 

TT  "T 

operating  point  is  confined  in  the  linear  region  of  the  trans¬ 
mittance-exposure  characteristic  of  the  film, then  the  ensemble 
average  of  the  transmittance  at  (x,0)  is  s(x)  +  t  .  De¬ 
noting  by  p^j  the  probability  that  t(x1)  =  i  and  t(x2)  =  j , 
it  folio  s  that 

1  1 

3  <t<xi>  t(x^’  'SS  ij  PiJ  =  Pl1 


where  (  0  denotes  ensemble  average. 

If  B  is  a  random  variable  that  denotes  the  coordinate 
of  the  start  of  the  cell  (Fig.  A-l  )  that  contains  xx  (x2  >  x1) 
and  B  is  assumed  uniformly  distributed  in  the  interval  (Xl  -  {jX;L) 
then  the  distribution  function  of  B  is 


where 


Ffi(x)  =  0 
Fb(x)  =  f 

fb(x)  "  1 


x:-(, 

T 


X  <  xx  -  l 

<  X  <  xx 
X  0  x: 

(A-3) 


fb<x> 


Prob 


x 


} 


Now 
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FIG.  A-l  THEORETICAL  POINT  APERTURE  MICRODENSITOMETER  TRACE  OF 
RANDOM  CHECKERBOARD  MODEL  OF  FILM  GRAINS 
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^tt^Xl,X2^  =  Pr°b  {t^X2^  =  1  *  fc^Xi)  =  1}*  Pr°b  {t(x.)  =  X} 

=  (Prob  |x2  <  B  +  l  |  +  Prob  jx2  >  B  +  f  ^ 

•  Prob  |t(x2)  B  lj  )  Prob  |t(x1)  =  1  } 

=  (1  -  Fg  (x2--t)  +  Fg  (x2-( )  ( tQ+s  (x2  ) ) 

*  (t0+s(xi)  ) 

Hence  from  Eq.  (A-3)  it  follows  that  for  xg  x1 

^tt(x1,X2)  =  (to+s(xJ)  ^1  -  -  +  ‘■~1  1  (t0+s(xa))j 


0  <  X2“X1  <  l 


=  (to+s(xl))  (t0+s(xa)) 


Vxi  >  l 


(A-4) 


The  autocorrelation  function  for  both  cases,  x  >  x  and 

3  2  1 

x  ,  can  be  conveniently  written  as 


/  \  /  I 

ftQ  +  s(min(x1,x2) )  I  ll - - 

xx-x2))) 


x2“xl 


(t  +  s(max(: 


I  Xg-xJ  <  / 

=  (tc  +  s(x1))  (tQ  +  s(x2))  |  x2-xj  >  / 

(A-5) 


. 
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Consider  now  the  following  mathematical  model  which  in¬ 
cludes  both  additive  and  multiplicative  noise  in  a  representation 
of  transmittance  in  the  presence  of  signal: 

t*  (x)  =  tQ  +  s  ( X  )  +  U.r(x)  +  -\s(x)t'(x) 

(A-6) 

where  a  and  A  are  adjustable  constants  and  t(x)  and 
t' (x)  are  zero-mean  signal-uncorrelated  additive  and  multipli¬ 
cative  noise  components  respectively.  Then 

<t' (x)  >  -  to  +  S (x )  (A-7 ) 

0t,t,  =  to2  +  +  s(x2))  +  s(xi)s(x2) 

+  H20tt(x1,x2)  +  iA0t,  v(x1,x2)s(x1) 

+  u'\s(x2)0TT,  (xx,x2)  +  A2s(xi)s(x2)0t,t,  (xx,x2) 

(A-8) 

where 

^>t,t(x1,X2)  =  <T'(x1)  T  (x2)  > 

0  ,  (*  ,X  )  =  <t(x  )  t'  (x  )  > 

TT1  12  1  2 

0  T  i  j  i  (X!»X2)  =  <T'  (Xx)  T'  (X2)  > 

In  order  for  /tt(xi,x2)  of  Eq.  (A-5)  to  be  the  correlation 
function  for  the  process  t' (x)  ,  it  is  required  to  find  a 

set  of  values  for  u  and  ^  such  that  0..(x.x  )  =  0. ,.,(x  .x  ). 
In  looking  for  such  numbers,  the  following  restrictions  that 
apply  to  t’ (x)  must  not  be  violated: 
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(i )  t ' (x )  must  be  r 

(ii)  \,r,(o,o)yo 

(iii)  tTT(xl,x9)  =  tc 

Restrictions  (i)  -  (iii)  preclude  finding  real  values 
of  1 1  and  >  such  that  (x  ,x  )  =  (x  ,x  ).  For  ex- 

ample  in  attempting  to  equate  Eq .  (A-8)  with  Eq .  (A-4  ) ,  it  is 
found  that  \s  =  -1,  thus  violating  (i).  Hence  to  the  extent 
that  the  random  checkerboard  model  represents  film  granularity, 
film-grain  noise  cannot  be  resolved  into  signal-uncorrelated 
independent  additive  and  multiplicative  noise  processes. 


sal 

(1 


•  v  (■  • 
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APPENDIX  B 


SUPPRESSION  OF  THE  ZERO-  ORDER 
LIGHT  AT  DC 

The  contoured  apertures  discussed  in  Sec.  6.4  are  not 
adequate  for  measuring  the  power  spectrum  at  zero  spatial  fre¬ 
quency  because  their  diffraction  patterns  do  not  have  nulls  at 
the  origin.  A  diffraction  pattern  with  a  null  at  DC  can  be 
obtained  with  an  aperture  function  of  odd  symmetry  along  one 
or  more  axes.  For  example,  consider  a  one-dimensional  system 
with  odd  aperture  function  g(x).  Let  the  aperture  width  be 
D.  Then  the  diffraction  amplitude  is  proportional  to 


E(u)  =  fD/2  g(x) 

-D/2 


e-ig7:uxdx  =  -j  fD/2  g (x ) 
-D/2 


sin  2rrux  dx 

(B-l) 


which  is  zero  at  the  origin. 

One  of  the  simplest  aperture  functions  which  has  this 
property  is  the  half-wave  plate  shown  in  Fig.  B-l  .  The  half¬ 
wave  plate  is  a  one-cycle  phase  grating  and  produces  a  relative 
diffraction  intensity  given  by 


Hh(u,v)  =  LW  sine2  ^  sine2  vW  sin2  (B-2) 

which  is  shown  in  Fig.  B-2  .  The  sidelobe  response  is  shown 
in  Fig.  B-3  .  When  the  half-wave  plate  is  used  in  spectral 
density  measurements,  there  is  a  loss  in  resolution  in  the 
measurement  of  the  spectrum  along  the  frequency  axis  parallel 
to  the  direction  of  the  phase  grating.  However,  the  loss  of 
resolution  will  generally  not  be  serious  enough  to  affect  the 
measurement  of  ordinary  film  noise. 
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AMPLITUDE 

TRANSMITTANCE 


A-  006-1 -00 -001 1 


FIG.  B-l  HALF-WAVE  RETARDATION  PLATE  :  AMPLITUDE  AND  PHASE  CHARAC¬ 
TERISTICS 
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FIG. 


-Z  DIFFRACTION  PATTERN  OF  HALF-WAVE  RETARDATION  PLATE. 

(THEORETICAL) 
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FIG.B-3  LOG  DIFFRACTION  PATTERN  OF  HALF-WAVE  RETARDATION  PLATE. 

(THEORETICAL) 


When  the  half-wave  plate  is  used,  the  apparent  spectrum 
W(u,v)  is  given  by 


W (u, v ) 


-  —  /W^2  /W/2  dydrj  /L^2  /L/2  f  (x-$  ,y-r]  ,u ,  v  )dxd£ 

-W/2  -W/2  0  G 


+  /°  /°  f  (x-^y-n^v^xd^ 

-L/2  -L/2 


"  f°  .  fL^2  f (x-?,,y-q,u,v)dxd^ 
-L/2  0 


-  /L/2  /°  f (x-f,,y-r1,u,v)dxd^ 
0  -L/2 


(B-3) 


where 


f(x-e,Y-n,u,v)  -  Mx-^y-a)  <x'5  >  +  vh'"l)l 


With  the  transformation 


a  =  x-« 

P  =  y-i 

7  =  x 

<5  =  y 


Eq.  (B-3)  may  be  written  as 


w(u,v)  =  r  r  *(a,P)  e-^77 (ua+vfs ) 


—  fr  —nr 


do  dp 


(B-4) 


where  hH(a,P)  =  hH(n,o)  ^(o^P)  and  is  shown  in  Fig.  B-4 
From  Eq.  (b-4  ) ,  the  apparent  covariance  £(ct,p)  is 


(a,p)  -  V“  ,P)  <Ma,p) 


(B-5) 
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FIG.  B-4  SELF-CORRELATION  FUNCTION  OF  HALF-WAVE  PLATE  OF  AREA 

A  =  LW  . 

i 
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where  (a,p  )  is  the  true  covariance  of  the  process.  It  is 
evident  that  for  a  square  aperture  with  side  D,  there  will  be 
greater  distortion  of  the  covariance  for  correlation  intervals 
in  the  direction  of  the  phase  grating,  due  to  the  steeper  slope 
and  phase  reversal  property  of  ht,(n,0''.  This  distortion,  when 
viewed  in  the  frequency  domain,  is  equivalent  to  averaging 

over  frequency  with  weights  proportional  to  H^u-CjV-'j) 

-  D  sine"  (u-f  )d/2  sinc"(v-))D  sin"  ;-(u-*)d/2.  Observe  that  the 
weight  given  to  W ( ‘ ,  *) )  at  £  =  u  is  zero,  which  is  quite  the 
opposite  of  what  the  ideal  aperture,  which  would  give  infinite 
weight,  would  do.  The  significant  averaging  over  frequency  ex¬ 
tends  over  the  two  lobes  on  either  side  of  *  -  u  which  furnish 
a  null-to-null  width  in  frequency  of  4d-1.  For  an  ordinary 
rectangular  aperture  of  the  same  width,  the  significant  aver¬ 
aging  over  frequency  is  restricted  to  the  main  lobe  of  the 
diffraction  pattern  which  has  a  null-to-null  width  of  2D  1 

The  loss  of  resolution  inherent  with  the  use  of  the 
half-wave  plate  is  not  serious  provided  that  D/2  is  much 
larger  than  the  largest  significant  correlation  interval  of 
interest.  In  the  measurement  of  ordinary  film  noise,  the  cor¬ 
relation  intervals  are, at  most,  of  the  order  of  hundreds  of 
microns  while  D  may  be  several  centimeters. 

From  a  theoretical  viewpoint,  the  half-wave  plate  is  a 
solution  to  the  problem  of  zero-order  elimination  at  the  origin. 
However,  from  a  practical  viewpoint  there  are  very  serious  de¬ 
sign  problems  to  overcome.  For  example,  if  the  phase  difference 
is  accomplished  by  using  clear  glass  plates  of  different  thick¬ 
nesses,  the  required  thickness  difference  between  the  two  halves 
is  >/2(n-l)  which,  for  >  =  6328  A  and  a  typical  value  of 
n  ~  1.6,  is  ~  5*?  x  10-'1  mm  -  a  value  difficult  to  maintain  to 
high  tolerance  over  a  large  area.  The  tolerance  and  parallelism 
requirements  to  achieve  a  deep  null  are  very  severe.  Even  if  such 
a  null  were  furnished  there  would  remain  the  difficulty  of  measuring 
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W'OjO)  with  a  very  narrow  pinhole  or  slit  in  order  not  to  lose 
the  null  by  averaging  over  the  sides  of  the  adjacent  zero-order 
lobes.  That  is,  the  width  of  the  slit  must  be  confined  to  fre¬ 
quencies  for  which  the  zero-order  light  is  essentially  insignifi¬ 
cant  compared  with  the  power  spectrum.  Only  in  this  way  will  the 
null  be  preserved.  However,  this  leads  to  a  third  problem; 

Since  the  averaging  effect  of  the  slit  is  confined  to  a  very 
narrow  band  of  frequencies,  the  estimate  of  the  power  spectrum 
will  be  inherently  unreliable  because  of  very  little  smoothing  over 
frequency  .  Smoothing  could  be  improved  by  increasing  the  length 
of  the  slit,  but  this  would  cause  loss  of  resolution. 

An  interesting  technique  for  achieving  the  requited  phase 
shift  between  the  two  halves  of  the  aperture  while  avoiding  the 
problems  associated  with  the  use  of  plates  of  different  thick¬ 
nesses  is  to  use  two  plates  of  equal  thickness  (i.e.,  a  single 
plate  cut  in  two ^  and  then  rotate  one  plate  with  respect  to  the 
other.  Such  a  device  was  built  in  the  laboratory  specifically 
for  zero-order  cancellation  in  power  spectra  measurements  and 
is  shown  in  Fig.  B-5  .  The  principle  of  operation  can  be 
explained  as  follows.  The  lower  plate  is  fixed  and  positioned 
for  normal  incidence.  Rotation  of  the  upper  plate  about  an  axis 
orthogonal  to  the  optic  axis  changes  the  angle  of  incidence  for 
the  collimated  beam  and  hence  the  optical  path  length.  It  is 
not  hard  to  show  that  in  order  to  obtain  a  180°  phase  difference 
between  the  waves  of  the  upper  and  lower  plates,  the  upper  plate 
must  be  rotated  through  an  angle  A  where  a  can  be  obtained 
from 


2d 


n 


-  cos(a-a') 
COS  A  ' 


n+1 
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Laboratory  Built  Half-Wave  Phase  Grating 


where 


■=>'  .  sin'1  (  £±n«  ) 

n  =  refractive  index  of  the  glass  plate 
d  -  thickness  of  the  plate 
\  =  wavelength  of  radiation 


COS  ft 


When  the  angles  are  small,  as  they  usually  are, 
~  1  -  rr /2l  so  that 


] 


sin  -  ft, 


(B-6) 


For  the  grating  used,  d  =  1.8  cm,  n  =  1.515  and 
0  "0.7  deg.  This  order  of  sensitivity  was  easily  controlled 
with  a  micrometer  as  shown. 

The  diffraction  pattern  furnished  by  the  actual  half-wave 
plate  (Fig.  B-6  )  corresponded  closely  to  the  theoretical  dif¬ 
fraction  pattern  shown  in  Fig.  B-2  .  The  two-dimensional  dif¬ 
fraction  pattern  was  observed  in  the  Fourier  plane  and  recorded 
on  film  through  the  use  of  a  microscope  objective  which  imaged 
the  magnified  diffraction  pattern  on  the  back-plate  of  a  camera. 
Figures  B-7  and  B-8  show  the  diffraction  patterns  for  square 
and  diamond  apertures;  the  null  at  the  center  is  clearly  visible 
and  its  orientation  can  be  seen  to  be  unaffected  by  a  rotation 
of  the  aperture.  The  orientation  of  the  null  depends  only  on  the 
positioning  of  the  glass  plates  in  the  collimated  region. 

To  obtain  the  DC  value  of  the  spectrum  from  the  diffrac¬ 
tion  pattern  furnished  by  the  half-wave  plate,  it  is  not  possible 
to  use  the  formula  given  in  Eq.  (6.2-10)  because  in  the  present 
case  the  slit  is  not  used  to  smooth  the  zero-order  light;  in  fact, 
in  order  to  take  advantage  of  the  null  at  the  origin,  the  equiv¬ 
alent  bandwidth  of  the  slit  in  the  direction  parallel  to  scan 
must  be  much  smaller  than  2(ttD)_1  where  D  is  the  length  of  a 
side  of  the  square  aperture.  In  other  words  the  slit  is  required 
to  distort  the  diffraction  pattern  as  little  as  possible. 
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FIG.  B-6  EXPERIMENTALLY  OBTAINED  DIFFRACTION  PATTERN  OF 
HALF-WAVE  PLATE  (COMPARE  WITH  FIG.  B-2) 


A)  Diffraction  Pattern  of  10  x  10  nun2  Square  Aperture. 
Width  of  Main  Lobe  is  127  Microns.  Diffraction  Pat¬ 
tern  Corresponds  to  sinc?uD  sinc2vD. 


B)  Diffraction  Pattern  of  10  x  10  mm2  Square  Aperture. 
Width  of  Main  Lobe  is  127  Microns.  Diffraction  Pat¬ 
tern  Corresponds  to  sinc2ud  sinc2vD  with  Phase  Grating 
in  Place.  Null  is  Clearly  Seen  at  Center. 


Fig.  B— 7  Fourier  Plane  Photographs  Illustrating 
Use  of  Phase  Grating  with  a  Square 
Aperture 
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A)  Diffraction  Pattern  of  10  x  10  mm2 
Diamond  Aperture 


B)  Diffraction  Pattern  of  10  x  10  mm? 
Diamond  Aperture  with  Phase  Grating 
in  Place . 


Fig.  B-8  Fourier  Plane  Photographs  Illustrating 
Use  of  Phase  Grating  with  a  Diamond 
Aperture 
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The  calibration  technique  for  extracting  the  DC  value 
of  the  spectrum  can  be  derived  as  follows. 

From  Eqs .  (B-2)  and  (6.2-5).  the  x-y  recorder  when 
operating  in  the  log  mode  will  produce  a  record  given  by 

Ri  =  10  lo<?  (Af)V1  (u>v)*9 (u> v)  (B-7) 

where 

Tv  =  6'j28  A 

f  =  focal  length  of  lens 

A  =  slit  area 

s 

H„(u,v)  is  given  by  Eq.  (B-2) 

ri 

g(u, v)  is  given  by  Eq.  (3.6-6), 

when  the  system  is  devoid  of  film.  When  the  film  sample  is 
placed  in  the  collimated  region,  the  recorder  will  then  produce 
a  record,  Rg ,  given  by 

Rg  =  in  log  (7vf)2As-1  g(u,v)*  ^Hn(u,v)  +  WT(u,v) 

where  W^(u,v)  is  the  apparent  spectrum.  If  the  y  axis  of 
the  recorder  is  displaced  by  an  amount  K  so  that  the  two 
records  are  made  to  coincide  at  the  peaks  of  the  diffraction 
patterns  (i.e.  ,  v  =  0  u  =*  .  75  1  =  uQ)  ,  then  .the  mathematical 
equivalent  of  this  procedure  can  be  written  as 

10  log  (M)2As-1  g(u,v)*jt2HH(u,v)+wT(u,v)j  +  K 

=  10  log  (>f) 2As"1HH(u,v)*g(u,v)  u  =  uq 

v  =  0 

Except  when  t2  is  much  smaller  than  unity  it  can  be 
assumed  with  little  error  that 
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tJHH(u>v)#g(u>v)  »  WT(u,v)*g(u,v)  u  =  uq  (B-P) 

v  =  0 

—  2 

so  that  to  a  good  approximation  K  =  log  t  .  Hence  super¬ 
posing  the  two  diffraction  patterns  is  equivalent  to  compensating 
for  the  nonzero  density  of  the  film.  If  the  film  were  noise¬ 
less,  the  compensation  for  t2  would  enable  the  two  diffraction 
patterns  to  overlap  completely.  When  there  is  no  overlap,  the 
difference  between  the  diffraction  patterns  after  compensation 
is  a  measure  of  the  film  noise. 

Assume  that  at  the  origin  there  exists  a  displacement 
A  between  the  two  records.  Then  A  is  given  by 

A  =  10  log  (>«f )2Ag"1g(u,v)‘»^HH(u,v)  +  t~?WT(u,v)^ 

-  10  log  (Af)2Ag-1g(u,v)HH(u,v)*g(u,v)  u  _  Q  (B-8) 

v  =  0 

If  the  slit  is  properly  chosen  and  the  aperture  is  sufficiently 
large,  then  to  a  good  approximation 


(M£~ 

-JT—  WT(u,v)*g(u,v)  I  Q  *  Wt(0,0) 
s  I  v  =  0 


and,  from  Eq.  (B-8) 


t2  ( A  f)  2  /  4-  \ 

Wt(0,0)«-5S -  ^10  iU-l J  HH(u,v)^g(u,v)  |u  = 


u  =  0 
v  =  0 


(B-9) 


If  the  slit  is  rectangular  with  area  A  =  £d  with  l  »  (Af)(7TD) 

s 

and  f  y>  d,  then 
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h(u,  v)  *g(u, v) 


|u=  0 
v  =  0 


1 


where 


N 

Y(N)  =  / 

-N 


sin46 

e2 


ae 


( see  Table  B-l) . 


Hence,  Eq. 


(B-Q)  may  be  written  as 


wT(o,o) 


(B-10) 


The  procedure  and  theory  discussed  above  was  applied  to 
measuring  the  spectrum  of  Kodak  2479  film  developed  in  D-1Q. 

The  parameters  of  the  system  and  film  were:  t2  =0.31,  f  =  1  m, 
>  =  6328  A,  aperture  D  =  30  mm,  and  slit  length  l  =  1  mm. 

The  nominal  slit  width  d  was  1  a  but  an  experiment*  indicated 
d  =  l.c  U  was  more  likely.  Also  the  diffraction  pattern  of 
the  real  half-wave  plate  was  not  quite  as  given  by  Eq.  (B-2) 
because  there  existed  a  narrow  obscuration  between  the  plates 
whose  effect  was  to  furnish  a  pattern  described  by 

HH(u,v)  =  (D-6  )2  sinc?  ^  (D-6)  sin?  ~  (D+6) 


*  The  ratio  of  the  zero-order  peak-to-null  in  the  absence  of 
film  was  27  db.  This  number  is  a  function  of  the  slit  width, 
which  in  this  case  would  have  to  be  1.5  u.  A  la  slit  would 
furnish  a  30-db  peak-to-null  ratio. 
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TABLE  B-l 


N 

.30000000,-01 

.31000000,-01 

.32000000,-01 

.33000000,-01 

.34000000,-01 

.35000000,-01 

.36000000,-01 

.37000000,-01 

.38000000,-01 

.39000000,-01 

.40000000,-01 

.41000000,-01 

.42000000,-01 

.43000000,-01 

.44000000,-01 

.45000000,-01 

.46000000,-01 

.47000000,-01 

.48000000,-01 

.49000000,-01 

.50000000,-01 

.51000000,-01 

.52000000,-01 

.53000000,-01 

.54000000,-01 

.55000000,-01 

.56000000,-01 

.57000000,-01 

.58000000,-01 

.59000000,-01 

.60000000,-01 

.61000000,-01 

.62000000,-01 

.63000000,-01 

.64000000,-01 

.65000000,-01 

.66000000,-01 

.67000000,-01 

.68000000,-01 

.69000000,-01 

.70000000,-01 


.17993507,-04 

.19853020,-04 

.21836373,-04 

.23947545,-04 

.26190530,-04 

.28569305,-04 

.31087852,-04 

.33750140,-04 

.36560176,-04 

•39521900,-04 

.42639320,-04 

.45916393,-04 

.49357103,-04 

.52965413,-04 

.56745300,-04 

.60700680,-04 

.64835643,-04 

.69154053,-04 

.73659946,-04 

.78357220,-04 

.83249880,-04 

.88341903,-04 

.93637163,-04 

.99139700,-04 

.10485346,-03 

.11078234,-03 

.11693036,-03 

.12330144,-03 

.12989956,-03 

.13672862,-03 

.14379257,-03 

.15109540,-03 

.15864102,-03 

.16643335,-03 

.17447638,-03 

.18277393,-03 

.19133011,-03 

.20014872,-03 

.20923375,-03 

.21858908,-03 

.22821867,-03 


With  D  =  30  nun  and  6  measured  as  0.6  mm,  the  effect  of  the 
obscuration  on  the  null  of  the  diffraction  pattern  was  negligible. 

The  experimental  results  (see  Fig.  B-9)  yielded  a  value 
of  w  (0,0)  which  was  much  larger  than  the  extrapolated  grain 
noise  spectra  of  Sec.  7.3.  The  increase  was  sufficiently  large 
to  cause  serious  doubt  as  to  whether  the  experimental  technique 
did  indeed  furnish  the  DC  value  of  the  grain-noise  spectrum. 

The  results  raised  the  suspicion  that  despite  the  careful 
attempts  to  reduce  the  phase  noise  with  a  liquid  of  matching 
refractive  index,  the  residual  phase  noise  at  the  origin  played 
a  dominant  role.  To  test  this  assertion,  a  sample  of  2479  film 
was  developed  under  identical  circumstances  except  for  being 
bleached* .  and  tested  for  noise  at  DC.  The  grain  noise  of 
bleached  film  is  negligible, and  any  significant  noise  observed 
at  the  origin  must  be  attributed  to  phase  noise.  The  results 
are  shown  in  Fig.  B-10  and  strongly  imply  that  the  grain 
noise  at  the  origin  cannot  be  measured  in  the  presence  of  the 
residual  phase  noise  of  the  film. 

The  results  of  Sec.  7.4  demonstrated  that  the  noise 
associated  with  the  random  phases  of  film  was  a  peaked,  low- 
frequency  process.  The  use  of  a  liquid  gate  is  an  effective 
means  of  eliminating  most  of  the  phase  noise  except  for  vari¬ 
ations  in  the  refractive  index  of  the  gelatin  which,  it  seems,  can¬ 
not  be  disregarded  when  attempting  a  grain-noise  measurement  of  the 
origin.  The  results  of  Figs.  B-9  and  B-10  demonstrated  that 
the  grain  noise  must  indeed  be  very  small, as  witnessed  by  the 
insignificant  change  in  A  under  circumstances  that  either 
include  or  negate  the  granularity.  The  measurement  of  grain 
noise  at  the  origin  requires  a  technique  by  which  the  phase 
noise  is  rendered  insignificant.  Under  such  circumstances. 


*  A  process  that  removes  the  exposed  silver  grains  in  the 
emulsion, thus  in  effect  rendering  the  film  transparent  and 
grain-free . 
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FIG.  B-9  MEASUREMENT  OF  NOISE  AT  ORIGIN  OF  KODAK  2479 

DEVELOPED  IN  D- 19 


RELATIVE  AMPLITUDE  (db) 


I _ 

FIG.  B-l 
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SEE  TEXT. 
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0  MEASUREMENT  OF  FILM  NOISE  AT  ORIGIN  OF  BLEACHED  KODAK  2479 
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the  half-wave  plate  could  be  used  to  advantage  for  the  meas¬ 
urement  of  spectra  in  the  presence  of  a  zero-order  light  that 
is  several  orders  of  magnitude  above  the  spectra  to  be  measured. 
The  excellent  correspondence  between  the  theoretical  and  actual 
performances  of  the  half-wave  plate  establishes  it  as  a  valuable 
tool  for  spectral  analysis . 
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